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PREFACE 


This  report  was  prepared  by  Dr.  Thomas  E.  Owen,  Dr.  Jorge  0.  Parra,  and 
Mr.  James  C.  Biard  of  the  Department  of  Geosciences,  Electronic  System  Divi¬ 
sion,  Southwest  Research  Institute,  under  contract  number  DACA39-86-K-0017 . 
represents  the  first  task  in  investigating  and  applying  high-resolution  shear 
wave  seismic  reflection  exploration  techniques  to  the  problem  of  search  and 
detection  of  subversive  Intrusion  tunnels  of  military  concern.  The  study  was 
done  during  FY  87  under  funding  provided  by  LABCOM  Headquarters  to  the  U.S. 
Array  Engineering  Waterways  Experiment  Station  (WES)  as  a  part  of  the  research 
and  development  effort  to  solve  the  Korean  tunneling  problem.  Dr.  Donald 
Eccleshall  of  the  U.S.  Army  Ballistic  Research  Laboratory,  Army  Material  Com¬ 
mand  (AMC),  is  the  technical  reviewer  of  the  overall  tunnel  detection  program 

This  report  is  presented  in  two  volumes.  Volume  I  contains  the  body  of 
the  report  including  the  technical  results,  discussions,  and  conclusions. 
Volume  II  contains  the  appendices  in  which  the  complete  theoretical  analyses, 
computer  numerical  calculations,  and  computer  programs  are  presented. 

Mr.  Robert  F.  Ballard,  Jr.,  Earthquake  Engineering  and  Geophysics  Divi¬ 
sion  (EEGO),  Geotechnical  Laboratory  (GL),  WES,  was  the  Contracting  Officer's 
Representative  and  Program  Manager  and  Coordinator  for  tunnel  detection 
research.  General  supervision  was  provided  by  Dr.  Arley  G.  Franklin,  Chief, 
EEGD,  and  Dr.  William  F.  Marcuson  III,  Chief,  GL. 

Commander  and  Director  of  WES  during  the  publication  of  this  report  was 
Colonel  Dwayne  G.  Lee,  CE.  Dr.  Robert  W.  Whalin  was  the  Technical  Director. 
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INTRODUCTION  AND  SUMMARY  OF  PROJECT 


A.  Background 

During  the  past  decade,  military  priority  has  been  given  to  the 
search,  detection,  mapping,  and  neutralization  of  subversive  intrusion  tunnels 
constructed  by  North  Korea  in  the  DMZ  and  into  South  Korea.  Three  such  tun¬ 
nels  have  been  found  to  date  and  others  are  expected  to  exist.  The  fact  that 
approximately  eight  years  have  elapsed  since  the  third  and  largest  tunnel  was 
discovered,  reinforced  by  evidence  that  such  tunnel  construction  has  been  con¬ 
tinued  by  North  Korea,  implies  that  additional  tunnels  could  now  be  finished  in 
construction  and  lie  dormant  for  future  subversive  use  by  North  Korea.  There¬ 
fore,  while  current  military  efforts  have  emphasized  the  detection  of  the  tun¬ 
nel  construction  process  using  passive  monitoring  to  detect  and  record  seismic 
waves  generated  by  such  underground  operations,  there  is  now  a  need  for  other 
search  methods  capable  of  detecting  dormant  tunnels. 

Extensive  military  drilling  efforts  are  currently  being  used  to 
investigate  suspected  tunnel  construction  activities  inferred  from  the  infor¬ 
mation  collected  and  analyzed  from  passive  seismic  signal  monitoring.  This 
procedure  has  encouraged  the  use  of  relatively  deep  borehole  drilling,  combined 
with  hole-to-hole  electromagnetic  probing,  as  a  means  of  exploring  for  tunnel 
targets.  As  a  supplemental  method  of  tunnel  exploration  applicable  to  locating 
dormant  tunnels,  surface-based  seismic  reflection  techniques  are  judged  to 
offer  good  potential  for  efficient  operation  without  the  need  for  deep  drill¬ 
ing.  In  this  regard,  high-resolution  seismic  reflection  exploration  techniques 
offer  greater  promise  for  detecting  the  subversive  tunnel  targets  than  any 
other  surface-based  geophysical  technique. 

B.  Summary  of  Project  Efforts 

This  project  has  been  directed  toward  the  study  and  evaluation  of 
shear-wave  reflection  seismic  exploration  techniques  for  detecting  cavities 
and  tunnels.  Shear  waves  offer  certain  specialized  advantages  for  the  tunnel 
detection  application  because  the  polarization  of  shear  waves  may  be  usefully 
aligned  with  the  generally  known  directional  orientation  of  the  North  Korea 
intrusion  tunnels  and,  because  of  their  lower  velocity,  shear  waves  offer 
inherently  greater  resolving  power  in  comparison  with  compressional  waves. 
Horizontally  polarized  shear  waves,  in  particular,  offer  an  added  advantage  in 
that  the  Interference  and  masking  effects  caused  by  the  directly  transmitted 
waves  are  significantly  less  than  those  associated  with  compressional  waves  or 
vertically  polarized  shear  waves. 

Volume  I  -  Study  and  Design  of  Techniques,  as  reported  herein,  pre¬ 
sents  a  comprehensive  study  of  horizontally  polarized  shear-wave  exploration 
concepts  specialized  to  the  military  tunnel  search  and  detection  problem.  The 
primary  efforts  have  been  directed  toward  the  formulation  of  analytical  models 
representative  of  high-resolution  shear-wave  reflections  and  scattering  from 
cylindrical  tunnel  targets  in  dissipative  media  and  their  numerical  evaluation 
to  produce  characteristic  synthetic  seismograms.  Special  attention  was  given 
to  comparing  the  effectiveness  of  surface-coupled  and  borehole-coupled  seismic 
sources  and  detectors  because  of  the  significant  differences  in  the  time  and 
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logistics  required  to  apply  these  two  methods  in  the  field.  Finally,  an  opti¬ 
mum  shear-wave  seismic  reflection  survey  system  and  special-purpose  data  analy¬ 
sis  techniques  are  specified  for  future  development  and  application  to  military 
requirements  for  tunnel  search  and  detection. 

Tractable  analytical  solutions  for  the  scattering  of  plane  and  cylin¬ 
drical  horizontally  polarized  (SH)  waves  from  a  cylindrical  cavity  located  in 
homogeneous  dissipative  media  were  developed  and  computer  model  evaluations 
were  carried  out  for  cases  of  primary  Interest.  The  plane-wave  analysis  was 
useful  in  providing  insights  into  the  Intrinsic  scattering  of  SH  waves  polar¬ 
ized  parallel  to  the  axis  of  the  cylindrical  cavity.  In  particular,  these 
results  present  the  SH-wave  energy  scattering  cross-section  of  a  cylindrical 
tunnel  cavity  as  a  function  of  frequency,  the  angular  amplitude  scattering 
pattern  about  the  cavity  axis,  and  the  steady-state  frequency  response  of 
the  cavity  reflections  as  observed  by  a  seismic  detector  array  in  a  lossless 
homogeneous  whole-space  medium.  The  cylindrical  wave  analysis  has  provided 
synthetic  seismograms  illustrating  reflections  from  the  tunnel  cavity  in  a 
lossless  homogeneous  half-space  and  in  a  two-layer  half-space  in  which  dif¬ 
ferent  dissipative  absorption  conditions  may  be  specified  for  each  layer. 

The  energy  scattering  cross-section  of  the  cylindrical  cavity  for 
incident  plane  SH  waves  is  greater  than  the  projected  area  of  a  unit  length  of 
the  cavity  whenever  the  wavelength  of  the  incident  wave  is  equal  to  or  less 
than  the  circumference  of  the  cylindrical  cavity.  At  longer  wavelengths,  the 
scattered  signal  amplitude  decreases  rapidly  in  proportion  to  the  -2.5  power 
of  frequency.  At  shorter  wavelengths,  the  energy  scattering  cross-section 
becomes  independent  of  frequency  and  approaches  a  value  of  twice  the  projected 
area  of  a  unit  length  of  the  cavity  at  very  high  frequencies.  These  results 
are  significant  in  that  they  indicate  the  most  practical  frequency  range  for 
detecting  scattered  signals  from  a  tunnel  cavity  to  be  the  two-octave  range 
extending  from  a  lower  limit  of  approximately  fj  =  v/uD  to  an  upper  limit  of 
approximately  £2  =  Av/nD.  That  is,  for  incident  SH  waves  having  wavelengths 
greater  than  about  one  circumference  of  the  cylindrical  cavity,  the  scattered 
energy  rapidly  becomes  negligible  whereas,  for  incident  waves  having  wave¬ 
lengths  less  than  about  one  fourth  of  the  cavity  circumference,  the  increase 
in  scattered  energy  becomes  less  significant  as  frequency  increases,  whereas 
the  scattering  and  absorption  of  practical  geological  host  media  become  more 
and  more  significant.  In  the  granite  rock  host  formations  of  Korea,  the 
optimum  SH-wave  frequency  range  for  detecting  the  tunnels  of  interest  is 
approximately  400  Hz  to  1,600  Hz. 

The  synthetic  seismograms  derived  in  this  study  are  most  important 
for  the  lossy  two-layer  half-space  case.  In  this  case,  reflections  from  the 
cylindrical  cavity  form  a  hyperbolic  pattern  characteristic  of  the  localized 
scattering  target  and  are  most  readily  detectable  when  the  surface  layer 
(representing  a  low  shear  strength  surface  soil  layer)  is  either  absent  or 
when  the  seismic  SH-wave  source  and  detectors  are  installed  in  shallow  drill 
holes  that  allow  them  to  be  coupled  into  competent  ground.  Conversely,  when 
the  detectors  are  installed  at  the  surface  of  the  soil  layer,  the  absorptive 
effects  of  this  lossy  layer  and  added  SH-wave  multiple  reflections  in  the 
direct  source-to-detector  waves  impose  excessive  attenuation  of  the  cavity 
reflections  as  well  as  unwanted  masking  effects.  Therefore,  a  principal  con¬ 
clusion  drawn  from  this  analysis  is  that  the  SH-wave  source  and  detectors 
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must  be  operated  In  shallow  boreholes  that  provide  access  to  reasonably  high- 
quality  ground  capable  of  supporting  shear-wave  propagation  in  the  desired 
frequency  range. 

A  field  data  acquisition  system,  together  with  field  procedures  spe¬ 
cialized  to  SH-wave  operation,  has  been  specified  for  future  development  and 
application  in  Korea.  The  source  transducer  recommended  for  this  application 
is  an  assymetrical-force  arc  discharge  device  designed  for  borehole  operation 
and  capable  of  generating  repetitive  SH-wave  pulse  signals  having  spectral 
frequency  content  up  to  about  2,000  Hz.  The  recommended  seismic  detectors 
consist  of  three-component  geophones  configured  to  clamp  in  shallow  boreholes. 
The  source  and  detectors  are  designed  to  be  accurately  oriented  to  optimize 
SH-wave  data  acquisition.  Other  components  of  the  SH-wave  seismic  survey  sys¬ 
tem  include  a  48-channel  digital  data  recording  system  and  a  truck-mountable 
shallow  borehole  drill  to  provide  a  survey  system  capable  of  complete  self- 
contained  field  operations. 

Data  processing  techniques  have  been  outlined  taking  into  account 
the  results  of  the  analytical  studies  and  derived  synthetic  seismograms  as 
well  as  the  practical  aspects  of  high-resolution  seismic  signal  analysis  in 
realistic  geological  formations  when  the  sought-after  target  is  a  localized 
tunnel  cavity.  The  computer  models  developed  as  part  of  this  project  will  be 
valuable  aids  in  further  studies  of  these  specialized  data  processing  tech¬ 
niques.  They  may  also  be  incorporated  into  the  data  reduction  process  to 
provide  a  powerful  means  for  recognizing  and  extracting  predictable  tunnel 
reflection  signatures. 

The  overall  results  of  this  study  and  design  of  SH-wave  seismic 
reflection  exploration  techniques  are  highly  encouraging  with  respect  to  the 
successful  search  and  detection  of  tunnel  targets  typical  of  those  known  to 
exist  in  Korea.  However,  to  be  successful,  the  necessary  features  of:  (1) 
the  unprecedented  high-resolution  SH-wave  data  acquisition  system;  (2)  the 
specialized  data  acquisition  field  procedures  and  data  quality  control;  and 
(3)  the  tunnel  target-specialized  data  processing  and  analysis  techniques  must 
be  recognized  and  carefully  implemented  and  perfected. 


II.  THEORETICAL  ANALYSIS  OF  SHEAR-WAVE  MODELS 


A.  Status  of  Modeling  of  Seismic  Reflections  from  Tunnel  Cavities 

Scattering  of  elastic  waves  by  cylindrical  cavities  has  been  the 
subject  of  numerous  investigations  in  the  last  decade.  Most  of  these  studies 
have  dealt  with  the  case  of  an  infinite  medium  and  a  semi -inf inite  lossless 
half-space  with  regard  to  applications  in  the  field  of  earthquake  engineering. 
The  main  subject  of  interest  has  been  in  the  analysis  of  stress  amplitudes 
along  the  tunnel  circumference  rather  than  to  investigate  the  problem  of 
detectability  of  deep  cavity  targets. 

The  analysis  of  time  harmonic  elastic-wave  propagation  from  a  line 
(or  point)  source  in  a  half-space  or  from  a  cavity  in  an  infinite  elastic 
medium  are  now  a  well-known  part  of  the  literature  on  elastic-wave  propagation. 
(Miklowitz  (1960),  Ewing  et  al.  (1957),  and  Pao  (1973)  give  a  good  collection 
of  references.]  The  more  difficult  problem  of  propagation  of  elastic  waves  in 
a  half-space  containing  a  cylindrical  cavity  has  been  developed  by  Gregory 
(1967).  The  application  of  matched  asymptotic  expansions  to  obtain  the  scat¬ 
tered  field  when  the  wavelength  is  large  compared  with  the  linear  dimensions 
of  the  cavity  has  been  implemented  by  Datta  et  al.  (1978). 

Other  references  to  seismic  reflection  from  tunnel-like  cavities 
that  have  been  helpful  and  relevant  to  the  present  study  include  Lewis  &  Kraft 
(1976),  Lewis  et  al.  (1976),  White  (1965),  Onda  and  Komaki  (1968),  Sato  (1969), 
Kong  (1986),  El-Akily  and  Datta  (1980),  deHoop  (1960),  Apsel  (1979),  Achebach 
(1975),  Akl  and  Lamer  (1970),  Bouchon  and  Aki  (1977),  Garvin  (1956),  Ichikawa 
et  al.  (1975),  Lamb  (1904),  Love  (1944),  Newlands  (1952),  and  Wong  et  al. 
(1974). 


The  contributions  of  the  present  study  have  focused  exclusively  on 
horizontally  polarized  (SH)  shear  waves  with  emphasis  on  the  detection  of 
reflections  from  a  cylindrical  cavity  in  a  dissipative  medium.  The  case  of  a 
two-layer  half-space  composed  of  separately  specified  dissipative  layer  mate¬ 
rials  has  provided  a  simple  but  realistic  model  of  the  typical  field  conditions 
to  be  encountered  in  Korea.  Tunnel  detection  feasibility  and  the  recommended 
field  data  acquisition  techniques  have  been  derived  from  this  model  as  well  as 
from  the  intermediate  model  formulations  utilizing  plane  waves  and  cylindrical 
Waves. 


B.  Rationale  for  Horizontally  Polarized  Shear  Waves 

Commercial  seismic  exploration  methods  rely  primarily  on  the  use  of 
compressional  waves  generated  by  explosive  shots  or  low-frequency  vibrators 
and  are  generally  directed  toward  mapping  relatively  large  and  deep  subsurface 
geological  structures  and  spatial  trends  as  compared  with  resolving  small 
localized  anomalies.  For  this  reason,  conventional  seismic  methods  are  seldom 
required  to  achieve  small  target  detection  or  unusually  high  resolution. 
However,  when  seismic  detection  of  small  anomalies  such  as  man-made  tunnels 
located  at  typically  100  diameters  below  surface  is  considered,  the  seismic 
detection  sensitivity  and  resolution  of  the  method  must  be  increased  substan¬ 
tially  above  that  of  conventional  seismic  systems.  Further,  since  natural 
geologic  anomalies  comparable  in  size  with  the  tunnel  targets  of  interest  may 
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be  present  in  most  metamorphlc  geological  formations,  the  data  processing 
requirements  for  accurate  tunnel  detection  and  Interpretation  are  also  more 
demanding  than  conventional  seismic  exploration  requirements.  Nevertheless, 
conventional  reflection  seismic  methods  are  especially  well  developed  for 
their  industrial  exploration  purposes  and  can  serve  as  the  technical  point  of 
departure  for  developing  the  specialized  methods  needed  for  tunnel  detection. 

Detection  of  man-made  tunnel  cavities  in  the  size  range  of  two  meters 
in  diameter  requires  seismic  signal  wavelengths  an  order  of  magnitude  shorter 
than  those  used  in  conventional  seismic  exploration.  Thus,  the  spectral  con¬ 
tent  of  the  seismic  wave  source  must  extend  an  order  of  magnitude  higher  in 
frequency  than  that  provided  by  conventional  explosive  and  vibrator  sources. 

The  lower  velocity  of  shear  waves  in  comparison  with  compressional  waves  offers 
an  inherent  advantage  in  that  the  high-frequency  content  of  the  source  signal 
need  be  extended  only  by  a  factor  of  about  5  or  6  to  provide  a  resolution  limit 
comparable  with  that  obtained  using  compressional  waves.  The  polarization 
characteristics  of  shear  waves  offer  an  additional  degree  of  freedom  over  com¬ 
pressional  waves  and,  since  tunnel  targets  are  elongated  and  directionally  ori¬ 
ented  anomalies,  this  polarization  parameter  may  be  used  to  good  advantage  in 
the  preferential  detection  of  such  targets. 

Horizontally  polarized  shear  waves  have  several  distinctive  advan¬ 
tages  over  compressional  waves  and  vertically  polarized  shear  waves  with 
respect  to  hlgh-resolutlon  seismic  reflection  exploration  for  tunnels.  These 
advantages  stem  from  the  horizontal  orientation  of  the  shear-wave  particle 
motions  and  their  Inherently  low  propagation  velocity  in  comparison  with  com¬ 
pressional  waves.  These  advantages  are  summarized  as  follows: 

(1 )  Horizontal  Polarization  Orientation 

•  Polarization  alignment  with  the  tunnel  axis  produces  stronger 
reflections  than  compressional  waves  of  the  same  wavelength 
because  the  incident  compressional-wave  energy  is  parti¬ 
tioned,  upon  reflection,  into  components  of  compressional 
waves  and  vertically  polarized  shear  waves,  whereas  incident 
SH  waves  are  reflected  without  a  change  in  wave  type. 

•  Polarization  alignment  of  SH  waves  with  geological  layer 
structures  in  the  tunnel  host  medium  offers  similar  advan¬ 
tages  to  alignment  with  the  tunnel  axis  since  SH-wave  energy 
transmitted  (or  reflected)  from  a  layer  Interface  is  not  par¬ 
titioned  to  other  wave  types  when  the  particle  motions  are  in 
the  plane  of  the  interface. 

(2)  Low-Velocity  Propagation 

•  The  shorter  wavelength  of  shear  waves  in  comparison  with  com¬ 
pressional  waves  provides  greater  detection  sensitivity  and 
spatial  resolution  of  tunnel  targets  than  compressional  waves 
of  comparable  frequency  spectrum.  This  effective  reduction 
in  source  signal  frequency  spectrum  for  SH  waves  offers  the 
prospect  of  less  dissipative  attenuation  along  the  two-way 
signal  propagation  path  and  places  less  demand  on  the  source 
transducer  to  generate  the  high-frequency  signals  necessary 
for  useful  target  detection. 
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(3)  Reduction  of  Coherent  Noise 


•  Compressional-wave  sources  at  or  near  the  surface  generate 
Rayleigh  waves  (surface  waves)  which  often  interfere  with 
body  waves  (compresslonal  or  shear  waves)  reflecting  from 
subsurface  anomalies.  Since  surface  waves  consist  of  verti¬ 
cal  and  radial  shear-particle  motions,  seismic  detectors 
designed  to  detect  compresslonal  waves  also  respond  to  sur¬ 
face  waves.  In  contrast,  seismic  detectors  designed  to 
detect  SH  waves  do  not  respond  to  surface  waves  associated 
with  the  SH-wave  source  transducer  or  other  compresslonal 
wave-noise  sources. 

•  Shear  waves  in  the  frequency  range  appropriate  for  tunnel 
detection  experience  strong  attenuation  effects  in  the  near¬ 
surface  weathered  layer.  Therefore,  direct  excitation  of  the 
detectors  by  the  SH-wave  source  transducer  will  be  advanta¬ 
geously  reduced. 

•  The  dissipative  attenuation  of  SH  waves  in  unconsolidated 
weathered  materials  near  the  surface  imposes  the  need  to 
drill  down  through  such  materials  to  place  the  SH  source  and 
detectors  in  competent  ground.  (Note:  This  requirement 
does  not  necessarily  require  drilling  to  bedrock  depths.) 
Although  this  requirement  places  extra  demands  on  the  field 
data  acquisition  process,  such  drilling  and  competent  trans¬ 
ducer  placements  are  also  demanded  by  high-resolution 
compressional-wave  techniques. 

•  Since  the  polarization  of  the  SH-wave  source  may  not  always 
be  in  alignment  with  the  tunnel  axis  to  give  the  strongest 
SH-wave  reflections,  a  seismic  data  collection  system  which 
records  two  orthogonal  horizontal  components  will  aid  in 
improving  the  reflected  SH-wave  signal-to-noise  ratio  and  may 
yield  useful  information  on  the  target  alignment  direction 
relative  to  the  seismic  survey  line.  (Note:  Misalignments 
between  the  SH-wave  polarization  orientation  and  the  tunnel 
axis  will  result  in  conversion  of  part  of  the  incident  SH- 
wave  energy  to  compresslonal  and  vertically  polarized  shear 
waves.  Since  such  wave-type  conversion  nay  be  unique  to  tun¬ 
nel  cavity  geometries,  a  seismic  data  collection  system  which 
records  three  orthogonal  components  of  particle  motion  is 
recommended  for  use  in  the  field.) 

C.  Organization  of  Theoretical  Model  Analyses 

The  shear-wave  model  studies  presented  in  this  report  are  organized 
in  three  parts.  The  first  part  addresses  the  questions  of  SH-wave  scattering 
from  a  cylindrical  cavity  in  an  unbounded  homogeneous  host  medium.  Ideas  and 
results  from  this  simplified  stage  of  analysis  aided  in  obtaining  a  useful 
perspective  of  the  overall  analysis  problem  and  were  valuable  in  planning  the 
more  extensive  parts  of  the  theoretical  study.  The  second  part  addresses 
SH-wave  scattering  from  a  cylindrical  cavity  in  a  homogeneous  half-space. 
Together,  the  analyses  presented  in  these  two  parts  cover  the  cases  of  a  plane 
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wave  source,  an  infinite  line  source,  and  a  point  source  with  emphasis  on  polar 
amplitude  scattering  patterns,  seismic  energy  reflection  cross-section,  and 
reflection  responses  derived  from  source  and  detector  geometries  associated 
with  reflection  seismic  surveys  for  relatively  deep  tunnel  cavities.  An  exten¬ 
sion  of  this  theoretical  analysis  includes  the  case  of  SH-wave  scattering  from 
a  cylindrical  cavity  in  a  lossy  unbounded  medium. 

The  third  part  of  this  investigation  presents  the  analysis  of  SH-wave 
scattering  from  a  cylindrical  cavity  in  a  lossy  two-layer  half-space.  This 
analysis,  not  reported  in  the  literature  yet,  covers  the  cases  of  an  infinite 
line  source  in  the  bedrock  with  the  detector  located  either  in  the  bedrock  or 
at  the  top  of  the  surface  layer.  To  reduce  the  computations  necessary  to  pro¬ 
duce  synthetic  seismograms  and  as  a  check  of  the  analysis,  the  high-frequency 
approximation  for  the  total  SH-wave  particle  displacement  was  also  derived  and 
evaluated.  In  addition,  the  analytical  solution  of  the  transient  SH-wave  dis¬ 
placement  in  the  presence  of  a  two-layer  half-space  was  developed  using  the 
Cagnlard-deHoop  method  of  inversion.  The  classical  first  motion  approximation 
(Helmberger,  1968)  for  a  disturbance  in  an  elastic  medium  was  also  implemented 
in  this  analytical  solution. 

Theoretical  highlights  and  graphical  displays  of  computer  calcu¬ 
lations  are  presented  in  each  part  to  describe  and  Illustrate  the  derived 
results.  The  complete  theoretical  analyses  together  with  supplemental  data 
displays  and  the  documented  computer  model  programs  are  presented  in  Volume  II 
-  Appendices  of  this  report. 

D.  SH-Wave  Scattering  from  a  Cylindrical  Cavity  in  an  Unbounded 

Homogeneous  Medium 


The  steady-state  solution  for  elastic-wave  scattering  from  an 
infinitely  long  cylindrical  cavity  in  an  unbounded  homogeneous  lossless  solid 
medium  is  developed  for  incident  plane  SH  waves  and  for  incident  spherical  SH 
waves.  The  spherical  SH  waves  are  generated  by  a  point  horizontal  force  within 
the  unbounded  medium.  Formal  analytical  expressions  are  developed  for  the 
elastic-wave  displacement  of  scattered  waves  at  any  point  around  the  cylindri¬ 
cal  cavity  for  incident  plane  waves  oriented  parallel  to  the  axis  of  the  cavity 
or  for  point-source  SH  waves  having  polarization  oriented  parallel  to  the  axis 
of  the  cavity.  In  addition,  theoretical  analvsls  of  the  steady-state  elastic- 
wave  scattering  from  an  infinitely  long  cylindrical  cavity  in  a  homogeneous 
lossless  solid  half-space  is  developed  for  an  SH-wave  point  source  and  for  an 
SH-wave  line  source  at  any  position  within  the  half-space.  Formal  analytical 
expressions  are  derived  for  the  elastic-wave  displacements  to  investigate  the 
Influence  of  the  free-surface  Interface  on  the  scattering  of  SH  waves  from 
cylindrical  tunnel  targets. 

1.  Formulation 


For  elastic  waves  in  an  isotropic  solid  body,  the  vector 
motional  displacement,  ti,  in  the  absence  of  body  forces,  satisfies  the  dis¬ 
placement  equation  of  motion: 


(\  -^1 )  V  V  •  u  +  pV  ^u 


(1) 


where  X  and  ^  are  Lame's  elastic  moduli  and  p  is  the  mass  density  of  the  body. 
In  the  most  general  case,  the  displacement  components  in  Equation  (1)  are  func“ 
tions  of  the  time,  t,  and  the  three  spatial  coordinates  x,  y,  and  z. 

For  a  tunnel  or  cylindrical  cavity  subject  to  a  special  distri¬ 
bution  of  external  forces,  various  approximations  can  be  made  regarding  the 
displacement  components  in  order  to  simplify  the  equation  of  motion  given  by 
Equation  (1).  For  the  long  cylindrical  tunnel  shown  in  Figure  II~1,  the  motion 
can  be  classified  as  either  anti-plane  strain  or  plane  strain  in  accordance 
with  how  the  external  forces  are  applied.  When  a  distributed  force  is  applied 
parallel  to  the  cavity  length  direction  z  such  that  the  dominant  wave  displace¬ 
ment  is  in  that  direction  and  has  constant  magnitude  along  the  z-axis,  the 
motion  is  in  the  category  of  anti-plane  strain.  If  the  force  is  applied  per¬ 
pendicular  to  the  characteristic  length  of  the  cavity  and  is  distributed  uni¬ 
formly  in  the  length  direction,  there  is  no  displacement  along  the  length 
direction.  In  this  case,  the  motion  is  defined  as  plane  strain. 


FIGURE  II-l.  GEOMETRY  OF  CYLINDRICAL  CAVITY  AND  COORDINATE  SYSTEM 

For  the  cylindrical  cavity  shown  in  Figure  II-l,  the  anti-plane 
strain  displacement  is  derived  as 


0 


w(x,y,z)  . 


(2) 


The  assumption  that  the  displacement  is  restricted  to  the  z-axis  reduces  the 
equation  of  motion  given  in  Equation  (1)  to  a  single  scalar^ave  equation 


^i7  2w(x,y,t) 


(3) 


where : 


y 


Under  the  assumption  of  anti-plane  strain,  the  dilatation  7  •  u  Is  zero  (Aki 
and  Richards,  1980)  and  the  resulting  wave  is  a  rotational  shear  wave.  Because 
the  displacement  vector  of  the  wave  is  always  parallel  to  the  z-axis,  which  for 
convenience  is  assumed  to  be  horizontal,  the  waves  are  horizontally  polarized 
shear  (SH)  waves. 

2,  Scattering  of  Plane  SH  Waves  by  a  Cylindrical  Cavity 

The  cylindrical  cavity  illustrated  in  Figure  II-2  has  a  radius, 
a,  and  is  infinite  in  length  in  the  cylindrical  coordinate  system  (r,(J>,z).  The 
lossless  elastic  medium  external  to  the  cylinder  has  a  density,  p,  and  a  shear 


figure  II-2.  GEOMETRY  OF  PLANE  SH-WAVE  INCIDENT  ON  A  CYLINDRICAL  CAVITY 
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modulus,  To  consider  the  scattering  of  shear  waves  by  the  cylindrical 
cavity,  an  incident  plane  SH  wave,  having  its  particle  motion  polarized  paral¬ 
lel  to  the  axis  of  the  cavity,  is  defined  as 

Ujj  =  0,  Uy  =  0,  and  U2  =  w^^^  =  ,  (4) 

where : 

0)  =  angular  frequency; 

k  =  oj/cg  “  wave  vector  amplitude;  and 
Cg  =  SH-wave  velocity. 

That  is,  the  SH  wave  propagates  in  Che  positive  x-direction  with  constant 
velocity,  Cg,  and  angular  frequency,  u).  Upon  impinging  on  Che  surface  of  the 
cavity,  the  incident  wave  is  scattered.  The  scattered  wave  is  also  of  the  SH 
type  and  is  represented  by 

Ux  =  0,  Uy  =  0,  and  Uz  =  w(s)(x,y,t)  .  (5) 

Since  the  geometry  of  the  cavity  is  expressed  in  cylindrical 
coordinates,  all  field  quantities  will  also  be  expressed  in  cylindrical  coor¬ 
dinates.  The  nonvanishing  displacement  and  stress  components  in  cylindrical 
coordinates  are 


®rz  “  I7  “  0,  at  r  =*  a.  (8) 

In  addition,  no  disturbance  should  be  present  at  r  ■  ®  other  chan  the  source 
plane  SH  wave. 
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Each  particle  on  the  surface  of  the  cylinder  acts  as  a  second¬ 
ary  source  generating  scattered  SH  waves  which  are  represented  by  the  unknown 
function  w(s)(x,y,t).  To  determine  this  secondary  displacement,  the  cylindri¬ 
cal  wave  functions 


H(l).(2)(^r) 

n 


cos 

sin 


(n(|> ), 


for  n  *  0, 1 , 2, . . .  , 


(9) 


must  be  considered.  These  functions  satisfy  the  Helmholtz  equation  In  plane 
polar  coordinates  obtained  by  assuming  a  time  dependence  w(x,y,t)  =■ 
w(x,y,w)e~j‘^^;  that  Is: 


dr^  r  5r  f2 


k2) 


w 


0  . 


(10) 


,The  Hankel  functions  of  the  first  kind,  Hjj  (kr),  and  of  the 
second  kind,  Hjj  "^(kr),  are  related  to  the  Bessel  functions  of  the  first  kind, 
Jn(kr),  and  the  second  kind,  Yjj(kr),  by  the  conjugate  relationship 


-  Jn(kr)  ±  jYn(kr) 


(11) 


The  wave  functions  given  In  Equation  (9)  are  obtained  by  the  method  of  separa¬ 
tion  of  variables.  If  w  =  R(r)0(<|»),  the  Helmholtz  condition  of  Equation  (10) 
can  be  separated  Into  two  ordinary  differential  equations 


and 


2  R  ^ 

dr2  dr 


^  +  (k2r2  -  n2) 


^  +  n2e  =  0  , 

d4i2 


(12) 


where  n,  an  Integer  In  this  case.  Is  the  separation  constant.  The  0  equation 

Is  harmonic,  giving  rise  to  harmonic  solutions  In  sin  n4i  and  cos  nA.  and  the  R 

equation  Is  Bessel's  equation  of  order  n  for  which  H^^)(kr)  and  H(2)(kr)  are 
•1  ^ .  n  n 

solutions. 


When  the  wave  functions  given  by  Equation  (9)  are  combined  with 
the  time  factor,  they  represent  cylindrical  waves  generated  by  a  line 

source  along  the  z-axls.  From  the  asymptotic  behavior  of  H^^^»^2)(kr) 


„(1).(2) 


(kr) 


-  (l+2n)Ti/41  (13) 
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the  function  H^^^(kr)  cos  represents  a  diverging  or  outgoing  cylindri¬ 

cal  wave;  that  is,  the  waves  generated  by  the  line  source  propagate  away  from 
the  origin.  Similarly,  the  function  cos  represents  a  con¬ 

verging  or  incoming  cylindrical  wave.  The  boundary  condition  at  r  =  ®  in  this 
problem  eliminates  the  latter  from  the  scattered  waves,  hence  only  the  function 

H^^)(kr)  is  valid  in  the  solution  for  the  scattered  waves, 
n 

Thus,  the  scattered-wave  function,  w(s)^  may  be  expressed  as 


w^®)(r,$,t)  =  ^  (Itr)  cos  ; 

n=0 


where  the  are  the  scattered-wave  coefficients  which  represent  the  axially 
polarized  shear  component  of  the  scattered-wave  field  in  the  host  medium.  The 
total  SH-wave  displacement  at  the  measurement  point,  P,  is  given  by  the  super¬ 
position  of  the  incident  and  scattered  functions. 


Uj  =  w(r,iti,t)  =  w*'^^  +  w(s)  . 


To  determine  the  unknown  coefficients,  from  the  boundary  conditions  at  the 

surface  of  the  cylindrical  tunnel,  the  Incident  plane  wave  function  given  in 
Equation  (4)  is  written  in  polar  coordinates  as 


y(i)  s  '5’  , 


which  may  be  expanded  in  complex  Fourier  series  is 


gjkr  cos  (|»  =.  ^  Cn(r)eJn^» 


in  which 


Cjj(r)  =  ■=—  /  eJ^^  cos  41  d<{i  , 


1 

__  J  cos  (J)  d()»  . 


The  integral  definition  of  the  Bessel  function 


4 

2  2n 


Jn(kr) 


J  eJ^*^  cos  ni(»  d^i  , 


n 


may  be  compared  with  C^(r)  to  yield 


and 


where : 


Cn(r)  -  , 


gjkr  cos  it*  -  ^  Cnj^JnCkr)  cos  m})  , 

n*0 


e 


n 


1. 

2. 


for  n  *  0 
f or  n  >  1 . 


(17) 


Thus,  the  incident  plane  SH  wave  is  expressed  in  polar  coordinates  as 


w(l-^  ■  ^  £nJ”'^n(l^^^  *^0®  n(}ie“J‘^^  . 
n*0 


Equation  (17)  may  be  regarded  as  a  superposition  of  an  infinite  number  of 
cylindrical  shear  waves,  resulting  in  a  pattern  of  wave  propagation  along 
the  x-axls. 

Omitting  the  time  factor,  the  total  wave  displacement  at 
P,  as  given  by  Equation  (15)  is,  therefore,  represented  by  the  infinite  series, 


cos  ni|>  .  (18) 

n=0 


At  the  surface  of  the  cylindrical  cavity,  the  boundary 
in  Equation  (8)  yields  the  coefficients,  as 


An  = 


^nJ" 


’^'n(ka)  . 

Hn(l)(ka) 


Using  the  Bessel  identities 


for  n  =  0,1,2, 


ka.Tn(ka)  =  nJn(ka)  -  kaJ^+lCka)  , 


condition  given  earlier 


U9) 


kaH'(ka)  -  nH^^^ka)  -  kaH^ii(ka)  , 


and 


Equation  (19)  may  be  reduced  to 


nJn(ka)  -  kaJn+l(ka) 

^  nH(*^)(ka)  -  kaH(|^?(ka) 

n  n+i 


(20) 


Thus,  the  SH-wave  displacement,  w,  at  a  measurement  point,  P,  near  the  cavity 
as  given  by  the  Fourier-Bessel  sum  as 


w  =  I  [enJ"'^n(><^r)  +Gj^(k,a,r)]  cos  nb  ,  (21) 

n=0 


where : 

Gjj(k,a,r)  =  Ajj(ka)jjH(  ^  ^(kr )  . 


To  evaluate  the  amplitude  and  phase  of  the  normalized  SH-wave 
displacement  for  the  vector-wave  amplitude  at  k  =  0,  it  is  convenient  to 
express  the  general  solution  of  the  scattered  SH-wave  displacement  in  terms 
of  small  argument  Bessel  functions.  Thus,  the  scattered-wave  coefficients, 
A^,  given  in  Equation  (19)  for  small  arguments,  (i.e.,  as  ka -»•  0),  are 


An(ka) 
ka  -►  0 


for  n  =  0, 


for  n  >  1 , 


(22) 


and  for  k  *  0,  the  function  =  0  for  n  =  0,1,...  ®. 

A  similar  analysis  may  be  performed  for  the  function  Gjj(k,a,r)  = 
Ajj(ka)H^^^^  (kr)  as  ka  0.  Employing  the  asymtotic  expression  for  the  Hankel 
function  for  small  arguments,  i.e.. 


Hn(kr)  “ 
ka  0 


for  n  =  0 


(n  -  1)1  1  . 

«  (kr/2)" 


for  n  >  1 , 


yields 


Gn(k,a, 


2  (|^)  [log  +  log  (|^)];  for  n  =  0, 

(^)"  “TTTn  ;  for  n  >  1, 

oHt) 


(23) 


(24) 
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and,  for  k  •  0,  the  function  Gjj(k,a,r)  »  0  for  n  ■  0,1,...  « .  The  results  for  | 
k  *  0  in  Equations  (22)  and  (24)  indicate  that  static  non-oscillatory  source  | 
forces  do  not  propagate  as  elastic  waves,  a  fact  that  serves  as  a  check  on  the  t 
derived  analytical  expressions.  ^ 

The  final  expression  for  the  SH-wave  displacement  at  any  mea-  a 
surement  position,  P,  in  the  vicinity  of  the  cavity  in  terms  of  the  function  i 
Gjj  is  given  by  i 


w  ■  I  Gjj(k,a,r)  cos  n^ ,  for  r  >  a 

n=0 

whe  re : 


0;  for  k  =  0 


Gjj(k,a,r) 


n  ~  ^^‘^n-tl  ^kr ) 

nH(^)(ka)  -  kaH(0(ka) 
n  n+i 


for  k  >  0. 


The  total  SH-wave  displacement  in  Equation  (25)  may  be  expressed  in  real  and 
imaginary  components  as 


w  «  gSH  ♦  jhj 


where : 


g-„  -  Real  (w), 
hsH  “ 


Then, 


in  which 


<l)*(kr)  *  tan~khgj^/gg^). 


A  normalized  amplitude,  Ah(kr),  and  relative  phase,  il>^^(kr),  may  be  defined 
by  combining  Equation  (25)  and  the  Incident  SH-wave  displacement  given  in 
Equation  (16).  Thus, 


Ah(kr)  =  (h|jj  + 


4.j^(kr)  -  tan-1  (hgj^/g  )  -  kx 


where  is  in  radians.  Equations  (25)  and  (26)  were  used  to  develop  the 

displacement  amplitude  programs  dispa,  dispb,  dispcd,  and  dispe,  as  well  as  the 
subroutine  plane-scatter  presented  in  Appendix  H. 

3.  Scattered  Energy 

The  rate  of  energy  flow  across  a  surface,  S,  having  a  unit 
normal  vector,  n  ,  is  given  as  (Pao  and  Mow,  1973) 


(27) 


in  which  E  is  the  total  energy  (strain  energy  and  kinetic  energy),  n  •  a  is 
the  stress  vector  at  the  surface  element,  dS,  and  u  is  the  particle  velocity. 
The  integrand  of  Equation  (27)  is  the  actual  rate  of  work  done  by  the  surface 
traction  per  unit  area  and,  by  the  law  of  conservation  of  energy,  this  rate  of 
work  equals  the  flow  of  energy  across  a  unit  area  per  unit  time. 


In  tae  case  of  steady-state  waves,  both  the  stresses  and  dis¬ 
placements,  and  Uj ,  are  harmonic  functions  of  time  which  can  be  written  as 


Oij(xk,t)  -  aij(xit)e“>'^ 


(28) 


and  Uij(xk,t)  »  uj (x}t)e"J‘^*^ 

where  and  uj  are  complex  functions.  Since  Real  (u)  =  (u  +  u*)/2  where  an 
asterisk  indicates  the  complex  conjugate  of  the  corresponding  quantity,  then 
the  energy  may  be  expressed  as 

E  *  -1/4  //gni(crijUj  +  cfijuj  +  o^juj  +aijUj)dS  .  (29) 

Thus,  the  time  average  of  energy  flow  is  given  by 

Ave(E)  »  Y /q  *  - -I-  ju)  //5ni(o£jU  -aijUj)dS  .  (30) 


For  an  Incident  plane  SH  wave,  the  time  average  energy  flux  per  unit  area  is 

Ave(e)  =  — Y -  “  "  ^  jwWxz''  "  *  (31) 

where: 


®xz  “  pke-}  (^x"^  ^2^ 

-* 

‘^xz  ”  pke-j(kx+rt/2) 
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w  “  jikeJ^x 
w*  ••  |ik.e“J^x  . 


After  substitution,  the  time  average  energy  carried  by  the  incident  wave  per 
unit  area  is 


Ave (e^^)  1 


(32) 


The  energy  flux  across  a  cylindrical  surface  r  =  R  (per  unit  length  of  the 
cylinder)  is 


Ave(E)  *  -  X  Joj/  [OfzW  -  Orzwlr=R 
^  o 

Substituting  the  scattered-wave  displacement,  w(®\  and  stress,  ii'to 

Equation  (33)  leads  to 


Ave(E)  (f-)AnA*kR{Hn(kR)H*(kR)  -  H^*(kR)Hn(kR))  , 

n*o  “ 


where  use  has  been  made  of  the  orthogonality  condition 


j  cos  m$  cos  ni^d(>=6uni»^^™^®* 
o 


As  kR  -»■  ®  , 


Hn(kR)  > 


-  (n+l/2)’t/2] 


and,  hence,  H^(kR)  -►  jHj^(kR). 

Thus,  at  a  large  distance,  R, 


(34) 


CO  ^ 

Avele(s)l  ^  vikcs[  ZA^Aq  +  I  A^A^l 

11-1 


(35) 
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A  natural  reference  against  which  to  compare  the  total  scattered  energy  is  the 
time  average  of  the  energy  flux  per  unit  area  of  the  incident  wave  given  by 
Equation  (32).  Thus,  defining  the  ratio  of  these  two  energies  by  Q  yields 


Ave [e(s  ^ 1 
Ave (e^^^ 1 


^  [2AqAq  +  ^  ^^n^  • 

n=l 


(36) 


Since  Q  has  the  dimensions  of  area  per  unit  length,  the  quantity  Ave[E(s)]  nay 
be  referred  as  the  energy  scattering  cross-section  of  the  cylindrical  cavity. 
Thus,  the  normalized  cross-section  (or  the  dimensionless  normalized  cross- 
section)  is  given  by 


where : 


I“nlAn(ka)| 

n=l 

I  ;  for  n  “  0 

«n  “  { 

1/2  ;  for  n  >  I  , 


(37) 


and  the  scattered-wave  coefficients  which  represent  the  axially  polarized 
shear  component  of  the  scattered  wave  field  in  the  host  medium  is  given  by 
Equation  (20). 

Equation  (37)  may  be  used  to  calculate  the  normalized  plane  SH- 
wave  scattering  cross-section  as  a  function  of  the  variable  ka  to  investigate 
the  cavity  scattering  response  at  large  distance  and  as  a  function  of  fre¬ 
quency.  To  calculate  the  angular  distribution  of  the  scattered  energy  at  dif¬ 
ferent  wavelengths,  the  energy  flux  per  unit  area  given  by  the  integrand  of 
Equation  (33)  may  be  used.  That  is. 


Y(^,r)  I  I  (AnA„kHn(kr)H„(kr)  -  A^A^kHn  (kr)H„(kr)]  ,  (39) 

n*0  m*0 
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where  the  scattered-wave  coefficients,  and  A^,  are  obtained  from  Equation 
(20)  and 


krHn(kr)  -  nHn(kr)  -  krH^+iCkr) 


krH^*(kr) 


nHn(kr)  -  krHjj+i(kr)  , 


in  which  Hjj(kr)  = 


Equation  (39)  may  be  used  to  calculate  the  normalized  plane  SH-wave  scattering 
cross-section  as  a  function  of  the  variable,  ka,  to  Investigate  the  change  of 
direction  of  the  scattered  energy  about  the  cavity. 

A.  Scattering  of  Spherical  SH  Waves  by  a  Cylindrical  Cavity 

The  geometry  of  the  point  source  SH-wave  scattering  problem  is 
Illustrated  in  Figure  II-3,  An  infinitely  long  cylindrical  cavity  of  radius, 
a,  occupies  the  space  r  <  a  and  -<»  <  z  <  ®  in  the  cylindrical  coordinate  system 
(r,(^,z).  The  region  r  >  a  external  to  the  cavity  has  a  density,  p,  and  a  shear 
modulus,  p.  The  SH-wave  point  source  in  the  elastic  medium  acts  along  the 
z-axis  (parallel  to  the  cylinder)  and  consists  of  a  body  force,  f(t  -  R/cg), 
at  the  position  S(-Xg,yg,Z3).  The  axis  of  the  cavity  contains  the  origin  of 
the  cartesian  coordinate  system  x,  y,  and  z. 

For  the  conditions  stated,  the  SH-wave  displacement,  u^^^, 
excited  by  the  point  source  in  a  lossless  unbounded  medium  is  a  solution  of 
the  Inhomogeneous  wave  equation  and  is  given  by  deHoop  (1960)  and  White  (1965). 


u(i) 


f(t-R/c3) 

AxpCg^R 


(AO) 


where : 


and 


Cg  *  SH-wave  velocity 


R  -  ((x-xg)^  +  (y-yg)^  +  (Z-Zg)^!^^^ 


i«  the  distance  from  the  SH-wave  source  to  the  measurement  point,  P,  at  coordi¬ 
nates  X,  y,  and  z. 
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FIGURE  II-3.  (EOMETRY  OF  SPHERICAL  SH-WAVE  INCIDENT  ON  A  CYLINDRICAL  CAVITY 


In  practical  seismic  exploration,  the  source  usually  generates 
short  pulses  and  the  objective  is  to  obtain  the  transient  solution  of  SH-wave 
scattering  from  buried  tunnel  targets.  Since  in  most  cases  such  transient 
solutions  are  difficult  to  obtain  directly,  the  steady-state  solution  of  the 
scattering  of  transient  spherical  SH  waves  may  be  synthesized  by  using  the 
Fourier  integral  theorem.  Thus,  the  Incident  SH-wave  displacement  may  be 
expressed  as 


u(^^(R,t)  “2^/  u^^HR,u))e*J'^^dw  (41) 

whe  re : 

u^^^(R,oj)  ■  /  u^^^(R,t)e"*’-}‘^^dt  (42) 


provided  that  /  | u^^ ^(R, t ) |dt  exists. 


The  quantity  u^^^  is  the  Fourier  transform  of  u^i).  If  Equation  (40)  is  sub¬ 
stituted  into  Equation  (42),  the  SH-wave  displacement  excited  by  the  point 
source  is  obtained  in  the  frequency  domain  as 


u(l)  -  ej-— 

4npCg  ^ 


(43) 


where  F(u))  is  the  body  force  in  the  frequency  domain. 

Introducing  the  known  integral  representation. 


eJkR 

R 


dp 


(44) 


21 


where: 


k-rgi  .  [(x-Xs)]2  +  (y-ys)2]^^^  . 


kr  -  (k2-p2)i/2  . 


The  Hankel  function  In  Equation  (44)  may  be  expanded  In  terms  of  cylindrical 
functions  using  the  addition  theorem  for  Bessel  functions, 


Ho(^^(krl*^-%l) 


I  H^^^(rskr)Jn(rki.)e^"^‘^~^®^;  for  r  <  rg 
n=-® 

I  l4^\rkr)Jn(r8kj.)e^”^^~^®\  for  r  >  rg 
n— * 


(45) 


Since  the  point  source  is  outside  the  cylindrical  cavity,  then  the  incident 
SH-wave  displacement  Is  expressed  In  terms  of  cylindrical  functions  as: 


u(i).ilM^j"  I  Hn(l)(rgkr)Jn(rkr)eJn(^“^s>dPe^^^*"^s^  . 

^Cg  "  ^  n-« 

(46) 

«0  «0 

/  I  cosP(z-Zs)dp  . 

^^Cg  0  n— » 


Introducing  the  Integral  summation  operator  defined  by 


r  !•••]  “  /  I  (•••)  cos  p(z-Zg)dp  ,  (47) 

0  n=-<» 


the  operator  T  Is  recognized  as  a  Fourier  Integral  In  the  z-directlon  along 
the  axis  of  the  cylinder  combined  with  a  Fourier  series  In  the  4)  direction. 
Thus,  the  incident  spherical  SH-wave  displacement  expressed  in  a  cylindrical 
mode  expansion  is  given  by 
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C(u)Hn(rskr)Jn(*^rr) 


for  n  =  0,±1,±2,...,  ±  ®  . 


(48) 


2 

where  C(a))  =  jF(tj  )/4iipCg  . 

Since  the  scattered  spherical  SH-wave  displacement  is  a  solution  of  the  Helm¬ 
holtz  equation,  this  displacement  may  be  expressed  in  terms  of  cylindrical 
wave  functions  referred  to  the  cylindrical  cavity  axis  as 


=  BuHn^^Ckj-r)  ,  for  n  =  0,±1 . ±  ®  ,  (49) 

to  represent  outward-traveling  waves.  Since  the  total  SH-vave  displacement 
is  considered  as  the  sum  of  the  incident  and  scattered  waves,  then 


Un  =  C(a))Hn(rski.)V‘«-rr)  +  B^H^^^krt)  . 


(50) 


At  the  surface  of  the  cylindrical  cavity,  r  =  a,  the  stress-free  boundary 
condition 


CT 


rz 


=  0 


leads  to 


B 


n 


C(u))Hn(rgkj.) 


Jn(akr) 

Hj^(akr) 


for  n  =  0,±1,...,  ±  ® 


(51) 


Thus,  substituting  the  scattered-wave  coefficients  given  by  Equation  (49)  into 
Equation  (46)  and  applying  the  integral-summation  operator  defined  by  Equation 
(47),  the  steady-state  solution  for  the  total  SH-vave  displacement  at  the 
measurement  point,  P,  is 


ut(a,)--i^^/  I  eJn<^^s){H„(rgkr)Jn(krr) 
0  n-« 


+  Gn(r,u.,kr)Hn(krc)}  cos  8 (z-Zg )d8  . 


(52) 


where : 


>  k|- ) 


Jn(kra) 

Hn(kra) 


(53) 
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and  the  corresponding  transient  solution  becomes 


J  CD 

Ut(t,r)  “  /  [e“J‘^tF(u))/  D(a),p)  cos  p  (z-Zg  )dp  ]du) 

8ii^c^  0 


(5A) 


where: 


D(a),p)  -  I  Gn(r,a),kr)Hn(krr)e"J('^^s)  . 
n=->» 


5.  Scattering  of  Cylindrical  SH  Waves  by  a  Cylindrical  Cavity 

10  solve  the  steady-state  solution  for  an  SH-wave  line  source, 
the  point  6  rce  representation  given  in  Equation  (40)  may  be  integrated  along 
the  z-axis.  That  is, 


,(1) ,  ^  ^ 

4llC^  -a> 

8 


Next,  using  the  Integral  solution 


CD 

/ 


iijHj^^k|r-rs|)  . 


the  incident  cylindrical  SH-wave  displacement  becomes 


jF(u)) 

Apc^ 

8 


)(k| r-rgl ) 


which  may  be  expressed  in  terms  of  cylindrical  functions  as 


jF((D) 

Apcg^ 


I  Hi*\rgk)Jn(rk)eJ"^^^s)  . 
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HHiggiigi 


As  in  the  SH-wave  point  source  solution  previously  developed,  the  scattered 
cylindrical  SH-wave  displacement  can  be  expressed  in  terms  of  cylindrical-wave 
functions  as 


a?  =  I  C„H„<l)(,k)eJ"<*-*s)  , 


to  roprssent  outward*“tr3V6ling  waves*  Thus^  the  total  SH-wave  displacement  at 
position,  P,  is 


n=- 


Applying  the  stress-free  boundary  condition  at  the  surface  of  the  cylindrical 
cavity  leads  to  the  solution  of  the  scattering  coefficients  for  a  line  source 
of  SH  waves 


c„  -  -  J  k,  . 


4p  ^ 
^  s 


Hn(ak) 


Thus,  the  steady-state  solution  for  cylindrical  SH  waves  scattered  from  a 
cylindrical  cavity  in  an  unbounded  elastic  medium  is 


I  e"J^<’'^s)(H„(rsk)Jn(rk)  +  QnHn(rk)]}  , 

V  Cg  n»-« 


whe  re : 


Jn(ka) 

Qn  -  -;rVrTHn(krg)  , 
Hn(ka) 


and  the  corresponding  transient  solution  is 


Uf(t,r)  * — J j  /  e'J^^^FCcj )N(w )dj  , 
4iipc^  -«o 
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where: 


N(u))  -  I  Hn(kr8){Jn(rk)  +  QnHn(rk) 


6.  Attenuation  and  Dispersion  Along  the  Ray  Path 

The  steady-state  solution  reported  for  elastic-wave  scattering 
from  an  Infinitely  long  cylindrical  cavity  in  an  unbounded  homogeneous  solid 
medium  may  be  extended  to  Include  attenuation  and  dispersion.  The  total  SH- 
wave  displacement  expressed  in  terms  of  a  complex  propagation  constant,  kg,  is 
given  by 


u^(a) ) 


F(a)) 

4itp 


{  jitHoCkg 


r-rs! )  ~  I 

n=-a> 


Jn' (kga) 

Hn’(ksa) 


X  Hn(rsks)Hn(r8ks)eJ"^^'^s^}  (57) 


where: 


a(oj)  ■  the  absorption  coefficient;  and 
c((»))  ••  the  phase  velocity. 

To  calculate  synthetic  seismograms,  the  attenuation  and  disper¬ 
sion  model  developed  by  Futterman  (1962)  was  implemented  in  the  general  steady- 
state  solution  of  the  total  SH-wave  displacement  stated  in  Equation  (57).  The 
corresponding  time  derivative  for  the  total  SH-wave  displacement  was  obtained 
by  the  Fourier  transform  method.  Futterman  states  that  if  the  medium  behaves 
in  a  linear  manner,  then,  as  a  consequence  of  the  principle  of  causality,  the 
presence  of  attenuation  also  requires  the  presence  of  dispersion;  that  is,  if 
the  imaginary  part,  a(u)),  of  the  complex  propagation  constant,  kg,  is  known, 
then  the  real  part,  u)/c(u),  can  be  determined.  In  particular,  if 


(58a) 

then, 

C(U))  -  Co[l  log  (58b) 
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rnmm 


where : 


Cq  =  the  phase  velocity  at  cJq;  and 

Qq  =  value  of  the  quality  factor,  Q,  at  Uq* 

The  Futterman  model  has  wide  applicability  and  covers  a  broad 
frequency  range.  It  is  limited  at  low  frequencies  by  a  finite  arbitrarily 
small  but  nonzero  cut~off  frequency,  u)q,  below  which  there  is  assumed  to  be  no 
attenuation.  That  is. 


kc  = 


oj/cq 

;  for  0)  <  (Jq 

(59a) 

cM  * 

;  for  0)  >  (jJq 

(59b) 

i)  and  (58b) 

into  Equations  (59a)  and  (59b), 

respec- 

tively,  the  parameter,  akg,  for  the  Futterman  model  is 


akg  “ 


ka 


;  for  ka  <  k^a 


^  ^  ^  ®  ;  for  ka  >  k^a  .  (60) 


2Qc 


koaJ 


The  numerical  solution  of  Equation  (60)  is  computed  by  the  complex  function  zk 
program  presented  in  Appendix  H. 

Figure  II-3  shows  the  geometry  used  in  analyzing  SH-wave  scat¬ 
tering  from  a  cylindrical  cavity  in  an  unbounded  lossy  medium  for  incident 
cylindrical  SH  waves.  The  total  SH-wave  displacement  at  measurement  point, 
P(x,y),  consists  of  the  sum  of  the  incident  wave  and  the  wave  scattered  from 
the  cylindrical  cavity.  The  total  SH-wave  displacement  is  expressed  by  Equa¬ 
tion  (57)  presented  earlier,  and  its  corresponding  time  derivative  is 


u’(t)  “  /  F(w  )T(o}  )u)e~J'^*-da)  , 


(61) 


where : 


T(u) 


-“S’  {jtiHo(ks 
4itp  c^ 


r-^sl  )  " 

n=-<»> 


Jn'(ksa) 
^n  * (^s® ) 


X  H„(ksrs)H„(kgr)e-^"^^^s>} 


*^3  "  7^  +  ja  (u) )  ;  and 
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F(u))  “  a  body  force  acting  parallel  to  the  z-axis  at  the  position 
(xs.yg)- 

The  body  force,  F(a)),  is  obtained  from  the  frequency  transform 


F(u))  -  /  f(t)ej‘^tdt 
0 


in  which  f(t)  is  the  applied  SH-wave  pulse  signal  defined,  for  the  present 
analysis,  by 


f(t)  =  pte-?t  sin  (i)QtU(t )  , 


where  U(t)  is  the  unit  step  function.  The  frequency  spectrum  of  f(t)  is  given 
by 


2pUo(p-ju)) 


[  +  Wo^]" 


where : 


P  =0,9  Wq/h  ’,  and 

cJq  =  the  peak  frequency  in  the  pulse  spectrum. 

Equation  (63)  was  used  to  compute  the  real  function  pulse  (t)  program  pre¬ 
sented  in  Appendix  H. 

E .  SH-Wave  Scattering  from  a  Cylindrical  Cavity  in  a  Homogeneous 

Half-Space 


1 •  Scattering  of  Spherical  SH  Waves  by  a  Cylindrical  Cavity 

Figure  II-4  Illustrates  a  semi -inf inite  homogeneous  isotropic 
elastic  half-space  having  a  shear  modulus,  p,  and  a  density,  p.  A  cylindri¬ 
cal  cavity  of  radius,  a,  is  located  in  the  lossless  half-space  at  a  depth,  h, 
from  the  free  surface.  An  SH-wave  point  force  in  the  elastic  medium  acts 
along  the  z-axis  (shear  polarization  parallel  to  the  cylinder)  and  represents 
a  body  force,  f(t  -  R/cg),  at  the  position  S(-Xs  ,ys ,  Zg ) .  The  axis  t|je 
cylindrical  cavity  contains  the  origin  of  the  x,y,z  cartesian  coordinate 
system. 


FREE  SURFACE 


FIGURE  II-4.  CROSS-SECTION  OF  A  CYLINDRICAL  CAVITY  IiN  A 
HOMOGENEOUS  HALF-SPACE 


The  displacement  field  of  the  SH-wave  point  force  satisfies  the 
inhomogeneous  wave  equation  and,  in  the  frequency  domain,  is  represented  in 
terms  of  cylindrical  functions.  For  the  region  y  <  h: 


(j(i)  .  J  ^  cog  p  (z-z^. 


4jipc‘  0 


)dp 


(65) 


where : 


r,  -  Ul  * 

r  =  (x2  +  y2)^/2. 


"  tan'  ■kx/y); 


't>s 


tan'^-Xg/yg); 


and 


(k2  -  p2) 


1/2 


The  analysis  of  scattering  of  spherical  SH  waves  by  a  cylin¬ 
drical  cavity  in  a  homogeneous  half-space  may  be  deduced  in  terms  of  the 
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corresponding  whole-space  solution  by  using  image  theory.  In  fact,  one  of  the 
boundary  conditions  for  this  half-space  problem  states  that  the  traction  is 
zero  at  the  surface  Interface.  This  effect  may  be  achieved  in  a  whole  space 
by  matching  each  source  in  the  lower  half-space  plane  (y  <  h)  with  a  source  of 
equal  magnitude  and  opposite  polarity  placed  at  the  source  image  position  in 
the  upper  half-space  (y  >  h).  The  SH-wave  displacement  distribution  in  each 
half  of  this  composite  whole-space  arrangement  will  then  be  the  same  as  the 
sought-after  displacement  distribution  in  the  half-space  problem. 

From  a  physical  point  of  view,  the  spherical  SH  waves  produced 
by  the  point  source  in  the  lossless  half-space  medium  are  reflected  from  the 
plane-surface  (y  “  h),  and  scattered  and  diffracted  by  the  surface  of  the 
cylindrical  cavity,  r  ”  a.  Therefore,  the  resultant  total  SH-wave  displace¬ 
ment,  u*^,  is  a  superposition  of  the  incident,  reflected,  scattered,  and  dif¬ 
fracted  waves  and  must  satisfy  the  traction-free  boundary  conditions: 


At  the  surface  of  the  earth  (y  =  h): 


'yz 


ay 


At  the  surface  of  the  cavity  (r  “  a): 


a 


rz 


au^ 

aT" 


ss 


0 


(66) 


To  obtain  the  solution  for  the  scattered  SH  waves  by  the  cylin¬ 
drical  cavity,  two  cylindrical  cavities  are  considered  to  exist  in  the  image- 
based  whole-space  geometry  shown  in  Figure  II-5.  The  cavity  of  radius,  a,  in 
the  lower  half-space  (y  <  h)  with  axis  at  (0,0)  occupies  the  space  r  <  a  and 
<  z  <  “  in  the  cylindrical  coordinate  system  (r,4),z).  The  cylindrical 
cavity  of  radius,  a,  in  the  upper  half-space  (y  >  h)  with  axis  at  (0,2h)  occu¬ 
pies  the  space  r^  <  a  and  <  zj,  <  »  in  the  cylindrical  coordinate  system 
(rh.<>h»zh)* 


30 


To  represent  outward'traveling  waves,  the  scattered  SH-wave 
displacements  must  be  of  the  form 


*  /  I  AnH(l)(rkr)eJn^  cos  P (z-Zo)dp 


-  /  I  p  (z-Zo)dp 

0  n**-® 

in  which 

rh  *  [Ah^+r^“4rh  cos  • 


In  the  whole  space,  the  derivatives  with  respect  to  y  of  the 
and  must  vanish  at  y  *  h  because  of  the  symmetry  of 

their  representation.  ft\us,  in  the  half -space,  y  <  h,  the  sum  of  the  incident 
SH-wave  displacements,  u^^)  +  and  the  scattered  SH-wave  displacements, 

u^®)  +  satisfies  the  traction-free  condition  at  y  ■  h.  The  resulting 

total  SH-wave  displacement  in  the  elastic  half-space,  u^,  is 


ut  »  u(i)  +  ^  u(s)  +  u^®  ^ 


and  must  satisfy  the  traction-free  condition  at  the  surface  of  the  cavity. 


Using  the  addition  theorem  for  the  Hankel  function 


HO)(rhkr)eJ"^h  ,  ^  Jm(krr)H(l)(2krh)eJn>l'  , 


the  total  SH-wave  displacement  in  the  wave  number  domain  may  be  written  as 


H  ^  {Hn(rskr)e"^"*^®  +  H^(r®kr)e"^'^^h}  J„(rkr  )ejn4. 

P  4npc2  h 


+  I  AnHn(rkr)eJn4*  +  I  I  VmCk^r )Hn+„(2krh)ej'^  . 
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The  use  of  the  traction-free  boundary  condition 


[ 


d  r 


] 

r=a 


0 


(73) 


leads  to  an  infinite  set  of  algebraic  linear  equations  for  the  determination 
of  the  unknown  wave  coefficients,  Upon  normalizing  Equation  (72)  by  the 

factor  ,  this  linear  system  of  equations  is  given  by 

4iip 


■X 


Jn(kra) 

Hi(kra) 


Hn+a,(2krh)} 


=  ^77^  {Hn(krrh)e“^"^^  +  Hn(krrs )e-jn4's}  ; 
H'(kj.a) 


for  n  =  0,±1,±2,...  i  <» , 


(74) 


where : 


and 


"m _  . 

iF(u))  ’ 

4itpc^ 


6 


nm 


1; 

0; 


for  n  *  m 
for  n  ^  m. 


Finally,  the  steady-state  solution  for  the  total  spherical 
SH-wave  displacement  at  the  measurement  point,  P,  expressed  in  terms  of  the 
scattered-wave  coefficients,  a^,,  is  given  by 


u*^(u) ) 


F(oj )  f 


4iipc^ 


{' 


RS 


F(a))  .  I  ^  r 
4iip  Cg 


I  aniitj[Hn,(rkr)eJ“<*  +  Hn(rhkr  )e^”^h] 


» 


X  cos  p(z-ZQ)dp  ,  (75) 
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where : 


Rg  -  [(x-xg)^  +  (y-yg)2  +  (z-Zg)^)^^^; 

Rh  ■  [(x-Xg)2  +  (y+yg-2h)2  +  (z-Zg))^^^; 

r  -  (x2+y2)^^^; 

Th  ■  (4h2+r2-Arh  cos 


'J’h 


sln”^  (—  sin  ; 


and 


^  -  tan“^  (x/y). 


The  corresponding  transient  shear-vave  displacement,  u*^(r,t), 
may  be  represented  as  the  sum  of  the  Incident  field  displacement,  u^^nr,t), 
and  the  scattered  field  displacement  as  follows 


u*^(r,t)  *  u(i)(r,t)  +  u(®Hr,t)  . 


The  Incident  field  Is  given  by 


u(i) 


1  f(t-Rs/Cg)  Kt-R^/Cg), 


4jipc^ 


{ 


R®. 


(76) 


i>Si 


The  scattered  field  can  be  expressed  In  the  form  of  a  two- 
dimensional  Fourier  transform  over  frequency  and  axial  wave  number  as 


ut(r,t)  -  r  r[F(a))T(u),p)]e^^^^~^o^e-Jwtdpdw  (77) 

C2n  -CO 


where  F(u))  Is  the  Fourier  transform  of  f(t)  and  T(uj,p)  Is  given  by 


T(a),p)  - - -  I  a„,xj[Hn,(rkr)ej®l'  +  H„,(rhkr)e^'^’']  •  (78) 
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The-  integral  representation  for  the  scattered  field  in  Equation  (77)  may  be 
evaluated  numerically  and  the  result  added  to  the  incident  field  given  by 
Equation  (76)  to  obtain  the  total  field.  A  numerical  method  of  Integration 
used  by  Tsang  and  Rader  (1979)  may  be  implemented  to  evaluate  the  Integral  in 
Equation  (77).  This  integration  method  is  referred  to  as  real-axis  integra¬ 
tion  because  a  Laplace  contour  is  used  in  the  frequency  integration,  while  a 
path  along  the  real  axis  is  used  in  the  wave  number  integration.  Using  this 
method,  singularities  in  T(<»),3)  are  located  away  from  the  integration  path  and 
aliasing  problems  are  reduced. 

In  this  method,  the  integral  over  frequency  is  performed  first 
for  each  complex-valued  frequency.  The  p  integration  is  given  by 

b(a),z)  =  . 


The  result  is  then  Laplace  transformed  over  frequency  by  evaluating  the 
integral , 


0O+j(i)J 

u^®^(r,t)  “  /  F(u))  b(u  ,2  )e~j‘^^dj  , 

^  -oa+jWJ 


where  wj  >  0. 

2,  Scattering  of  Cylindrical  SH  Waves  by  a  Cylindrical  Cavity 

The  steady-state  solution  for  elastic-wave  scattering  from  an 
infinitely  long  cylindrical  cavity  in  a  homogeneous  lossless  half-space  for 
incident  cylindrical  SH  waves  proceeds  in  a  manner  similar  to  that  derived  for 
incident  spherical  SH  waves. 

The  incident  cylindrical  SH-wave  displacement  is  expressed  in 
terms  of  cylindrical  functions  as 


(79) 


and  its  associated  image  incident  displacement  function  is 


u(i) 

h 


^  I  H„(rSk)J„(kr)eJ"^^-^h) 


(80) 
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s  s 

where  the  radial  distances  rg  and  r^  and  the  angles  $  and  '^h  have  been  defined 
in  Figure  11-5.  Alternatively,  the  scattered  cylindrical  SH-wave  displace¬ 
ments  and  are  represented,  respectively,  by 


and 


u<s)  -  I  AnHn(rk)eJ"^ 
n»-«o 


h  n*-® 


(81) 


in  which  r^  (4h^  +  r^  -  Arh  cos 

Using  the  addition  theorem  given  in  Equation  (71),  the  result¬ 
ing  total  SH-wave  displacement  may  be  written  as 


u^ 


7^  I  «j[Hn(rsk)e~^*^s  +  Hn(r^k)e'J"<*hj„(kr  )eJn'>] 
^wpCg  n-^ 


I  AjiHjj(rk)eJ™<‘ 

n=-® 


+  I  I  AnJn,(kr)Hn+^(2kh)eJ”^  . 


(82) 


After  applying  the  traction-free  boundary  condition  at  the  surface  of  the  air- 
filled  cavity,  the  following  Infinite  set  of  linear  equations  is  obtained: 


I  am  [*nm  +  !  Hn+„(2kh)]  »  [  Hn(krh)e'J"^h  +  Hn(krs)e’^"'^s]  ; 

1=-®  H^(ka)  H^(ka) 


for  n  *  0,±1,±2,...  ±  ®  ,  (83) 


where: 


*m 


F(u)) 

r 


for  k  ^  0 


(nj) 


0  ;  for  k  »  0 
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and 


^  nm 


1; 

0; 


for  n  -  m 
for  n  ^  m. 


Finally,  the  steady-state  solution  for  the  total  cylindrical  SH-wave  displace¬ 
ment  at  the  measurement  point,  P,  expressed  in  terms  of  the  scattered-wave 
coefficients,  a^,  is  given  by 


u(t)(a,)  ==  --—-j  (TtjHo(kri)  +itjHQ(kr2)l 
4npc^ 


I  an,xj[Hn,(kr)eJ“<'  +  H^(krh)ej”^'i]  (84) 

4Ttpc|  ni=-« 

where : 

ri  »  [(x-xg)^  +  (yXs)^!^^  * 

r2  “  t(x-X3)2  +  (y+ys“2h)^]  ; 

r  -  U2■^y2)‘/^ 

1/2 

r^  =  (^h^  +  r^  -  4rh  cos  (> )  ; 

>•  sin“^  (—  sin  $)  ;  and 
>•  tan”^  (x/y). 

The  corresponding  transient  shear-wave  displacement  becomes 


u(^^(r,t)  -u^^^(ri,t)  +  u^^)(r2.c)  +u(®nr,t)  . 


(85) 


The  scattered  field  displacement,  u^®^(r,t),  can  be  expressed  in  the  form  of 
one~dimenslonal  Fourier  transform  over  frequency  as 


a 


u(®^(r,t)  ■  “  f(t)  *  Tj^(t) 


(86) 


whe  re : 


Tjl((i))  -  -  -/  —  V  I  anjiij  [Hnj(rk)eJ"^  +  H^(krh)e"“*’^] 
m=— 


(87) 


and 


Ti(t)  =  ^  r  Ti(u>)e-j“tda,  . 


(88) 


The  Incident  field  displacements,  u^^nri,t)  and  u(^^(r2>t),  can  be  deduced  by 
applying  the  Cagnlard-deHoop  method  (deHoop,  1960)  for  solving  seismic  pulse 
problems.  Each  of  the  incident  SH-vave  displacements  given  in  Equation  (85) 
can  be  obtained  from  the  function  represented  by 


c 

u^^Hr.t)  *  /  f  (t-t  )g(r ,T  )dT  ,  (89) 

0 

where: 


(90) 


is  regarded  as  the  wave  function  corresponding  to  the  delta  function  time 
dependence  of  the  SH-line  source. 
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F.  SH-Wave  Scattering  from  a  Cylindrical  Cavity  in  a  Two-Layer 

Half-Space 

The  steady-state  and  transient  solutions  for  elastic-wave  scattering 
from  an  infinitely  long  cylindrical  cavity  in  a  two-layer  homogeneous  lossless 
solid  medium  is  developed  for  incident  cylindrical  SH  waves.  The  cylindrical 
SH  waves  are  generated  by  a  line  source  within  the  two-layer  half-space  medium. 
Formal  analytic  expressions  are  developed  for  the  particle  displacement  of 
scattered  waves  at  any  detector  position  in  the  host  medium  for  an  SH-wave  line 
source  oriented  parallel  to  the  axis  of  the  cavity.  Analytical  expressions  are 
derived  for  SH-wave  displacements  to  investigate  the  effect  of  the  surface 
layer  on  the  detection  of  SH  waves  scattered  by  the  cylindrical  cavity.  In 
addition,  to  compute  particle  velocity  seismograms,  an  expression  for  the  time 
derivative  of  the  total  SH-wave  displacement  is  presented  and  the  corresponding 
transfer  function  is  included  as  part  of  the  computed  sythetic  seismograms. 

1.  Formulation 


The  derived  theoretical  solution  is  in  the  form  of  an  integral 
equation,  the  unknown  being  the  total  SH-wave  displacement  on  the  surface  of 
the  cylindrical  cavity.  The  integral  equation  requires  the  two-layer  half¬ 
space  Green's  function  and  the  incident  SH-wave  displacement  for  a  line  source 
to  be  evaluated  at  the  surface  of  the  scattering  target.  Then,  the  Integral 
equation  is  reduced  to  a  matrix  equation,  which  can  be  solved  numerically  by  a 
direct  method.  Thus,  the  SH-wave  displacements  outside  of  the  cylindrical 
cavity  are  found  by  integrating  the  two-layer  half-space  Green’s  function  on 
the  surface  of  the  cavity.  Summed  with  the  Incident  SH-wave  displacement,  the 
combined  results  represent  the  total  SH-wave  displacement  at  any  detector 
position  in  the  two-layer  medium. 

2.  The  Two-Layer  Half-Space  Green's  Function 

A  stratified  earth  representative  of  a  weathered  surface  layer 
and  a  bedrock  semi-infinite  half-space  may  be  simulated  by  a  two-layer  model. 
This  model  consists  of  a  bedrock  cavity  host  medium  of  density,  P2»  shear 
modulus,  p2»  shear-wave  velocity,  C2»  overlain  by  a  surface  layer  of 

thickness,  h,  density,  P^,  shear  modulus,  m,  and  shear-wave  velocity,  cj,  in 
which  Cl  is  less  than  C2« 

An  SH-wave  line  source  is  located  at  the  position  P'(x',y'),  and 
the  radiated  SH-wave  particle  displacement  is  detected  at  the  position  P(x,y). 
Both  the  line  source  and  the  measurement  point  are  located  in  the  bedrock  layer 
as  shown  in  Figure  II-6. 

The  steady-state  Green's  function  for  SH-wave  particle  motions 
satisfies  the  wave  equations 


(V^+ki)Gi(x,y;x' ,y ' )  “  0  ;  for  0  <  x  <  h  , 


and 

(7  2+k2)G2(x,y;x' ,y ' )  "  (x-x '  )6  (y-y ' )  ;  for  x  >  h  ,  (91) 
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FIGURE  I 1-6,  TWO-LAYER  HALF -SPACE  HAVING  AN  SH-WAVE  LINE  SOURCE 
AND  A  DETECTOR  LOCATED  IN  THE  BEDROCK 


where : 


2  2 

V2  =  ;  and  ki  =  ^  ;  k2  =  ^  in  which  ci  <  c2  . 

ox'‘  dy^  C2 

The  Green's  functions  for  cylindrical  SH-wave  motions  may  be 
developed  using  the  spectral  representation  of  the  free-space  Green's  function 
developed  in  Appendix  A.  Thus,  the  Green's  functions  for  the  regions  of  dif¬ 
ferent  elastic  properties  are: 


For  the  region  0  <  x  <  h: 


For  the  region  x  >  h: 


(92) 


r.  =  ^  -■^2|x-x'l  ^  1  -^2X1  e  JX(y-y’) 

G2  =  Le  '  +  f2(X)e  ^  j  fT  dX  , 


(93) 


where: 


and 


n  -  a2-k,2)*«  ; 
'>2  -  . 


40 


The  determination  of  the  square  roots  is  such  that  Real  {Tj[,Y2}  >  0  over  the 
entire  path  of  integration,  <  X  <  The  factors  fi(X),  gi(X),  and  f2(X) 
are  arbitrary  functions  of  X  which  render  the  integrals  corveigent  and  satisfy 
the  following  boundary  conditions: 

At  the  surface  of  the  earth: 


dx  -'x=0 


=  0  . 


(94) 


At  the  layer  interface  (continuity  of  stress  and  displacement 

are  required): 


and 


x=h 


^  ^  x=h 


(G^]jj=h  “  ^^2Jx=h  • 


(95) 


These  relationships  lead  to  three  linear  algebraic  equations  to  be  solved  for 
the  unknown  coefficients  fi(X),  gi(X),  and  f2(X).  Solutions  of  the  Green's 
functions  for  the  surface  layer  and  the  bottom  semi-infinite  region  are: 


For  the  region  0  <  x  <  h: 


,  ,  M  1  ^  e  ^  -Yi(h-x), 

l(x,y:x’,y')  = J  ei(X)  —  {e  ^  +e  '  ) 


X  eJ^<y-y')dx  . 


(96) 


For  the  region  x  >  h; 


n  /  1  ^  f*f  “Y2|x-x'|  (.  -Y  2(x+x'-2h)i 

G2(x,y;x' ,y ' )  - J  |e  '  +e2(X)e  ^  | 


X 


eJX(y-y’) 

tT  " 


(97) 
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where: 


and 


in  which 


and 


ei(X) 


T(\) 


1  -  R(\)e 


e2(X) 


-R(X)  +  e 


+  I*' 


1  -  R(X)e"^l^ 


T(X) 


2^1 2y  2 


+  ^^2Y2 


R(X) 


lilYi  -  ^i2Y2 

^lYi  +  ^2^2  * 


(98) 


(99) 


(100) 


(101) 


To  represent  cylindrical  SH-wave  particle  motion  in  terms  of 
outgoing  waves  at  infinity,  the  Green's  functions  given  in  Equations  (96)  and 
(97)  are  expressed  in  polar  coordinates  (r,^).  For  convenience,  Y and  Y2  are 
replaced  in  Equation  (97)  by  -jai  and  -ja2,  respectively.  Thus,  the  bedrock 
Green’s  function,  G2,  may  be  rewritten  as 


r  f  f”/  ja2(x~x')  j.  'i  ja2(x+x’-2h), 

G2(x,y;x',y')  =  j  {e-^  ^  +  e2(X)e-^  ^  } 


ejx(y-y*) 

“2 


dX 


(102) 


where : 


«,  .  (k,2  -x2)>/2  i 

®2  “  (k2^  ~  X^)^/^  . 


Figure  II-7  shows  the  appropriate  geometry  for  developing  a 
Solution  for  cylindrical  SH-wave  motions.  This  solution  may  be  expressed  in 
terms  of  a  cylindrical  wave  function  having  its  origin  at  the  position  (H,d). 
Making  the  substitution  X  =  k2  sin  6  in  the  second  term  of  Fquation  (102) 
leads  to  the  Green's  function  in  polar  coordinates. 
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i.} 


(105) 


into  Equation  (103)  yields  the  cylindrical  mode  expansion  for  the  Green's 
function  in  the  cavity  host  medium.  That  is, 


G2(x,y,x',y')  - I  Jn(k2r  ’  )e"J"'J* '{ jTiHjj(rk2)eJn't' 

n=-a> 


+  Jit  I  ejn»>Jo(rk2)PnJ  (106) 


where: 


1  ®  .  2ja2(H-h)  +  j0(m+n) 

Pnm  *  J  £2^)  f-  dX  (107a) 

-«o  u2 


in  which 


6  *  tan“^(\/a2) 


or  where: 


f„«  -  ,  (I07b) 


where  the  Integration  path,  F,  is  described  in  Appendix  D. 


Following  the  procedure  used  in  deriving  the  bedrock  Green’s 
function,  the  surface-layer  Green's  function,  Gj,  in  terms  of  aj  and  a2  may  be 
written  as 


Gi(x,y;  x’,y') 


4n 


/  ei(X) 


gja2(x'-h) 

ai 


jgjai(x+h)  ^  gjai(h-x)j 


X  eJ^^y-y’)dX  .  (108) 


Making  the  appropriate  substitutions  to  convert  Equation  (108)  to  polar 
coordinates,  the  surface  layer  Green's  function,  G^,  becomes 


Gl(x,y;x',y')  -  I  jxVr ’k2 )e-J"4>  '  ^  eJ^^* Jn(rk2)Qnm  (109) 
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whe  re : 


^  -«o 


in  which 

9  “  tan~^  (X/02)  (110a) 


or  where: 


71  p 


X  ^gjk2«(9)(x+h)  ^  ^jk2a  (9  )(h-x)  j 


(110b) 


in  which: 

a (9 )  =  (k|/k|  -  sin^0)^^^. 

3,  Incident  SH-Wave  Displacement  for  a  Line  Source  in  a 
Two~Layer  Half-Space 

The  incident  SH-wave  displacements  for  a  line  source  located  at 
(Xs»ys^  seml“inf inite  bottom  layer  satisfies  the  wave  equations 


(7^  +  =  0;  for  0  <  x  <  h 


(7^  +  k2^u2^^^)  =  ---— j  6(x-xs)6(y-ys^»  x  >  h  , 

P2C2^ 


(111) 


where  k^  and  k2  are  defined  in  Equation  (91).  The  function  F(uj)  is  a  body 
force  acting  parallel  to  the  z-axis  at  the  position  (xg.ys)*  Since  the  wave 
equations  given  in  Equations  (111)  and  (91),  respectively,  have  the  same  forms 
and  they  are  subject  to  the  same  boundary  conditions,  the  incident  SH-wave 
displacements  may  be  deduced  directly  from  Equations  (96)  and  (97). 

For  the  region  0  <  x  <  h: 


^  f”ei(\)  ^\gjai(x+h)  ^  ^jai(h-x)i 

^  2^2^  ^  ®  1 


X  e*’' >dX  . 


(112) 


4.  Integral  Equation  Solution  for  SH-Wave  Scattering  from  a 

Cylindrical  Cavity  In  a  Two-Layer  Half-Space  for  a  Detector 
In  the  Bedrock 

The  total  SH-wave  displacement,  In  a  medium  Is  composed  of 
two  parts;  the  Incident  wave,  and  the  scattered  waves,  That  Is, 


ut  m  y(l )  +  ) 


(118) 


Each  of  these  wave  functions  satisfies  the  first  (Interior)  and  second 
(exterior)  Helmholtz  equations  as  given  In  Appendix  B.  For  a  cylindrical 
alr~fllled  cavity,  the  normal  derivative  of  the  total  SH-wave  displacement 
vanishes  at  the  cavity  surface,  S  (i.e.,  the  stress-free  surface  boundary  con¬ 
dition  Is  satisfied).  Under  this  condition,  the  total  SH-wave  displacement, 
u^,  at  the  surface,  S,  may  be  obtained  from  Equation  (B-7)  and  it  satisfies 
the  integral  equation 


du(l^(r)  a  f  (-/  ,x  dG(r,r')  , 

- -  =  — -J  u'-(r’)  - rt -  dS  ;  for  r  on  S 

Sr  dr  S  5r 


(119) 


Next,  the  Green's  function  given  by  Equation  (106)  and  the 
Incident-wave  displacement  as  given  by  Equation  (115)  together  with  the  total 
SH-wave  displacement  on  the  boundary  S(r'  ■  a),  assumed  to  be  of  the  form 


u^  »  ^  BjieJ-^*>  '  , 
i=-« 


(120) 


are  substituted  into  the  integral  equation  given  by  Equation  (119).  With 
dS  =  ad(>',  where  a  is  the  radius  of  the  cavity,  the  Integral  equation  takes 
the  form 


bu2^^^ 


dr 


1  d  - 


^  {jit  I  Jn'(k2a)e-Jn^'[Hn(rk2)ejn^  +  I  Jni(rk2)eJ"»PnJ}  .  (121) 


m= 


Utilizing  the  properties 


I  2it  (  for  ^  *  n 


1;  for  I  *  n  , 
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Eqiration  (121)  (at  r  “  a)  is  reduced  to 


5u2^^  ^ 
5  r 


“  I  '(k2a){% '(ak2)ei'^<*  +l  dm' •  (122) 

r®a  s-a» 


(122)  yields 


Substituting  the  incident  wave  of  Equation  (115)  into  Equation 


-  I  Jn’(k2a)eJn4»[H„(rsk2)e"j"‘>s  +  sj} 

2^2  n*“® 


fl  jTt(BjiJi’(k2a))[Hi‘(ak2)eji4>  +  ^  Jm’(ak2)ej'n>PiJ  . 
A  ■-«o 


(123) 


reduced  to 


After  some  algebraical  manipulations,  Equation  (123)  may  be 


*  J  ’(akol 

I  [f  B„J„'(ak2)J{6n„  +  Pnml 

•na-«  n  fc' 


_F((jO_  k,)e”^"'^s  +  S  1 

4np2C22  ®  ^  ’ 


(124) 


where : 


Sn  -  I 

A*-® 


Then,  making  the  substitution  of 


f  Jm'(a>^2) 


In  Equation  (124),  the  linear  system  of  equations  for  the  unknowns  may  be 
expressed  as 


1  ^  nm 


+  ^n'<ak2) 
Hn’(ak2) 


^nml 


F(u))  Jn*tk2a) 
4itp2C2^ 


(125) 


for  n  ■  0,±1,±2,...  ±  ®  . 


After  the  value  for  the  total  SH-wave  displacement,  u^,  at  the 
surface,  S(r*a),  Is  known,  the  scattered  waves  In  the  field,  r  >  a,  may  be 
derived  from  the  second  term  of  Equation  (II-9)  as 


.  dG2(r,r’) 

/  u^(r')  - X -  dS;  for  r  in  V 

S  ar' 


(126) 


The  integration  of  Equation  (126)  can  be  carried  out  in  the  same  manner  as 
before  with  the  result 


I  I  jitJn'(k2a)e-j™>’)d(>’ 

0  j{,=-<«>  n=-® 


X  [Hn(rk2)eJn^  +  I  Jm’ (r'^2)eJ'^Pnm] 

|QS««0 


or 


u(8)  -  jx  I  [f  BnJn’(k2a)][H„(rk2)ejn^  +  I  VrW2)eJ“>Pnm]  •  (^27) 


In  terms  of  the  wave  coefficients,  bn»  Equation  (33)  Is  expressed  as 


(jitbn)[H„(rk2)ejn^  +  I  Vr»^2)eJ'^Pnm]  •  <128) 
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a.  Steady-State  and  Transient  Solutions  for  the  SH-Wave 
Particle  Motion 


Superposition  of  the  incident  and  scattered  SH-wave  dis¬ 
placements,  and  given  by  Equations  (113)  and  (128),  respectively, 

leads  to  the  steady-state  solution  for  the  total  SH-wave  displacement  for  the 
case  where  the  line  source  and  the  detection  point  are  both  located  in  the 
bedrock.  That  is. 


U2^(oj)  a  {TijHo(k2ri)  +  j/ 

4iip2C2^  “2 


where : 


E  (jiibn)[Hn(rk2)eJ"^  +  I  •Jm('-k2)eJ‘^PnJ 
n=-«“  m=-® 


(129) 


1  /  ^ 

ri  =“  (r^+rg^  -  2rrs  cos  (ij>-1)s)]  '  ; 


r  -  ((x-H)2  +  (y-d)2)^/2  . 

*■3  “  [(xs-H)2  +  (yc“<i)^l^^^  ; 


*  tan 


-1  (y-d)  . 

x-H  ’ 


=  tan 


-1  -^ys-d)  . 
(x,-H)  ' 


®1  *  (ki2-\2)^/2  . 


ai  -  (k22-x2)l/2 


To  include  attenuation  and  dispersion  in  Equation  (129),  the  complex  propaga¬ 
tion  constants  for  the  surface  layer  and  bedrock  are 


kl 


U) 

Cl((i)) 


+ 


ja  I  (u) ) 
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and 


^"2  -  ^ 


in  which  the  absorption  coefficients  and  the  phase  velocities  are 


a 


2QjCj 


and 


-1 


for  j  *  1,2. 

4 

To  evaluate  the  Integral  Pnm  In  Equations  (124)  and  (129),  a  numerical  tech¬ 
nique  is  developed  in  Appendix  C.  In  addition,  for  the  case  of  scattering  from 
deep  cavities,  an  approximate  solution  of  the  integral  using  the  saddle- 
point  method  is  derived  in  Appendix  D. 

The  transient  SH-wave  displacement,  u^(t),  corresponding 
to  the  steady-state  solution  given  in  Equation  (35).  may  be  represented  as  the 
sum  of  the  incident  time-dependent  displacement,  u2^^Ht),  and  the  scattered 
displacement,  U2^®Ht).  That  is. 


U2^(t)  ••  U2^^^(t)  +  U2^®nt) 


(130) 


The  solution  of  the  incident  transient  SH-wave  displacement,  U2^^^(t),  for  a 
line  source  located  in  the  bedrock  at  (xs,yg)  is  derived  in  Appendix  E.  Thus, 
for  a  transient  SH-wave  pulse  acting  along  the  z-axis  parallel  to  the  cylin¬ 
drical  cavity  and  consisting  of  a  body  force,  f(t),  the  Incident  transient  SH- 
wave  displacement  is  expressed  as 


.-xax,.  T(i)(t)  . 

2'tip2C2^ 


(131) 


and  its  time  derivative  is 


du2^^ )(t ) 
dt 


d  ^  f (t) 
2i;p2C22 


*  T(i)(t))  , 


(132) 
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where.: 


T(i)(t) 


H(t-ts) 


+  ^  Im 
1*0 


(p)  dp 
CLi  dt  j 


H(t-tji+i ) 

P=‘Pl+l(t) 


“  I 
1=0 


R^'*~kp)  dp 
02  dt 


H(t-tii) 

p=pjl(t) 


(133) 


The  first  term  of  the  transfer  function,  T^^-^Ct),  represents  an  incident  SH 
wave  traveling  directly  from  the  source  to  detector  with  a  time  of  arrival 
given  by 


tg  =■  [(x-Xg)2  +  (y-yg)2]2/C2 


(134) 


The  last  two  terms  of  Equation  (133)  represent  reflections  in  the  bedrock  and 
trapped  reflections  within  the  surface  layer.  Both  terms  are  characterized 
by  the  reflection  coefficient  R(p)  and  the  arrival  times  tj^  and  tj^+j  (for 
1  =  0,1,2,...  N).  These  arrival  times  are  given  by 


(x+X3-2h)  ^  2hl 

C2  cos  9  2(^  )  Cl  cos  0lU) 


and 


(x-hcs-2h)  ^  2h(l-H) 

C2  cos  0  2^'*^^^^  C2  cos  9i(^'*’^^ 


(135) 


(136) 


where  92^'^^*  02^"^^^  are  the  angles  of  the  ray  trajectories  in 

the  surface  layer  and  in  the  bedrock. 

The  scattered  field  displacement,  u2^®Ht),  can  be 
expressed  in  the  form  of  a  one-dimensional  Fourier  transform  over  frequency  as 


U2(3)(t)  -  /"e">tF(u))T(s)((j)d(j  -f(t)  *T(s)(t) 


(137) 
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whe  re : 


t(s)(u).-- - (jiibn)[Hn(rk2)eJnt  +  Jn,(rk2)ej“«l>Pn  J  (138) 

^P2^2  n**^ 


and 


T(s)(t)  - j“T(s)(a,)e-J‘*>tdu, 


(139) 


The  corresponding  time  derivative  for  the  scattered  field  displacement, 
Is  given  by 


du2^®^(t) 

;  dt 


(f(t)  *  T(s)(t))  . 


(140) 


b.  Check  of  the  Steady-State  Solution 

In  order  to  check  the  general  SH-wave  displacement  solu¬ 
tion  given  by  Equation  (129),  the  half-space  solution  for  scattering  SH  waves 
from  a  cylindrical  cavity  may  be  derived  directly  from  Equation  (129)  by  let¬ 
ting  h  0.  In  fact,  in  that  case  the  function  e2(?^)  •*'  1  and  the  second  term 
of  the  Incident  SH-wave  displacement  becomes: 


iijHo(k2r2) 


(141a) 


where: 


r2  “  [(yys)^  (x+Xs)^]^^^  . 


In  addition,  the  Integral,  Pnm>  scattered  SH-wave  displacement  term  of 

Equation  (129)  and  In  the  linear  system  of  equations  given  by  Equation  (124) 
as  h  ->•  0  are  reduced  to  the  Hankel  function.  That  Is,  for  h  +  0  and  e(X)  +  1, 


nm 


(-J  i  ;  S  *  -  W2k2»)  • 


(141b) 


Using  the  addition  theorem  for  the  Hankel  function. 


Hn(Rk2)e±Jn'I^  -  f  Jm(rk2)Hn-hi>(2k2H)e± . 
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the  function,  Sjj,  In  Equation  (124)  becomes 

Sn  =  I  =  Hn(rs«k2)e‘J"‘^s“ 


where : 


cos  5 


’=‘s  “  [(ys'd)^  +  (xs-H)2]^'^^  ; 


=  sin“^(sin  (tig  1  and 

r  e 


■  (^7h>  • 


In  the  same  manner,  the  second  term  of  Equation  (129)  is  reduced  to 


Dn  “  I  eJ“«>J„(rk2)H„+^(2k2H)  = 

m— 

where : 

-  4rH  cos  ; 

r  >  [(y-d)2  +  (x-H)2i^/^  ; 

<>  =  tan“^  ;  and 

x-H 


(|>H  -  sin“l[-^  sin  ♦] 


Thus,  the  final  solution  for  scattering  of  cylindrical  SH  waves  by  a 
cal  cavity  in  a  half-space  is 


I  ^m[  ^  nm 


Jn'(ka) 
■"  Hjj*(ka) 


Hn+in(2kH)] 


F(a)) 


4npc' 


{Hjj(krg^)e  +  Hn(krg)e 


for  n  ■  0,± 1,2,...  t  ® 

S5 


(142) 


(143) 


cylindrl- 


(144) 


and 


u(w)  -  •^•^'2  +  TijHo(kr2)] 


+  I  (iijb^)[H^(kr)ej">^  +  H„(krH)eJ'^»]  .  (145) 

m*"® 

This  result  Is  consistent  with  that  derived  using  the  image  theory  developed 
earlier  for  a  cylindrical  cavity  in  a  homogeneous  half-space. 

c.  The  Steady-State  Approximate  Solution 

For  small  air-filled  cavities  located  at  large  distances 
from  the  layer-bedrock  interface,  the  steady-state  SH-wave  solution  given  by 
Equation  (129)  can  be  simplified.  In  fact,  when  the  total  SH-wave  displace¬ 
ment  satisfies  the  stress-free  boundary  condition  at  the  surface  of  the 
cavity,  the  effect  of  the  surface  layer  may  be  neglected  for  a  detector  at  the 
position  r  *  a.  In  this  case,  the  effect  of  the  second  term  in  the  scattered 
SH-wave  displacement  is  much  smaller  than  that  of  the  first  term.  Thus,  under 
this  condition  the  scattered  SH-wave  coefficients  at  the  surface  of  the  cavity 
may  be  determined  using  the  following  expression 


bn 


F(a))  I 

4np2C2^  ^'^^2®^ 


[Hn(rsk2)e 


-jn(l>i 


+  ^  e  ^'^^®Jji(rsk2)Pin]} 
Jl  =->» 


(146) 


To  obtain  the  scattered  SH-wave  displacement  at  the  detec¬ 
tor  location,  the  coefficients,  bji.  given  by  Equation  (146)  are  substituted  in 
Equation  (128)  to  obtain  the  steady-state  scattered  SH-wave  displacement 


"  ~  '%p2c'2^  ^  'H^(k2r)~  ^”n(rsk2)e  +  Sn][  Hn(rk2)ejn^  +  Dj,]  (147) 


where: 


Sn  =  I  e  5 

i=-® 


“  I  ®  )Pnm  * 

m=-® 
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r  =  [(x-H)2  +  (y-d)2]l/2  . 


(j)  =  tan"l(^)  ;  and 


>s  =  tan'H^^)  • 


d.  The  High-Frequency  Limit  Approximation 

To  reduce  the  computations  necessary  to  produce  synthetic 
seismograms  and  as  a  check  on  the  theoretical  analysis,  a  high-frequency  limit 
approximation  for  the  scattered  SH-wave  particle  displacement  is  derived. 

This  solution  may  be  obtained  by  substituting  the  function  given  by  Equa¬ 
tion  (D-17)  into  Equation  (128)  to  yield 


=  I  (jubn)  (Hti(rk2)ej'i'>  -  Ro^n(*^2ra)e^^‘^® 


+  I  a,Ro^(Ro‘^  -  Ro)HniR2ra^^^]e^"^3^ 


(148a) 


where : 


s;2v  -  1 
^2v  +  1  ’ 


rg  =  (r^  +  4(H-h)^  -  4r(H-h)  cos  <>  ] 


,-i  y-~  d  . 

^  X  -  h  ’ 


r  1  -t-  hJlv/(H-h).  1/2 
1  +  hJl/v(H-h)J 


•((g  -  sln"^  sin  (P)  ; 

t  a 


57 


^(A)  m  [^2  +  _  2sjjr  cos 


and 


In  which 


♦  -  sln"^  ("  sin  (J>)  ; 


8j^  *  2(H  -  h  +  htv )  ; 


^2  -  m/ti2  ; 


and 


V  =  ki/k2  . 


The  numerical  solution  of  Equation  (148a)  for  the  scattered  SH-wave  particle 
displacement  is  produced  by  the  computer  program  scat_pulse  in  Appendix  H. 

The  scattered-wave  coefficients,  b^,  may  be  obtained  from 
Equation  (146).  Here,  the  source  function,  S^,  in  the  high-frequency  limit 
approximation  is  given  by 


Sn  -  -Roe 


HnU2rs^°^l  +  I  Hnlk2rs(^>) 

A-1 


(148b) 


where: 


fJ  »  (Rq-I  -  Ro)Ri 
rg^°^  “  f^s^  4(H-h)2  -  4(H-h)rg  cos 
rg(A)  B  (rg2  +  g^2  _  4(H-h)sji 
(l;g(o)  “  sln'l  ^77^ 

=  sin“^  ( ~(X)  • 
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In  a  similar  manner,  the  high-frequency  limit  approximation  for  the  incident 
SH-wave  displacement  may  be  derived  by  substituting  the  function  Pnm  in  Equa¬ 
tion  (VI-17)  into  Equation  (115)  to  obtain 


U2^^  ^ 


F(u)) 

%p  2^2^ 


(itj)  (  Ho(k2|  r-rgj  )  +  I  e  ^“*’sj^(l5_2rg) 


X  I  ejn<'Jn(rk2){-RoHn-hDl2k2(H-h)l 

n=-<*‘ 


^  ^ ”Ro  ^Ro^n+Tii^R2SJl 

A=1 


(149) 


or 


U2(i)  =  (^j)  /  Ho(k2|r-rs|  )  -  Rq  I  e^"^'*'’^'^®^Jni(k2rs)Hn(k2rb) 

4iip2C2^  ^  m=-« 


+  I  aji(Ro-l-Ro)RoA  J  e^"^‘‘'^^^^"^s\(k2rs)Hn(k2rb^^^)}  •  (150) 

i.=l  n=-®  ^ 


Finally,  Equation  (150)  is  further  reduced  to  yield  the  SH-wave  particle 
displacement 


U2(l)  -  (uj)  i  Ho(k2|r-rsl  )  -  Rollo(R2rt) 

4np2C2^  V 


+  I  aji(Ro''^-Ro)RiHo[k2rt(i>)l  .  (151) 

A-l  1 


where : 


fb 


-  (r2  +  A(H-h)2  -  Ar(H-h) 


cos  0  1 


1/2 
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4»l,  -  sin'l  sin  41) 

^  “  Ir^  +  sjj^^  -  2sj^r  cos 

4,jj(A)  .  8ln“^  ^"TTIT 

b 

1  /o 

Tf.  -  Irg2  +  rb^  -  2rsr5  cos  (4»b-^s)]  ^ 


and 


r^U)  -  {tg^  +  +  2rsrb^-^^  cos  [(i^b^-^  ^-(J)sl}  • 


The  numerical  solution  of  Equation  (151)  Is  produced  by  the  computer  program 
lnc__pulse  In  Appendix  H.  To  simplify  the  calculations,  the  geometric  rela¬ 
tionships  given  In  Equations  (148)  and  (151)  may  be  deduced  directly  from 
Figure  11-9.  In  terms  of  the  detector  and  Its  Image  coordinates,  the  dis¬ 
tances  and  angles  included  in  the  SH-wave  scattering  equation  are  given  by 


r<o)  i  rg  =•  [(y-d)2  +  (x+H-2h)2] 


4,(0) 


*>a 


tan”^( 


■  y.~d  ..-j 
x+H-2h^ 


rCf)  =  *  { (y-d)^  +  (x+H+2(Jlv-l  )hl2} 


and 


,a) 


M) 


tan 


-1 


..y-d 


x+H+2(Jlv-l)h 


■) 


For  the  Incident  SH-wave  displacement,  the  new  distances  are 


rgCo)  =  rt  =  ((y-yg)2  +  (x+Xs-2h)2]^/2 


and 


rs^^  -  ”  t(y-ys)^  +  [x+Xs+2(Jlv-l)h]2}  . 
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(a)  Scattered  Waves  (b)  Direct  Waves 


5 


I 


Integral  Equation  Solution  for  SH-Wave  Scattering  from  a 
Cylindrical  Cavity  In  a  Two-Layer  Half-Space  for  the  Detector 
In  the  Surface  Layer 

After  the  value  for  the  total  SH-wave  displacement,  u^,  at  the 
surface  of  the  cavity  (r  a)  Is  known,  the  scattered  wave  In  the  field 
0  <  X  <  h  may  be  derived  from  the  second  term  of  Equation  (B-9)  for  the  sur¬ 
face  layer  Green’s  functions,  Gj,  as 


ui^s)  «  j  ut(r') 

S 


5Gi(r.r*) 

dr’ 


dS; 


for  r  In  V. 


where  u^(r')  Is  given  by  Equation  (120).  The  integration  of  Equation  (152)  is 
carried  out  in  the  same  manner  as  was  done  for  the  case  of  a  detector  in  the 
bedrock.  That  is, 


I  I  jiiJni'(ak2)e~j“l* '  d(|> ' 

o  m="® 


X  I  eJ"^Jn(rk2)Qnm  (153) 

n=-® 

or 

OD  00 

ui(8)  -nj  I  f  B„J„’(ak2)  I  ejnf Jn(rk2)Qnm  •  (15^) 


In  terms  of  the  wave  coefficients,  bn,»  Equation  (154)  is  expressed  as 


ujCs)  .  I  (jxb„,)  I  ej'»^Jn(rk2)Qnm  ,  (155) 


where  Qnm  1®  given  by  Equation  (110). 

a.  Steady-State  and  Transient  Solutions 

Superposition  of  the  incident  and  scattered  SH-wave  dis¬ 
placements,  ujd^  and  ui(®^,  given  by  Equations  (117)  and  (155),  respectively, 
leads  to  the  steady-state  solution  for  the  total  SH-wave  particle  displacement 
for  the  line  source  in  the  bedrock  and  the  detection  point  in  the  surface 
layer  as 
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uit(u))  =  — I  jiiJn,(rsk.2)e  I  jj^(rlc.2)Qr 

2^2  m="®>  n*“®* 


+  I  (jitbm)  I  eJn^J„(rk2)Qnn,  .  (156) 


where : 


r  =  t(x-H)2  +  (y_d)2]l/2  . 

re  =  [(xs-H)2  +  (y  -d)2]l/2  . 


x-H  ’ 


-1  ^ys~d> 

'  (x,-H)  * 


To  include  attenuation  and  dispersion  in  Equation  (156),  the  complex  propaga¬ 
tion  constants  for  the  surface  layer  and  bedrock  are  given  by 


C2((»)) 


ja  2  («  ) 


in  which  the  absorption  coefficients  and  the  phase  velocities  are 


Cj(u)  -  cj[l  -  l»g  (^)1 


for  j  *  1,2 
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For  the  use  of  scattering  from  deep  cavities,  an  approximate  solution  of  the 
integral,  Qnm»  using  the  saddle-point  method  is  derived  in  Appendix  F. 

The  transient  SH-wave  displacement,  ui*^(t),  corresponding 
to  the  steady-state  solution  given  in  Equation  (156),  may  be  represented  as 
the  sum  of  the  incident  time-dependent  displacement,  ui^^^(t),  and  the  scat¬ 
tered  displacement,  ui^s)(t).  That  is. 


ui^(t)  *  ui(^)(t)  +ui(®)(t)  . 


(157) 


The  solution  of  the  incident  transient  SH-wave  displacement,  ui^^^(t),  for  the 
line  source  in  the  bedrock  at  (xg,yg)  and  the  detector  in  the  surface  layer  is 
derived  in  Appendix  F.  Thus,  for  a  transient  SH-wave  pulse  acting  along  the 
2-axls  parallel  to  the  cylindrical  cavity  and  consisting  of  a  body  force, 
f(t),  the  incident  transient  SH-wave  displacement  is  expressed  as 


liCi)  .  -tlil- .  i(i)(t)  . 

2iip  2C2^ 


(158) 


and  its  time  derivative  is 


dui ^ 

It 


(t) 


1  d 

2iip2C2^ 


(f(t)  *  Ti(i)(t)]  , 


whf  re : 


(159) 


Tjd) 


N 

2  I  14 


1-0 


T(p)R^(p)  dp-j 
a2(p)  dt-l 


p=Pjt(t) 


H(t-tji) 


in  which  pji,(t)  is  a  solution  of  the  nonlinear  equation  given  by 


t  «  p(y-ys)  (xg-h)  +  (— -  P2)^^^  h(2i+l)  . 

Cz  Cj 


Each  term  of  Equation  (159)  represents  reflections  within  the  surface  layer. 
These  terms  are  characterized  by  the  reflection  coefficient  R(p)  and  the 
arrival  times  tjj^  (for  A  -  0,1,2,...  N).  The  equation  to  determine  the  arrival 
times  for  the  detector  at  the  surface  is 


Xs~h  ^  h(21-H) 

C2  cos  0  2^^  Cj  cos  9 j ^ 


(160) 
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tiHfcl 


where  the  angles  Qj  and  02  are  the  angles  of  the  ray  trajectories  in  the  sur¬ 
face  layer  and  in  the  semi -inf inite  half-space,  respectively. 

The  scattered  field  displacement,  Ui^®^(t),  can  be 
expressed  in  the  form  of  one -dimensional  Fourier  transform  over  frequency  as 


ui(s)(t)  =  2^“^  e"JwtF(u,)Ti(s)(t^)da,  -  f(t)  *  Ti<s)(t) 


(161) 


where: 


- - 2  ^  I  e3®>Jn(rk2)Qnm 

^xp2C2  m*"^ 


and 


Ti<s)(t)  “ . 


The  corresponding  time  derivative  for  the  scattered  field 
displacement,  ui^®^(t),  is  given  by 


dui(t) 

dt 


^  tf(t)  *  Ti(s)(t)I  . 


(162) 


b.  The  High-Frequency  Limit  Approximation 

The  high-frequency  lijit  approximation  of  the  scattered 
SH-wave  displacement  for  a  detector  at  the  surface  of  the  earth  is  obtained  by 
replacing  the  approximate  solution  of  as  given  by  Equation  (F-15)  into 
Equation  (155).  The  resulting  particle  displacement  is 


ui(s)  =  2  ^  bjlToRo  I  I  Jn(rk2)Hn+n(k2rji )  (163) 

1  *0  m*-®  n®-® 


where: 


ri  »  (2H-h)  +  (2i+l)hv 

X  ■  _ - _ 

°  +  1 
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R 


o 


-  1 
+  1 


fi  ■  (2i  -*•  l)vhT^/^ 
I-  ^  2H  -  h  -I 
1  ,  (21  +  l)h/v 
2H  -  h 


and 


V  *  • 

The  solution  given  by  Equation  (163)  can  be  further 
reduced  by  Introducing  the  Identity 


Hn(k2Ri)e^“^^  *  I  Jn(rk2)Hn+m('^2rjl )  •  (164) 

n=-*“ 


Therefore,  by  substituting  Equation  (164)  into  Equation  (163),  the  SH-vave 
displacement  becomes 


1 


2 


I  a^TgRo-^ 

A-O 


I  (jTibni)Hnj(k2rjt)e'^“*'^ 


(165) 


where: 


fi  «  [r^ 


2  _ 


2rrji 


cos 


4>] 


1/2 


r  -  [h2  +  (y-d)2]^/^  ; 


4i«  *  sin“l  f-ir-  sin  ; 
^A 

=•  V^lhl 


and 


V  “  ki/k2  . 


The  numerical  solution  of  Equation  (165)  for  the  scattered  SH-wave  displace¬ 
ment  is  calculated  by  the  computer  program  scat_pulse2  in  Appendix  H.  In  a 
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similar  manner, 
displacement  may 
Equation  (117). 


the  high-frequency  limit  solution  for  the  incident  SH-wave 
be  obtained  by  substituting  the  approximation  of  Qnn  into 
That  is, 


uj(i)  = 


— I  I  jttJn(rsk2)e'^'^s 

4jrp2C2^  jl-0  m*-* 


X  I  eJn'>Jn(rk2)Hn+n(1^2rjl)  •  ^66) 

n— ® 


The  substitution  of  the  identity  given  by  Equation  (164)  into  Equation  (166) 
leads  to  a  reduced  form  of  the  high-frequency  limit  solution  for  the  direct 
SH-wave  displacement.  That  is, 


il).  — I  2bAToRoHotk2rs^^>l 
4itp2C2^  A-0 


(167) 


where: 


N  *  the  number  of  reflections  ; 


=  [tg^  +  rji^  -  2rsrji  cos  (4>Ji-1>s)l 


1/2 


((xs-H)2  +  (ys-d)2l^^^  ; 


2  1/2 

(r2  +  rji  -  2rrj^  cos  <p  ]  » 


and 


“  sln“i  sin  ♦]  • 


The  numerical  solution  of  Equation  (167)  for  the  direct  SH-wave  displacement 
is  calculated  by  the  computer  program  lnc__pulse2  in  Appendix  H. 
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Finally,  the  hlgh-f rquency  limit  approximation  of  the 
total  SH-wave  displacement  as  measured  for  a  detector  at  the  surface  of  the 
earth  is  given  by 


2  1  2bxToRojnHo[k2r8^>] 

4iip2C2‘^  1*0 


+  I  2bAToRi  I  (jitbo)Ho(k2?A)eJ'^^  .  (168) 

i.“0  ID=-® 


where  the  scattered-wave  coef '^Iclents ,  b^j,  are  obtained  from  Equation  (1A6). 

In  this  equation,  the  optical  limit  approximation  of  the  function,  Sjj,  is  given 
by  Equation  (lA8b). 


In  order  to  simplify  calculations,  more  appropriate 
expressions  for  the  distance,  rjj^  and  and  the  angle,  4»j^,  may  be  deduced 

directly  from  Figure  11-10.  These  are 


Fo  -  {Iy-d]2  +  (H-h(l-v)]2}^/^  ; 

ri  -  {[y-dl2  +  (H-h4vh(l+2il)]2}^/^  ; 

-  tan-1 

‘I'Jl  '  ^H-h^h'('H-'2f')^ 


Fg^°^  *  l(y~ys)^  (xs-h+vh)2]  ; 


and 


fgU)  -  {(y-ys)2  +  (xs-h-K.h(l+2ll)]2}  . 

6.  Numerical  Applications  using  the  Transient  Analytical  Solutions 

To  demonstrate  the  applicability  of  the  Cagnlard-deHoop  tech¬ 
nique  and  as  a  check  on  the  numerical  evaluations,  SH-wave  pulses  associated 
with  a  line  source  in  the  bedrock  are  presented  for  the  surface  layer  models 
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shown  In  Figure  II-ll.  The  source  depth  is  15  m,  the  horizontal  source- 
detector  distance  is  20  m,  the  surface-layer  thickness  is  h  =*  10  m,  and  the 
constitutive  parameters  in  the  sediment  and  the  bedrock  are  cj  =  500  m/sec, 
p  1  >■  1.5  gr/cm^,  C2  *  2,500  m/sec,  and  P2  “  2.5  gr/cm^.  The  SH-wave  pulse 
signal  transmitted  by  the  line-source  is  defined  by 


f(t)  *  pte~Pt  gjn  (2TifQ)  u(t). 


where : 


P  -  1.8  fo 


and 


fg  =  1,000  Hz. 


Figure  11-12  shows  SH-wave  pulse  responses  computed  for  the 
two-layer  earth  model  Illustrated  in  Figure  11-11 (a).  The  detector  is  in  the 
bedrock  close  to  the  surface  layer  Interface  at  a  depth  of  15  m.  The  total 
seismic  signal  shown  in  Figure  11-12  was  computed  by  the  superposition  of  the 
effects  of  the  direct  SH-wave  pulse  having  an  arrival  time  of  tg  =  8  msec  and 
three  reflected  SH-wave  pulses  having  arrival  times  of  tg  *  8.9  msec,  ti  =  48 
msec,  and  t2  *  88  msec,  respectively. 

Because  of  the  small  arrival  time  difference.  At  “  0.9  msec, 
between  the  direct  pulse  and  the  first  reflection,  these  two  waveforms  appear 
in  Figure  11-12  as  a  single  pulse  having  an  arrival  time  of  8  msec.  As 
expected,  when  the  effective  source-detector  separation  becomes  larger,  the 
number  of  multiples  in  the  surface  layer  will  Increase. 

In  order  to  investigate  the  individual  reflected  waveforms,  the 
total  SH-wave  pulse  without  the  masking  effect  of  the  direct  SH-wave  pulse  is 
illustrated  in  Figure  11-13.  Here,  the  first  and  second  reflected  pulses 
exhibit  the  same  polarity  because  of  their  similar  geometrical  ray  path.  In 
fact,  the  first  SH-wave  pulse  is  reflected  back  to  the  detector  at  the  boundary 
between  the  surface  layer  and  the  lower  half-space,  while  the  second  pulse  is 
reflected  back  to  the  detector  at  .the  air-earth  interface.  Alternatively,  the 
third  SH-wave  pulse  is  reflected  once  in  the  surface  layer  where  it  changes 
polarity  and  then  it  follows  a  trajectory  similar  to  that  of  the  previous 
pulses. 


Figure  11-14  shows  computed  SH-wave  pulses  for  a  detector  at 
the  surface  as  shown  in  Figure  Il-ll(b).  In  this  case,  the  total  field  was 
obtained  by  the  superposition  of  pulses  having  arrival  times  of  ti  =  27  msec, 
t2  “  67  msec,  and  t3  *  108  msec,  respectively.  These,  waveforms  exhibit  alter¬ 
nating  polarities.  For  example,  the  first  pulse  follows  a  positive  trajectory, 
the  second  pulse  changes  polarity  by  reflecting  once  at  the  bottom  of  the  sur¬ 
face  layer,  and  the  third  pulse  exhibits  positive  polarity  by  reflecting  twice 
the  surface. 
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FIGURE  II-IA.  S’ifNTHETIC  SEldMIC  PULSES  ASSOCIATED  WITH  THE  SH-LINE  SOURCE 
FOR  THE  BEDROCK  AND  THE  DETECTOR  AT  THE  SURFACE  OF  A 
TWO-LAYER  LOSSLESS  HALF-SPACE 
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III.  NUMERICAL  EVALUATION  OF  MODELS 


The  progressive  theoretical  analyses  presented  in  Section  II  address  SH- 
wave  scattering  from  a  cylindrical  cavity  in  the  simple  case  of  plane  waves  in 
an  unbounded  medium  to  the  more  representat'ive  cases  of  cylindrical  and  spheri¬ 
cal  SH-wave  scattering  in  a  half-space.  The  effects  of  dissipative  absorption 
in  the  two-layer  medium  were  introduced  to  yield  practical  Information  on  the 
scattered  seismic  signal  amplitudes  required  to  be  detected  and  the  various 
masking  effects  imposed  by  the  low-quality  weathered  layer  at  the  surface.  The 
special  considerations  for  performing  the  numerical  evaluations  in  each  case 
are  discussed  below. 

A.  Plane  Wave  Source  in  a  Lossless  Unbounded  Medium 

The  numerical  solution  for  plane  wave  scattering  from  a  cylindrical 
cavity  in  a  lossless  unbounded  medium  was  developed  using  Equation  (26)  pre¬ 
sented  earlier.  The  plane  wave  case  allows  an  arbitrary  unit  amplitude  to  be 
assigned  to  the  incident  wave  so  that  the  scattered  wave  values  can  be  scaled 
as  fractions  of  the  incident  wave.  By  formulating  the  frequency  dependent 
terms  in  units  of  wavenumber  multiplied  by  the  cavity  radius,  the  calculations 
were  generalized  to  any  frequency  range  or  cavity  size  by  specifying  the  fre¬ 
quency,  the  velocity,  and  the  cavity  radius  in  proper  combinations.  Programs 
were  written  to  calculate  the  energy-scattering  cross-section  and  the  SH-wave 
displacement  as  functions  of  frequency,  source  and  detector  geometry,  and 
cavity  target  distance  from  the  detector  while  holding  other  parameters  such 
as  the  SH-wave  velocity  and  mass  density  of  the  medium  constant. 


Ove rvlew  of  the  Plane  Wave  Analysis  Computer  Programs 


Several  computer  programs  were  developed  to  provide  plane  wave 
scattering  results  for  horizontally  polarized  shear  waves  in  several  informa¬ 
tive  formats.  First,  the  total  energy-scattering  cross-section  of  the  cylin¬ 
drical  cavity  target  is  computed  as  a  function  of  frequency  to  illustrate  the 
effects  on  scattering  efficiency  as  governed  by  the  incident  wavelength  and  the 
cavity  diameter.  Next,  the  scattered  SH-wave  displacement  amplitude  at  the 
surface  of  the  cavity  is  computed  as  a  function  of  angle  about  the  cavity  axis 
relative  to  the  incident  wave  direction  and  for  several  frequencies  to  funda¬ 
mentally  characterize  the  angular  scattering  signature  of  such  targets.  This 
same  form  of  calculation  is  next  carried  out  for  several  detector  distances 
away  from  the  cavity  axis  to  Illustrate  the  geometrical  spreading  losses  and 
the  far-field  angular  scattering  signatures  of  the  cavity  target.  Finally,  the 
plane-wave  analysis  program  is  adapted  to  simulate  a  plane-wave  source  and  a 
line  of  detectors  positioned  within  the  whole-space  medium  to  represent  a  seis¬ 
mic  survey  traverse.  SH-wave  displacement  amplitudes  as  might  be  detected  at 
points  along  this  survey  line  are  calculated  as  functions  of  the  cavity  depth 
below  the  traverse  line  and  for  several  source  signal  frequencies.  All  of 
these  plane-wave  computations  are  for  steady-state  source  signals  whose  wave¬ 
lengths  are  significantly  larger  than,  about  equal  to,  and  significantly 
smaller  than  the  cylindrical  cavity  diameter. 


In  each  program,  the  numerical  calculations  consist  of  summa¬ 
tions  of  terms  involving  Bessel  and  Hankel  functions  representing  the  scattered 
wave  fields  as  functions  of  the  geometrical  and  frequency  parameters.  These 
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summations  are  performed  over  at  least  15  Iterations  and  up  to  a  maximum  of 
50  Iterations  until  the  difference  In  computed  values  between  two  successive 
Iterations  Is  less  than  0.01  percent  of  the  calculated  value.  In  all  of  the 
computer  programs  except  that  for  the  total  energy  scattering  cross-section, 
this  summation  Is  performed  by  the  subroutine  called  plane_scatter.  The  com¬ 
puter  codes  for  the  plane-wave  analysis  programs  and  subroutines  are  written 
In  FORTRAN  77  language  and  are  presented  In  Appendix  H,  Section  H.l. 

a.  Energy  Scattering  Cross-Section  Calculation 

The  program  called  csect  computes  the  total  energy  scat¬ 
tering  cross-section  of  a  cylindrical  cavity  when  Illuminated  broadside  by  a 
plane  SH-wave.  This  cross-section,  discussed  earlier  In  Section  II. B. 3,  Is 
calculated  for  values  of  wavenumber  times  cavity  radius  (l.e.,  ka)  ranging 
from  0.05  to  2.0  In  Increments  of  0.05. 

b.  SH-Wave  Displacement  Amplitude 

Four  programs  are  used  to  calculate  the  displacement 
amplitudes  of  plane  SH-waves  scattered  from  a  cylindrical  cavity  In  a  homo¬ 
genous  whole-space. 


The  first  program,  dlspa,  calculates  the  amplitude  of  the 
total  displacement  field  as  a  function  of  frequency  at  three  different  posi¬ 
tions  on  the  surface  of  the  cavity.  The  frequency  Is  contained  In  the  factor, 
ka,  for  values  ranging  from  0  to  2  In  steps  of  0.05  and  from  2.1  to  10  In  steps 
of  0.1.  The  locations  at  the  surface  of  the  cavity  at  which  the  calculations 
are  evaluated  are  on  radii  at  45,  90,  and  135  degrees  relative  to  the  incident 
wave  vector.  The  displacement  amplitude  is  a  complex  quantity  expressed  In 
terms  of  amplitude  and  phase,  which  is  unwrapped. 

The  second  program,  dlspb,  calculates  the  total  displace¬ 
ment  amplitude  as  a  function  of  distance  from  the  cavity,  along  three  different 
radial  lines  and  at  three  different  frequencies.  The  distance  is  expressed  in 
terms  of  the  radial  distance  divided  by  the  cavity  radius  (l.e.,  r/a)  calcu¬ 
lated  from  values  of  1  to  10  in  steps  of  0.5.  The  angles  of  the  radial  lines 
along  which  the  displacements  are  calculated  are  45,  90,  and  135  degrees  with 
respect  to  the  Incident  wave  vector.  The  frequencies  at  which  the  results  are 
obtained  are  expressed  in  terms  of  the  factor,  ka,  and  are  0.5,  1,  and  2.  The 
displacement  amplitude  Is  a  complex  quantity,  but  only  the  magnitude  of  this 
quantity  Is  used  as  output. 

The  third  program,  dispcd,  calculates  the  amplitude  of 
the  scattered  SH-wave  displacement  field  as  functions  of  angle,  detector  dis¬ 
tance  from  the  cavity,  and  frequency.  Depending  upon  whether  a  certain  sec¬ 
tion  of  code  noted  In  the  program  listing  is  commented  out  or  not,  the  program 
will  provide  either  the  scattered  amplitude  only  or  the  combined  direct  and 
scattered  amplitude.  The  scattered  amplitude  Is  calculated  at  angles  from  0 
to  180  degrees  In  steps  of  1  degree.  The  frequencies  are  expressed  in  terms 
of  the  factor,  ka,  for  values  of  0.5,  1,  and  2.  The  radial  distances  from  the 
cavity  axis.  In  units  of  r/a,  are:  1,  2,  5,  10,  20,  50,  and  100.  The  dis¬ 
placement  amplitude  is  a  complex  quantity,  but  only  the  magnitude  of  this 
quantity  Is  used  as  output. 
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The  fourth  program,  dlspe,  calculates  the  amplitude  of 
the  displacement  field  at  a  number  of  detector  locations  on  the  x-axis  for  a 
cavity  located  on  the  y-axis  at  several  distances  below  the  x-axis.  The  inci¬ 
dent  wave  propagation  vector  may  be  directed  at  an  angle  with  respect  to  the 
y-axis  (i.e.,  vertical),  and  may  vary  from  0  to  180  degrees.  The  frequencies 
at  which  the  displacement  is  calculated  are  expressed  in  terms  of  the  factor, 
ka,  and  are  0.5,  1,  and  2.  The  distances  of  the  tunnel  origin  from  the  x-  and 
y-axls  origin  are  entered  in  terms  of  r/a  and  are  20,  40,  60,  80,  and  100.  The 
positions  on  the  x-axls  at  which  the  displacement  are  calculated  range  from  -50 
to  +50  meters  in  steps  of  0.5  meters.  The  angle  between  the  incident  propa¬ 
gation  vector  and  the  y-axis  may  be  either  preset  or  calculated,  according  to 
which  segments  of  the  program  code  are  commented  out  or  not.  If  the  preset 
version  is  chosen,  then  either  backscatter  or  through  transmission  results  may 
be  produced.  If  the  calculated  version  is  chosen,  the  propagation  vector  is 
specified  to  lie  on  a  line  with  a  y-axis  intercept  of  the  cavity  position  and 
an  x-axis  Intercept  of  -54  meters.  This  orientation  of  the  incident  wave  rela¬ 
tive  to  the  cavity  and  the  x-axlff  is  similar  to  actual  practice  in  seismic 
exploration  where  the  location  of  the  source  point  is  at  the  end  of  the  detec¬ 
tor  line.  The  displacement  amplitude  is  a  complex  quantity,  but  only  the  mag¬ 
nitude  of  this  quantity  is  used  as  output. 

B .  Line  Source  in  a  Half-Space 


The  numerical  evaluation  of  the  case  of  a  line  source  in  a  half¬ 
space  was  incrementally  developed  beginning  with  the  lossless  unbounded  medium 
case,  followed  by  the  addition  of  a  half-space  boundary,  a  surface  layer, 
lossy  media  factors,  and  finally  in  a  form  capable  of  providing  calculated 
results  for  any  combination  of  these  conditions.  For  example,  by  properly 
selecting  the  input  parameters,  the  program  can  be  made  to  provide  results 
for,  say,  an  unbounded  lossy  medium  or  a  two-layered  lossy  medium,  to  cite  two 
cases.  The  cases  of  a  lossless  and  a  lossy  simple  half-space  and  the  case  of 
a  two-layered  lossy  half-space  where  the  source  is  located  in  the  lower  layer 
and  the  detector  is  located  either  in  the  lower  layer  or  at  the  upper  boundary 
of  the  surface  layer  are  considered  in  detail. 

1 .  Overview  of  the  Program 

The  program,  synsels,  presented  in  Appendix  H,  Section  H.2.1, 
is  the  realization  of  the  numerical  solution  of  the  equations  for  scattering 
of  a  cylindrical  SH  seismic  wave  from  a  cylindrical  cavity  derived  in  Section 
II.  Since  the  direct  and  scattered  signals  received  at  the  detector  are 
decoupled  in  the  approximate  solution  developed  for  deep  cavities,  the  program 
solves  for  each  signal  component  separately. 

The  basic  solution  method  used  for  each  signal  component  is  the 
same.  The  time  independent  complex  amplitude  of  the  wave  is  found  as  a  func¬ 
tion  of  frequency  and  then  Inverse  Fourier  transformed  into  the  time  domain, 
where  it  represents  the  response  generated  by  a  delta  function  pulse  emanating 
from  the  line  source.  This  response  contains  time  delay  information  related  to 
the  propagation  of  the  wave  from  the  source  along  the  various  propagation  paths 
to  the  detector.  This  delta  function  response  is  then  convolved  with  the 
source  pulse  waveform  to  obtain  the  composite  signal  at  a  given  detector  posi¬ 
tion.  The  same  general  method  is  applied  for  both  the  direct  wave  and  the 
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scattered  wave  at  each  detector  position,  and  the  results  are  written  out  to  a 
file. 


Sample  input  and  output  files  for  program,  synseis,  are  pre¬ 
sented  in  Appendix  H,  Section  H.2.1. 

a.  Direct  Wave  Calculation 


The  direct  wave  signal  is  calculated  by  the  subroutines 
called  inc_j)ulse  and  inc_pulse2.  These  two  subroutines  are  essentially  iden¬ 
tical.  The  inc_pulse  subroutine  is  written  to  provide  a  solution  for  the  case 
for  the  source  and  detector  located  in  the  lower  layer,  and  the  lnc_pulse2 
subroutine  is  written  to  solve  the  case  for  the  source  in  the  lower  layer  and 
the  detector  at  the  top  of  the  surface  layer.  As  a  result,  the  description 
that  follows  is  applicable  to  both  subroutines. 

The  subroutine  is  given  the  source  and  detector  posi¬ 
tions  as  well  as  the  values  of  the  other  parameters  to  be  used  in  the  calcu¬ 
lation.  For  zero  to  n  reflections  from  any  layer  Interfaces  that  might  be 
present,  the  time-independent  direct  wave  amplitude  along  a  given  propagation 
path  is  calculated  at  128  wavenumber  values  from  0  to  12.7  m~^  in  steps  of  0.1 
m“^.  The  time  delay  for  the  arrival  of  the  wave  is  calculated  from  the  geome¬ 
try  and  the  velocity  in  the  medium  and  subtracted  from  the  frequency  domain 
signal  by  constant  phase  subtraction  to  prevent  aliasing  of  the  phase.  A 
cosine  window  function  is  then  applied  to  the  frequency  domain  signal  for  a 
wavenumber  range  of  7.8  to  12.8  m~^  to  prevent  spurious  ringing  in  the  time 
domain  signal.  The  frequency  domain  signal  is  then  inverse  Fourier  transformed 
to  obtain  a  256-polnt  time  domain  waveform.  This  waveform  is  then  convolved  in 
the  time  domain  with  the  256-point  source  pulse  waveform  to  produce  the  signal 
present  at  the  detection  point.  This  signal  is  then  added,  often  being  offset 
by  the  initial  time  delay  with  respect  to  the  total  waveform,  and  the  process 
is  repeated  for  all  reflections  that  are  to  be  included  in  the  calculation. 

b.  Scattered  Wave  Calculation 


The  scattered  wave  signal  is  calculated  in  the  sub¬ 
routines  called  scat_j)ulse  and  scat_pulse2.  These  two  subroutines  are  essen¬ 
tially  identical.  The  scat_pulse  subroutine  is  written  to  provide  a  solution 
for  the  case  for  the  source  and  detector  located  in  the  lower  layer,  and  the 
scat_pulse2  subroutine  is  written  to  solve  the  case  for  the  source  in  the 
lower  layer  and  the  detector  located  at  the  top  of  the  surface  layer.  As  a 
result,  the  description  that  follows  is  applicable  to  both  subroutines. 

The  subroutine  is  provided  with  the  source,  detector,  and 
cavity  locations  as  well  as  other  relevant  parameters.  The  algorithm  makes 
repeated  use  of  a  number  of  Hankel  function  values,  so  these  are  precomputed. 
The  routine  then  begins  an  Iteration  over  the  number  of  direct  and  reflected 
terms  that  are  selected  for  the  calculation.  For  a  given  propagation  path, 
the  time-independent  scattered  wave  signal  as  a  function  of  frequency  is  cal¬ 
culated  at  128  points  at  wavenumber  values  from  0  to  12.7  m“^.  At  each  of 
these  frequencies,  a  summation  is  performed  to  obtain  the  complex  amplitude. 

A  minimum  of  ten  terms  of  summation  are  calculated,  with  up  to  fifty  terms 
possible,  until  the  accuracy  of  the  sum  in  two  consecutive  terms  is  better 
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than  the  accuracy  limit  specified  in  the  program,  which  is  1  x  10”^.  After 
the  signal  at  all  128  frequencies  have  been  calculated,  the  time  delay  is  sub" 
tracted  and  the  cosine  window  is  applied  as  in  the  direct  wave  case.  The 
frequency  domain  signal  is  then  inverse  Fourier  transformed  and  convolved  with 
the  source  pulse  waveform.  This  result  is  then  added,  with  its  adjusted  time 
offset,  to  the  total  waveform,  and  the  process  is  repeated  for  all  selected 
propagation  paths. 

2.  Line  Source  in  a  Homogeneous  Half-Space 

For  the  case  of  a  homogeneous  half-space,  the  program  is  given 
an  input  file  in  which  the  density  of  layer  1  (the  surface  layer)  is  0  and  the 
number  of  reflections  allowed  is  1.  The  depths  of  the  line  source,  the  cylin¬ 
drical  cavity,  and  the  detectors  specified  in  the  input  file  are  referenced  to 
the  depth  of  the  second  layer.  If  desired,  the  depth  of  the  first  layer  may 
be  set  to  zero  to  simplify  the  task  of  depth  calculation.  If  the  surface 
layer  depth  is  set  to  zero,  care  must  be  taken  not  to  set  the  detector  depth 
to  exact  zero  as  well,  since  this  may  cause  the  program  to  produce  erroneous 
results. 


a.  Lossless  Medium 


To  generate  a  solution  for  a  lossless  half-space  medium, 
the  attenuation  cutoff  frequency  for  the  second  layer  should  be  set  to  a  very 

large  value,  say  1.0  x  10^®.  This  will  prevent  onset  of  attenuation  and  dis¬ 

persive  effects  in  the  practical  range  of  seismic  range  frequency.  The  quality 
factor  for  the  second  layer  may  also  be  set  to  a  high  value,  if  desired,  to 

suppress  any  dispersive  effects  in  case  the  cutoff  frequency  was  not  set  high 

enough. 


b.  Lossy  Medium 

To  obtain  a  solution  for  a  lossy  homogeneous  half-space, 
the  cutoff  frequency  and  quality  factor  should  be  set  so  as  to  simulate  the 
ground  conditions  desired.  It  should  be  kept  in  mind  that  the  cutoff  frequency 
is  given  in  terms  of  angular  frequency,  u.  To  simulate  a  lossy  medium  over  all 
frequencies  of  interest,  set  the  cutoff  frequency  to  a  value  less  than  the  low¬ 
est  frequency  of  interest.  When  the  cutoff  frequency  is  set  to  be  smaller  than 
the  value  of  the  frequency  step,  the  attenuation  condition  defined  by  the  Q 
factor  applies  to  all  frequencies.  The  quality  factor  should  be  set  to  larger 
values  for  less  lossy  media  and  to  smaller  values  for  more  lossy  media. 

3.  Line  Source  in  a  Two-Layer  Lossy  Half-Space 

To  simulate  a  line  source  in  a  two-layer  lossy  half-space, 
non-trivial  values  should  be  given  to  all  parameters  in  the  input  file.  If 
only  one  layer  is  to  be  lossy,  the  other  layer  may  be  given  a  large  cutoff 
frequency  or  quality  factor,  as  discussed  earlier.  The  layer  depth  parameter 
should  be  non-zero  in  any  case. 

a.  Source  and  Detector  in  Lower  Layer 

For  the  case  in  which  the  source  and  detectors  are  both 
located  in  the  lower  layer,  it  is  only  necessary  to  set  the  detector  array 
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y-orlgin  to  a  negative  value  greater  than  the  layer  depth.  If  the  yorlgin 
Is  less  than  the  layer  depth,  unpredictable  results  will  occur. 

b.  Source  in  Lower  Layer  and  Detector  in  Surface  Layer 

For  the  case  in  which  the  detectors  are  located  at  the 
top  of  the  surface  layer,  the  y-orlgin  should  be  set  to  zero.  The  present 
version  of  the  program  will  not  properly  treat  cases  in  which  the  detectors 
are  located  at  other  depths  in  the  surface  layer. 

C.  Graphing  Programs 

The  graphics  programs  to  produce  output  plots  of  the  results  from 
the  above  programs  are  presented  in  Appendix  H,  Sections  H.3.1  and  H.3.2. 
These  programs  are  written  in  C  language  for  ease  of  interfacing  with  the 
graphics  package,  Starbase,  on  the  Hewlett-Packard  9000  model  520  computer 
equipped  with  the  HP-UX  operating  system. 

The  graphics  program  for  the  display  of  synthetic  seismograms, 
called  vusels,  utilizes  an  input  file  to  specify  the  presentation  of  the  data. 
A  sample  input  file  is  included  in  Appendix  H,  Section  H.3.3,  with  a  short 
explanation  of  the  options  available. 
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IV.  MODEL  CALCULATIONS  AND  RESULTS 


A.  Plane  SH  Waves  Scattered  from  a  Cylindrical  Cavity  In  a  Lossless 
Unbounded  Medium 


The  theoretical  analysis  of  plane  SH-wave  scattering  from  a  cylin¬ 
drical  cavity  representative  of  a  horizontally  oriented  tunnel  developed  in 
Section  II. B  is  used  to  provide  preliminary  Insights  into  the  problem  of 
reflection  seismic  detection  of  such  tunnel  targets  at  depths  of  about  100 
meters  or  more.  In  the  results  to  follow,  simple  but  useful  and  informative 
measures  of  the  SH-wave  scattering  target  strength  and  directional  scattering 
signature  of  cylindrical  cavities  are  presented  for  the  conditions  of  geologi¬ 
cal  host  medium  and  tunnel  size  of  primary  interest.  The  numerical  model  is 
then  used  to  simulate  the  steady-state  detection  responses  that  would  be 
observed  along  a  surface  seismic  reflection  traverse  and  in  a  hole-to-hole 
seismic  scan. 

1 .  Energy  Scattering  Cross-Section 

Figure  IV-1  shows  the  geometry  used  to  analyze  Che  scattering 
of  a  plane  SH  wave  incident  on  a  cylindrical  cavity  in  an  unbounded  medium. 

The  total  SH-wave  displacement  at  measurement  point,  (x,  y),  consists  of  the 
sum  of  Che  incident  wave  and  the  wave  scattered  from  the  cylindrical  target. 


<j)  =  TAN-^^ 

X 

r  =  (x^  +  y^)^^ 


FIGURE  IV-1.  GEOMETRY  OF  PLANE  SH-WAVE  INCIDENT  ON  A 
CYLINDRICAL  CAVITY  IN  AN  UNBOUNDED  MEDIUM 


This  total  displacement  is  expressed  by 
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The  total  elastic  wave  energy  scattered  by  the  cylindrical  cavity  can  be  com¬ 
puted  from  the  second  term  of  Equation  (IV-1)  by  summing  the  products  of  each 
complex  displacement  term  times  its  complex  conjugate.  This  scattered  SH-wave 
energy  computed  for  r“a  may  be  divided  by  the  physical  (broadside)  cross- 
sectional  area  of  one  unit  length  of  the  cavity  to  yield  the  normalized  SH-wave 
energy  scattering  cross  section  for  an  incident  plane  SH  wave.  Thus,  the  SH- 
wave  scattering  efficiency  may  be  expressed  by  the  apparent  energy  scattering 
cross  section  of  the  cylinder  target  as 


^s 


—  r  A  a' 


n 


I  AnA^] 
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(IV-2) 


where: 


^n-{ 


1;  for  n  ■  1 
2;  for  n  >  1  ; 


(i) 


k  *  —  ;  and 


ka  ■  ^  ;  D  ■  cavity  diameter. 


Figure  IV-2  shows  the  normalized  SH  wave  energy  scattering 
cross  section  for  a  2-meter  diameter  cylindrical  cavity  in  a  medium  having  a 
shear-wave  velocity  of  Vg  “  2,500  m/s. 


FIGURE  IV-2.  NORMALIZED  ENERGY  SCATTERING  CROSS-SECTION  OF  A 
CYLINDRICAL  CAVITY  FOR  AN  INCIDENT  PLANE  SH  WAVE 
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In  Fir:re  IV-I,  the  parameter  ka  is  proportional  to  frequency  and,  for  the  value 
of  ka  =  t ,  the  frequency  is  1,250  Hz  and  the  shear  wavelength  is  equal  to  the 
diameter  of  the  cylindrical  cavity.  The  normalized  cross  section  in  this  case 
is  approximately  1.55,  a  value  which  is  about  55  percent  larger  than  the  area, 

D,  corraspondimg  to  the  per  unit  length  cross-sectional  area  of  the  tunnel. 

2.  Polar  Amplitude  Scattering  Patterns 

Figure  IV-3  shows  the  SH-wave  amplitude  scattered  from  a  2-meter 
diameter  cylindrical  cavity  when  the  detector  point,  (x,  y)  is  at  a  radius  of 
100  meters  away  from  the  cavity  axis. 

r/a»  100.  ka  »  .5 - 

a  =  I  METER  ka  =  1.0 - 

V  =  2,500  m/s  ka  *=  2.0 - 

s 
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(a)  Total  SH-Wave  Polar  Amplitude 
98 


(b)  SH-Wave  Reflected  Polar  Amplitude  Component 
FIGURE  rV-3.  PLANE  SH-WAVE  POLAR  AMPLITUDE  SCATTERING  FROM  A  CYLINDRICAL  CAVITY 
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Figure  IV-3(a)  shows  the  total  detected  amplitude  consisting 
of  the  combined  incident  and  scattered  waves  as  observed  at  a  distance  of 
100  meters  away  from  the  tunnel.  Figure  lV-3(b)  shows  the  scattered  component 
only,  Illustrating  the  directional  SH-wave  scattering  pattern  of  the  cylindri¬ 
cal  cavity  and  the  fact  that  the  scattered  amplitude  is  in  the  range  of  5-7.5 
percent  of  the  Incident  plane  wave  amplitude  at  point  (x,  y)  for  angular 
directions  of  “  180  ±30  degrees.  The  values  of  ka  Illustrated  in  these  fig¬ 
ures  correspond  to  frequencies  of  approximately  200  Hz,  400  Hz,  and  800  Hz  for 
the  condition  where  Vs  “  2,500  m/sec.  Frequencies  above  800  Hz  will  exhbiit 
backscatter  amplitudes  up  to  about  10  percent  of  the  incident  wave  amplitude. 

3.  Seismic  Survey  Traverse  Simulations 

Figure  IV-4  illustrates  the  amplitude  profile  of  the  total 
plane  wave  scattered  signal  projected  onto  a  representative  seismic  detector 
line  traverse  for  several  cavity  depths.  For  the  three  SH-wave  frequencies 
used  in  the  analysis,  the  Interference  between  the  source  wave  and  scattered 
wave  component  is  apparent,  with  more  spatial  fluctuations  in  the  near-source 
half  of  the  detector  spread.  The  directional  scattering  aspects  of  the  cylin¬ 
drical  target  are  only  slightly  noticeable  in  the  shallowest  target  depth  case 
but  can  be  anticipated  to  be  more  prevalent  at  higher  frequencies. 

These  preliminary  theoretical  results  show  that  the  scattered 
amplitude  of  SH  waves  at  distances  up  to  100  meters  away  from  the  cylindrical 
cavity  have  practical  detection  amplitudes.  The  magnitude  of  the  scattered 
waves  diminish  as  the  inverse  square  root  of  the  distance.  For  the  typical 
shear  wave  velocity  of  Vg  ■  2,500  m/s,  spectral  components  above  about  400  Hz 
are  strongly  reflected.  The  polar  amplitude  scattering  patterns  indicate  that 
the  cylindrical  cavity  has  a  directivity  which  may  serve  as  a  useful  means  for 
identifying  such  scattering  targets. 

4.  Hole-to-Hole  SH-Wave  Scattering  Simulations 

By  using  a  different  choice  of  the  detector  locations  used  in 
the  seismic  traverse  simulations,  the  plane-wave  analysis  program  may  be  used 
to  simulate  the  steady-state  response  that  would  be  observed  in  seismic  scan¬ 
ning  with  the  tunnel  cavity  oriented  transversely  between  the  boreholes. 

Figure  IV-5  shows  model  calculations  for  horizontal  ray  path  hole-to-hole  SH- 
wave  scattering  from  a  2-meter  diameter  tunnel  as  observed  when  the  cavity  is 
located  5,  10,  15,  and  25  meters  away  (toward  the  source)  from  the  borehole 
containing  the  detector. 

The  scattered  SH-wave  displacement  amplitude  is  symmetrical 
about  the  horizontal  ray  which  passes  through  the  center  of  the  cylindrical 
cavity,  and  the  scattered  signal  profile  has  the  characteristic  "W"  shape 
which  corresponds  to  the  diffraction  interference  of  waves  scattered  by  a 
symmetrical  obstacle.  The  frequency  dependence  of  the  forward-scattered  SH 
waves  for  values  of  ka  of  0.5,  1.0,  and  2.0  are  only  very  slightly  stronger 
IV-4  for  the  reflection  simulation  case. 

B.  Synthetic  Seismograms  of  SH  Waves  Scattered  from  a  Cylindrical 

Tunnel  in  a  Lossy  Unbounded  Medium 

The  numerical  method  developed  in  Section  II. D  has  been  implemented 
to  compute  synthetic  seismograms  of  SH-wave  scattering  from  a  cylindrical 
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TO  HOLE  RESPONSES  ,  ,  SH-WRVE  HOLE  TO  HOLE  RESPONSES 


FIGURE  IV-5.  TOTAL  PLANE  SH-WAVE  AMPLITUDE  PROFILES  BETWEEN  BOREHOLES 


cavity  in  a  lossless  medium  and  in  a  lossy  medium  for  incident  cylindrical 
SH-waves.  Tlie  model  geometry  and  the  parameters  of  the  medium  are  given  in 
Figure  IV-6.  The  characteristics  of  the  SH-wave  pulse  signal  transmitted  by 
the  line  source  are  presented  in  Sections  II. D. 6  and  II. F. 6. 


(x,y) 


4>  =  tan  ’  - 

X 

<t>c  =  tan"’ 

-X, 


fo  =  1000  Hz 
Cj  =  2500  m/sec 
Q  =  100 


FIGURE  IV-6.  GEOMETRY  OF  CYLINDRICAL  WAVE  LINE  SOURCE  AND 
CYLINDRICAL  CAVITY  IN  AN  UNBOUNDED  MEDIUM 


Figures  IV-7  through  IV-13  illustrate  synthetic  seismograms 
calculated  for  the  SH-vave  source  and  a  2-meter  diameter  cylindrical  cavity  in 
a  medium  having  a  shear-wave  velocity  of  Cg  =  2500  m/s.  The  SH-llne  source  is 
specified  to  produce  a  pulse  having  a  peak  frequency  fg  =  1000  Hz  and  is 
located  at  the  source  position  (-52m, 0).  The  applied  SH-wave  pulse  signal 
given  by  Equation  (63)  and  its  frequency  spectrum  given  by  Equation  (6A)  are 
illustrated  in  Figure  IV-IA  (a)  and  (b),  respectively.  The  seismic  detector 
line,  having  detector  spacings  of  2  m,  is  transverse  to  the  axis  of  the  cylin¬ 
drical  cavity.  The  synthetic  seismograms  were  calculated  for  both  a  perfectly 
elastic  medium  and  a  lossy  medium.  The  seismograms  associated  with  the  lossy 
medium  were  computed  for  a  shear-wave  quality  factor  of  Q  =  100,  which  cor¬ 
responds  to  a  spatial  attenuation  rate  of  0.25  dB/wavelength. 

Figure  IV-7  shows  a  synthetic  seismogram  containing  the  direct 
SH  wave  and  the  scattered  SH  wave  from  a  tunnel  target  located  at  a  depth  of 
50  m  in  a  lossless  medium.  As  expected,  the  time  of  arrivals  and  waveforms 
essoclated  with  the  direct-traveling  SH  wave  are  detected  first,  and  the  scat¬ 
tered  SH  waves  are  detected  later  in  the  time  section.  The  tunnel  target 
reflections  form  a  pattern  having  a  hyperbolic  shape  which  is  characteristic 
of  reflections  from  the  localized  air-filled  tunnel  target. 

Figures  lV-8  and  IV-9  show  only  the  direct-traveling  SH-wave 
seismic  pulses  radiated  by  the  SH-line  source  in  a  lossless  and  in  a  lossy 
medium,  respectively.  These  results  indicate  the  degree  to  which  the  ampli¬ 
tudes  of  the  SH-wave  pulses  are  reduced  by  a  medium  having  an  attenuation  rate 
of  0.25  dB  per  wavelength. 
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SH~WRVE  SYNTHETIC  SEISMOGRRM 


FIGURE  IV-7.  SYNTHETIC  SEISMOGRAM  OF  THE  COMBINED  DIRECT  AND  REFLECTED  SH-WAVE  PARTICLE  DISPLACEMENT  FOR  A  2-i 
DIAMETER  CYLINDRICAL  CAVITY  AND  AN  SH-LINE  SOURCE  IN  A  LOSSLESS  MEDIUM  HAVING  AN  SH  VELOCITY  OF  2,500  m/s. 
THE  CAVITY  IS  LOCATED  AT  A  DEPTH  OF  50  m  BELOW  THE  SURFACE,  AND  THE  SH-LINE  SOURCE  IS  LOCATED  AT  THE 
POSITION  (-52  m,  0).  THE  TRACES  CORRESPOND  TO  THE  OUTPUTS  OF  DETECTORS  SPACED  2  m  APART. 
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SH-WflVE  SYNTHETIC  SEISMOGRRM 


FIGURE  IV-9.  SYNTHETIC  SEISMOGRAM  OF  DIRECT-TRAVELING  SH-WAVES  FROM  AN  SH-LINE  SOURCE  LOCATED  IN  A  LOSSY 
MEDIUM  AT  THE  SOURCE  POSITION  (-52  m,  0).  THE  SH-WAVE  VELOCITY  AND  THE  QUALITY  FACTOR  IN  THE  MEDIUM 
ARE  c  »  2,500  m/s  AND  Q  -  100,  RESPECTIVELY.  THE  DETECTOR  SPACING  IS  2  m.  THE  AMPLITUDE 
PLOTTING  SCALE  USED  IN  THESE  TRACES  IS  THE  SAME  AS  THAT  IN  FIGURE  IV-8. 
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FIGURE  IV-11.  SYNTHETIC  SEISMOGRAMS  OF  SCATTERED  SH-WAVES  ONLY  FROM  A  2-m  DIAMETER  CYLINDRICAL  CAVITY 
LOCATED  IN  A  LOSSY  MEDIUM  AT  A  DEPTH  OF  50  m.  THE  SH-WAVE  VELOCITY  AND  THE  QUALITY  FACTOR  OF  THE 
MEDIUM  ARE  c  -  2,500  m/s  AND  Q  -  100,  RESPECTIVELY.  THE  DETECTOR  SPACING  IS  2  m. 
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FIGURE  IV-13.  SYNTHETIC  SEISMOGRAMS  OF  SCATTERED  SH-WAVES  ONLY  FROM  A  2-m  DIAMETER  CYLINDRICAL  CAVITY 
LOCATED  IN  A  LOSSY  MEDIUM  AT  A  DEPTH  OF  100  m.  THE  SH-WAVE  VELOCITY  AND  THE  QUALITY  FACTOR  OF 
THE  MEDIUM  ARE  c  =  2,500  m/s  AND  Q  »  100,  RESPECTIVELY.  THE  DETECTOR  SPACING  IS  2  m. 


Figures  IV-10,  IV-ll(a),  and  IV-ll(b)  show  only  the  SH  waves 
scattered  from  the  cylindrical  cavity  located  at  a  depth  of  50  m  in  a  lossless 
medium  and  in  a  lossy  medium.  These  results  show  that  the  tunnel  reflection 
signature  in  the  lossy  medium  is  significantly  reduced  with  respect  to  that 
from  a  tunnel  target  in  an  ideal  elastic  medium.  Figure  IV-ll(b)  shows  the 
hyperbolic  reflected  SH-wave  pattern  plotted  with  an  increased  scale  factor  of 
800,  which  is  about  four  times  greater  than  the  scale  factor  used  to  display 
the  attenuated  synthetic  seismogram  Illustrated  in  Figure  IV-ll(a)  (l.e.,  the 
effective  increase  in  display  sensitivity  is  a  factor  of  200). 

Figures  IV-12,  IV-13(a),  and  IV-13(b)  illustrate  the  synthetic 
seismograms  of  scattered  SH-waves  from  a  2-m  diameter  cylindrical  cavity 
located  at  a  depth  of  100  •■.  As  expected,  when  the  depth  of  the  cavity  is 
Increased,  the  amplitude  of  scattered  SH-wave  reflections  decrease  because  of 
the  geometric  spreading  loss  as  well  as  because  of  the  excess  absorptive 
attenuation  as  shown  in  Figure  IV-13(a).  To  adequately  display  the  tunnel 
reflection  signatures  for  these  conditions,  an  effective  scale  factor  of  200 
relative  to  that  in  Figure  lV-13(a)  is  used  to  show  the  synthetic  seismogram 
for  the  100-m  deep  tunnel  in  Figure  IV-13(b). 

C.  Synthetic  Seimograms  of  SH  Waves  Scattered  from  a  Cylindrical  Cavity 

in  a  Two-Layered  Lossy  Half-Space 

The  numerical  method  developed  in  Section  II. D  is  used  to  compute 
synthetic  seismograms  of  SH  waves  scattered  from  a  cylindrical  tunnel  in  a  two¬ 
layered  lossy  half-space.  These  seismograms  have  been  computed  for  the  cases 
of  a  line  of  detectors  in  the  bedrock  and  a  line  of  detectors  at  the  top  of  the 
surface  layer.  This  model  represents  a  systematic  display  of  synthetic  seismo¬ 
grams  for  several  cases  in  which  the  surface-layer  thickness,  quality  factor,  . 
and  shear-wave  velocity  together  with  tunnel  target  depth  and  the  detector 
location  (either  burled  or  surface)  are  independent  parameters. 

The  model  geometry  is  shown  in  Figure  IV-15,  and  the  various  parame¬ 
ters  of  the  calculated  cases  are  listed  in  Table  IV-1.  All  of  the  synthetic 
seismograms  are  computed  for  the  SH-line  source  located  at  the  position  [-52m, 
-(h+1)]  and  for  a  2-m  diameter  cylindrical  tunnel  in  a  medium  having  a  shear- 
wave  velocity  of  C2  *  2500  m/s.  The  seismic  detector  line  having  detector 
spacings  of  2  m  is  transverse  to  the  axis  of  the  cylindrical  cavity.  The  char¬ 
acteristics  of  the  SH-llne  pulse  signal  transmitted  by  the  line  source  and  the 
number  of  detectors,  nd,  are  the  same  as  those  used  earlier  in  the  unbounded 
medium  case.  Appendix  G  presents  the  practical  case  in  which  the  ground  over- 
lying  the  tunnel  target  consists  of  a  weathered  surface  which  is  typically  more 
lossy  than  the  deeper  and  more  competent  medium  containing  the  tunnel  cavity. 

The  model  calculations  for  this  case  are  aimed  at  Illustrating  and  comparing 
the  tunnel  detection  results  obtained  when  the  detectors  are  buried  (coupled  in 
boreholes)  at  the  depth  of  the  competent  (lower  loss)  ground  with  those  obtained 
when  the  detectors  are  coupled  to  the  top  of  the  less  competent  ground  surface. 

Tables  IV-2,  -3,  and  -A  summarize  parameters  used  in  the  two-layer 
model  calculations  and  identify  the  corresponding  figure  numbers  of  the  respec¬ 
tive  synthetic  seismograms  Illustrated  in  Appendix  G.  In  the  ilustrations  of 
Appendix  G,  the  figures  are  subdivided  in  three  parts  identified  as  (a),  (b), 
and  (c).  The  Illustrations  identified  by  (a)  show  the  total  SH-wave  particle 
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y  >  0 
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FIGURE  IV-15.  GEOMETRY  AND  MODEL  PARAMETERS  FOR  AN  SH-LINE  SOURCE 
AND  A  CYLINDRICAL  CAVITY  IN  A  TWO-LAYER  HALF-SPACE 

(a)  Detectors  In  the  bedrock 

(b)  Detectors  at  the  surface 
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TABLE  IV- 1 


LIST  OF  PARAMETERS  FOR  THE  TWO-LAYER  LOSS  HALF-SPACE 
(Ref.  Figure  IV-14) 

(xg,yg)  =  SH-llne  source  coordinates 

(x,y)  *  Detector  coordinates 

Ax  “  Detector  spacing 

nd  *  Number  of  detectors 

2a  *  Cavity  diameter 

H  =  Cavity  depth 

h  =  Layer  thickness 

N  *  Number  of  reflections 

Cl  *  Surface  layer  SH-wave  velocity 

C2  “  Bedrock  SH-wave  velocity 

PI  *  Surface  layer  density 

P2  *  Bedrock  density 

Ql  *  Surface  layer  quality  factor 

Q2  *  Bedrock  quality  factor 
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TABLE  IV-2 


TWO-LAYER  MODEL  WITH  DETECTORS  IN  THE  BEDROCK 
Ys  “  y  “  "(h+l)  m 


(G-la,  G-lb,  G-lc) 
(G-2a,  G-2b,  G-2c) 
(G-3a,  G-3b,  G-3c) 
(G-4a,  G-4b,  G-4c) 


Tunnel  Depth  (H)  =  50  m 


h 

(m) 

Ql 

Cl 

(m/sec) 

C2 

(m/sec) 

PI  , 

(kg  /m-^) 

P2- 

(kg/m^) 

5 

50 

100 

200 

2500 

1500 

2700 

10 

50 

100 

200 

2500 

1500 

2700 

5 

10 

100 

200 

2500 

1500 

2700 

10 

10 

100 

200 

2500 

1500 

2700 

Tunne 1 


)th  (H)  =  100  m 


(G-5a, 

G-5b, 

G-5c) 

5 

50 

100 

200 

2500 

1500 

2700 

(G-6a. 

G-6b, 

G-6c) 

10 

50 

100 

200 

2500 

1500 

2700 

(G-7a, 

G-7b, 

G-7c) 

5 

10 

100 

200 

2500 

1500 

2700 

(G-8a, 

G-8b, 

G-8c) 

10 

10 

100 

200 

2500 

1500 

2700 

Xg  =  -52  m,  Ax  =  2  m,  nd  *  51,  2a  =  2  m,  and  N  =•  3. 
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TABLE  IV-3 


TWO-LAYER  MODEL  WITH  DETECTORS  AT  SURFACE 

yg  =  -(h+1)  m 
y  -  0 


Tunnel  Depth  (H)  »  50  m 


Figure  Numbers 

h 

(m) 

Ql 

Q2_ 

Cl 

(m/sec) 

C2 

(m/sec) 

PI  7 
(ltR/m3) 

P2, 

(kg/m^) 

(G-9a,  G-9b,  G-9c) 

5 

50 

100 

200 

2500 

1500 

2700 

(G-lOa,  G-lOb,  G-lOc) 

10 

50 

100 

200 

2500 

1500 

2700 

(G-lla,  G-llb,  G-llc) 

5 

10 

100 

200 

2500 

1500 

2700 

(G-12a,  G-12b,  G-12c) 

10 

10 

100 

200 

2500 

1500 

2700 

Tunnel 

Depth 

(H) 

“  100  m 

(G-13a, 

G-13b, 

G-13c) 

5 

50 

100 

200 

2500 

1500 

2700 

(G-14a, 

G-14b, 

G-14c) 

10 

50 

100 

200 

2500 

1500 

2700 

(G-15a, 

G-15b, 

G-15c) 

5 

10 

100 

200 

2500 

1500 

2700 

(G-16a, 

G-16b, 

G-16c) 

10 

10 

100 

200 

2500 

1500 

2700 

Xg  =■  -52 

I  m,  Ax 

=  2  m,  nd 

=  51. 

2a 

“2  m, 

and 

N  “  3. 

100 


TABLE  IV-4 


TWO-LAYER  MODEL  SHOWING  HIGHER  ORDER  REFLECTIONS  IN  THE  SURFACE  LAYER 

Ys  *  -(h+1)  m 
y  *  0 

H  “  50  o 


Detectors  in  the  Bedrock 


h 


Figure  Numbers 

N 

(a) 

Qi 

(G-17a,  G-17b, 

G-17c) 

0 

5 

50 

(G-18a,  G-18b, 

G-18c) 

1 

5 

50 

(G-19a,  G-19b, 

G-19c) 

2 

5 

50 

(G-20a,  G-20b, 

G-20c) 

3 

5 

50 

Detectors 

(G-21a, 

G-21b, 

G-21c) 

0 

5 

50 

(G-22a, 

G-22b, 

G-22c) 

1 

5 

50 

(G-23a, 

G-23b, 

G-23c) 

2 

5 

50 

(G-24a, 

G-24b, 

G-24c) 

3 

5 

50 

Q? 

Cl 

(m/sec) 

C2 

(m/sec) 

PI  1 

(kg/m^) 

P2- 

(kg/m^ 

100 

500 

2500 

1500 

2700 

100 

500 

2500 

1500 

2700 

100 

500 

2500 

1500 

2700 

100 

500 

2500 

1500 

2700 

at  Surface 


100 

500 

2500 

1500 

2700 

100 

500 

2500 

1500 

2700 

100 

500 

2500 

1500 

2700 

100 

500 

2500 

1500 

2700 

Xg  =  -52  m,  Ax  *  2  m,  nd  =  51,  and  2a  =  2  m. 
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motion  seismograms  plotted  to  the  largest  practical  amplitude  (effective  scale 
factor  of  200)  in  an  effort  to  display  the  combined  direct  and  scattered  compo¬ 
nents  of  the  detected  signals.  The  illustrations  identified  by  (b)  show  only 
the  direct  arrival  SH-wave  particle  motion  seismograms  plotted  to  a  scale  which 
is  one-half  that  of  the  ilustrations  identified  by  (a).  The  illustrations 
identified  by  (c)  show  only  the  scattered  SH-wave  particle  motion  seismograms 
plotted  to  a  scale  which  is  ten  times  that  of  the  illustrations  Identified  by 
(a). 


Discussions  are  presented  below  based  upon  selected  examples  of  the 
calculated  model  cases  presented  in  Appendix  G  to  evaluate  the  effects  of 
attenuation,  thickness,  and  SH-wave  velocity  in  the  surface  layer  and  the 
location  of  the  detectors  on  the  tunnel  detection  reponse. 

1 .  Effect  of  Surface  Layer  Attenuation 

Examples  of  synthetic  seismograms  from  Appendix  G  are  discussed 
below  to  demonstrate  the  effects  of  surface  layer  attenuation  on  the  detection 
of  the  tunnel  cavity.  For  this  purpose,  the  SH-wave  line  source  is  in  the  bed¬ 
rock  layer  and  the  tunnel  cavity  depth  is  selected  to  be  H  =  50  m  and  the  sur¬ 
face  layer  thickness  is  selected  to  be  5  m.  The  layer  attenuation  parameters 
to  be  compared  correspond  to  quality  factors  of  Qj  =  50  (0.5  dB/wavel ength )  and 
Ql  =  10  (2.5  dB/wavel ength ) ,  and  the  detectors  are  located  in  the  top  of  the 
bedrock  layer  (Figures  G-1  and  G-3  in  Appendix  G)  and  at  the  top  of  the  surface 
layer  (Figures  G-9  and  G-ll  in  Appendix  G). 

Figures  IV-16  and  lV-17  show  the  comparative  effects  of  surface 
layer  attenuation  on  the  tunnel -scattered  SH  waves.  In  Figure  IV-16,  where  the 
SH-wave  source  and  detectors  are  in  the  competent  bedrock  layer,  the  tunnel 
detection  responses  are  not  significantly  affected  by  the  attenuation  quality 
of  the  surface  layer.  In  comparison,  the  responses  in  Figure  lV-17,  where  the 
source  is  in  the  bedrock  layer  but  the  detectors  are  at  the  top  of  the  surface 
layer,  show  a  very  substantial  attenuation  of  the  scattered  SH-waves  from  the 
tunnel  when  the  quality  factor  of  the  surface  layer  is  changed  from  50  to  10. 
The  detected  waveforms  observed  at  the  top  of  the  surface  layer  also  lack  the 
finer  details  contained  in  the  waveforms  detected  below  the  surface  layer, 
especially  in  the  more  realistic  Qi  =•  10  case.  Figures  G-5  and  G-7,  G-13  and 
G-15  in  Appendix  G  may  be  used  to  compare  corresponding  surface  layer  attenu¬ 
ation  effects  on  SH  waves  scattered  from  a  tunnel  target  at  100  m  below 
surface. 


2.  Effect  of  Surface  Layer  Thickness 

Examples  of  synthetic  seismograms  from  Appendix  G  are  discussed 
below  to  demonstrate  the  effects  of  surface  layer  thickness  on  the  detection 
of  the  tunnel  cavity.  For  this  purpose,  the  SH-wave  line  source  is  in  the 
bedrock  layer  and  the  tunnel  cavity  depth  is  selected  to  be  H  =  50  m  and  the 
surface  layer  quality  factor  is  selected  to  be  =  50.  The  surface  layer 

thicknesses  to  be  compared  are  h  *  5  m  and  h  =  10  m,  and  the  detectors  are 

located  in  the  top  of  the  bedrock  layer  (Figures  G-1  and  G-2  in  Appendix  G) 

and  at  the  top  of  the  surface  layer  (Figures  G-9  and  G-10  in  Appendix  G). 
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FIGURE  IV-16.  SYNTHETIC  SEISMOGRAMS  OF  SCATTERED  SH-WAVES  ONLY  FROM  A  2-m  DIAMETER  CYLINDRICAL  CAVITY 
LOCATED  IN  A  TWO-LAYER  LOSSY  HALF-SPACE  AT  A  DEPTH  OF  50  m.  THE  SURFACE  LAYER  THICKNESS  IS  5  m.  THE 
CONSTITUTIVE  PARAMETERS  OF  THE  SURFACE  LAYER  AND  THE  BEDROCK  ARE  ci  -  200  m/sec,  pi  =  1500  kg /m  , 

C2  -  2500  m/sec,  AND  p2  =  2700  kg/m^.  THE  QUALITY  FACTOR  OF  THE  BEDROCK  IS  Q2  =  100.  THE 
SOURCE  AND  ARRAY  OF  DETECTORS  HAVING  A  DETECTOR  SPACING  OF  2  m  ARE  IN  THE  BEDROCK. 
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Figures  IV-18  and  IV-19  show  the  comparative  effects  of  surface 
layer  thickness  on  the  tunnel -scattered  SH  waves.  In  Figure  IV-18,  where  the 
SH-wave  source  and  detectors  are  in  the  competent  bedrock  layer,  the  tunnel 
detector  responses  are  not  significantly  affected  by  the  surface  layer  thick¬ 
ness.  In  comparison,  the  responses  in  Figure  IV-19,  where  the  source  is  in 
the  bedrock  but  the  detectors  are  at  the  top  of  the  surface  layer,  show  a  sig¬ 
nificant  (53  percent)  increase  in  propagation  time  of  the  scattered  waves  as 
well  as  a  noticeable  loss  of  the  finer  details  of  the  detected  waveform,  all 
of  which  is  attributable  to  the  added  surface  layer  thickness.  The  more  real¬ 
istic  case  where  the  quality  factor  of  the  surface  layer  is  Qi  =  10  instead 
of  Qj  =  50  may  be  examined  by  comparing  Figures  G-3  and  G-4,  G-11  and  G-12 
in  Appendix  G  to  show  even  greater  effects  of  layer  thickness  on  the  tunnel 
detection  response. 

3.  Effect  of  Surface  Layer  SH-Wave  Velocity 

Examples  of  synthetic  seismograms  from  Appendix  G  are  discussed 
below  to  demonstrate  the  effects  of  different  SH-wave  velocity  conditions  in 
the  surface  layer  on  the  detection  of  the  tunnel  cavity.  For  this  purpose, 
the  SH-wave  line  source  is  in  the  bedrock  layer  and  the  tunnel  cavity  depth 
is  selected  to  be  H  =  50  m  and  the  surface  layer  quality  factor  is  selected 
to  be  =*  50.  The  thickness  of  the  surface  layer  is  h  =  5  m,  and  the  SH-wave 
velocities  to  be  compared  are  C[  =  200  m/s  and  c^  =  500  m/s,  respectively,  for 
the  SH-wave  detectors  located  in  the  top  of  the  bedrock  layer  (Figures  G-1  and 
G-20  of  Appendix  G)  and  at  the  top  of  the  surface  layer  (Figures  G-9  and  G-24 
of  Appendix  G). 


Figures  IV-20  and  IV-21  show  the  comparative  effects  of  SH-wave 
velocity  in  the  surface  layer  on  the  tunnel -scattered  SH  waves.  In  Figure  IV- 
20,  where  source  and  detectors  are  in  the  competent  bedrock  layer,  the  tunnel 
detector  responses  are  not  significantly  affected  by  differences  in  the  surface 
layer  velocity.  (The  multiple  reflections  within  the  surface  layer  have  been 
numerically  suppressed  in  Figure  IV-20(a)  but  not  in  Figure  IV-20(b)].  In  com¬ 
parison,  the  responses  in  Figure  IV-21,  where  the  source  is  in  the  bedrock 
layer  but  the  detectors  are  at  the  top  of  the  surface  layer,  show  a  significant 
(21  percent)  decrease  in  propagation  time  of  the  scattered  waves  and  a  slight 
but  noticeable  increase  in  the  curvature  of  the  hyperbolic  pattern  characteris¬ 
tic  of  the  tunnel  target  as  a  localized  reflector.  In  this  case,  since  only 
the  SH-wave  velocity  of  the  surface  layer  was  changed  without  a  change  in  the 
quality  factor,  the  details  in  the  detected  waveforms  shown  in  Figures  IV-20 
and  lV-2i  are  essentially  unaffected. 

4.  Direct-Wave  Interference 


The  direct  wave  transmitted  from  the  SH-wave  source  into  the 
surface  layer  will  cause  SH  waves  to  be  trapped  in  the  low-velocity  layer, 
giving  rise  to  SH-wave  groundroll  interference.  The  trapped  SH  waves,  known 
as  Love  waves,  propagate  as  normal  modes  which  are  dependent  upon  the  layer 
thickness  and  may  be  excited  when  the  source  is  located  either  in  the  layer  or 
in  the  bedrock.  However,  because  of  the  elastic  parameter  contrasts  at  the 
layer  interface,  the  magnitude  of  Love  waves  excited  by  a  source  below  the 
low-velocity  surface  layer  is  significantly  less  than  that  produced  by  a  source 
located  in  the  layer.  Examples  of  Love-wave  groundroll  interference  are  illus¬ 
trated  in  Figures  G-18  through  G-20  and  Figures  G-22  through  G-24  of  Appendix  G. 
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SH-wnVE  SYNTHETIC  SEISMOGRAM 
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LOCATED  IN  A  TWO-LAYER  LOSSY  HALF-SPACE  AT  A  DEPTH  OF  50  m.  THE  SURFACE  LAYER  THICKNESS  IS  5  m.  THE 
CONSTITUTIVE  PARAMETERS  OF  THE  SURFACE  LAYER  AND  THE  BEDROCK  ARE  pj  *  1,500  kg/m^,  =  2,500  m/sec, 
AND  P2  -  2,700  kg/m^.  THE  QUALITY  FACTORS  OF  THE  SURFACE  LAYER  AND  BEDROCK  ARE  Qi  =  50  AND 
Q2  -  100,  RESPECTIVELY.  THE  SOURCE  IS  IN  THE  BEDROCK,  AND  THE  ARRAY  OF  DETECTORS 
HAVING  A  DETECTOR  SPACING  OF  2  m  IS  AT  THE  SURFACE  OF  THE  EARTH. 


V.- 


SPECIALIZED  TECHNIQUES  FOR  SH-WAVE  TUNNEL  DETECTION 


A.  Coherent  Interference  Reduction 


A  striking  result,  although  not  unexpected,  in  the  various  SH-wave 
synthetic  seismograms  presented  in  this  study,  including  those  computed  for 
the  lossy  two-layer  model,  is  the  relatively  weak  amplitude  of  the  seismic 
pulse  signal  reflected  from  a  deep  tunnel  target  in  comparison  with  the  direct 
arrival  signals  that  travel  near  the  surface  betwe  n  the  source  and  detectors. 
These  large  non-target-related  signals  are  coherent  with  the  source  pulse  and 
the  desired  tunnel  reflections  since  they  emanate  from  or  are  excited  by  the 
same  source  that  produces  the  tunnel  reflection  signals.  Furthermore,  if 
there  are  other  reflecting  anomalies  in  the  tunnel  host  medium,  these  reflec¬ 
tions  will  also  be  coherent  with  the  desired  tunnel  reflection  signals.  Such 
coherent  Interference  signals  are  undesirable  for  two  reasons:  (1)  the  strong 
direct  arrival  signals  require  that  the  field  data  recording  system  have  a 
very  wide  dynamic  range  in  order  to  accurately  respond  to  and  simultaneously 
record  the  very  strong  direct  arrivals  and  the  very  weak  tunnel  reflections 
without  overload  distortion  or  loss  of  the  desirable  weak  signal  content;  (2) 
any  strong  interference  signals  that  arrive  at  the  detectors  simultaneously 
with  the  tunnel  reflections  will  mask  the  desired  signals,  thereby  reducing 
the  tunnel  response  signal-to-noise  ratio. 

Most  seismic  digital  data  recording  systems  have  a  dynamic  range 
of  15  digital  bits  (90  dB,  or  1  part  in  31,600)  which,  if  carefully  managed 
through  specialized  system  design  and  field  data  acquisition  techniques,  is 
sufficient  to  allow  the  weak  tunnel  reflection  signals  to  be  accurately 
recorded.  To  »'‘complish  this  goal,  the  tunnel  reflection  signal-to-nolse 
ratio  mist  be  maximized  using  every  practical  approach  and  resource  available 
since  strong  local  signals  at  the  detectors  and  anomalous  reflections  that 
can  mask  the  tunnel  reflections  will  generally  always  be  present.  The  tunnel 
target  response  signal-to-noise  ratio  can  be  maximized  in  two  ways;  namely, 
by  reducing  the  Interfering  and/or  dynamic  range  limiting  noise  and  by  enhanc¬ 
ing  the  tunnel  reflection  signal  by  discriminating  against  the  masking  noise. 
Techniques  for  applying  each  of  these  methods  are  discussed  below. 

1 .  Reduction  of  Direct  Arrival  Signals 

Direct  arrival  seismic  signals,  broadly  referred  to  as  ground 
roll,  travel  at  or  near  the  surface  either  as  surface  waves,  near-surface  body 
waves,  or  subsurface  interface  waves.  The  direct  propagating  path  of  the 
strongest  of  these  signals  may  be  considered  to  be  essentially  horizontal 
(parallel  to  the  ground  surface)  where  as  the  reflection  signals  from  deeper 
anomalies  occur  along  essentially  vertically  oriented  paths. 

The  strongest  ground  roll  interference  is  in  the  form  of  sur¬ 
face  (Rayleigh)  waves  which  are  two-dimensional  waves  locked  to  the  surface 
Interface  and,  in  a  practical  sense,  extend  only  a  few  wavelengths  below  the 
surface.  They  exist  and  propagate  at  approximately  the  same  speed  as  shear 
waves  by  Interchanging  wave  motions  and  energy  between  compresslonal  and  ver¬ 
tically  polarized  shear  displacement  motions  whenever  either  of  these  wave 
types  are  initiated  by  the  source  or  produced  by  wave-type  conversion  of 
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body-wave  source  signals  reflected  at  the  surface.  Thus,  the  SH-wave  tech¬ 
nique  evaluated  In  this  study  has  the  Inherent  advantage  of  not  directly  pro¬ 
ducing  surface  waves,  provided  that  Its  source  transducer  does  not  generate 
either  compresslonal  waves  or  vertically  polarized  shear  waves.  Two  other 
advantages  of  the  SH-wave  technique  which  discriminate  against  surface  wave 
Interference  and  that  have  been  emphasized  for  the  purpose  of  detecting  the 
tunnel  targets  of  Interest  are:  (1)  the  seismic  wave  frequencies  that  reflect 
efficiently  from  a  2-meter  diameter  tunnel  are  about  400  Hz  or  higher  and 
these  relatively  high  frequency  seismic  waves  are  strongly  attenuated  In  near¬ 
surface  geological  materials;  and  (2)  the  SH-wave  source  and  detectors  perform 
effectively  only  when  they  are  coupled  to  the  more  competent  materials  below 
the  surface  (typically  several  wavelengths  at  the  Rayleigh  wave  velocities 
characteristic  of  the  weathered  ground  surface  for  frequencies  of  400  Hz  and 
higher).  Therefore,  the  preferred  SH-wave  techniques  for  tunnel  detection 
Intrinsically  discriminate  In  several  ways  against  surface  wave  interference. 

However,  even  a  pure  SH-wave  source  will  generate  direct 
arrival  Interference  in  the  form  of  SH  body  waves.  If  the  source  Is  a  broad¬ 
band  Impulse  transducer,  as  It  must  be  to  fulfill  the  performance  requirements 
developed  In  this  study,  it  will  also  produce  unnecessary  spectral  components 
in  the  source  pulse  signal  which  are  below  400  Hz.  Since  these  low  frequency 
source  signal  components  are  not  attenuated  as  strongly  as  the  higher  fre¬ 
quency  components,  they  will  cause  unwanted  strong  interference  at  the  detec¬ 
tors.  However,  since  this  Interference  Is  not  In  the  useful  frequency  range 
for  detecting  the  small  tunnels  of  interest,  they  may  be  removed  by  means  of 
high-pass  filters  in  the  analog  preamplifier  circuits  following  each  seismic 
detector  transducer.  Filters  having  sharp  low-cut  response  and  high  stop-band 
attenuation  are  needed  to  effectively  remove  this  source  of  low-frequency 
Interference  which  may  occur  both  as  near-surface  direct  arrival  signals  and 
as  deeper  reflections.  Therefore,  an  effective  SH-wave  tunnel  detection  sys¬ 
tem  must  be  designed  to  effectively  reject  seismic  signals  at  frequencies 
below  about  400  Hz. 

Another  method  for  reducing  the  requirement  for  recording  the 
strong  direct  arrival  signals  is  simply  to  blank  the  input  to  the  recording 
system  for  an  accurately  known  time  delay  at  the  beginning  of  each  recorded 
detector  signal  trace.  Since  the  demand  for  wide  dynamic  range  is  imposed  by 
the  signals  that  arrive  during  the  earliest  trace  time  following  the  source 
pulse,  blanking  or  strongly  attenuating  this  part  of  the  trace  will  eliminate 
the  strong  interference.  This  blanking  process  must  generally  be  designed  as 
part  of  the  first  stage  of  the  analog  preamplifier  in  each  detector  channel  in 
order  to  prevent  overloading  the  amplifier  circuits  and  it  must  operate  with¬ 
out  Introducing  any  other  forms  of  distortion  or  electrical  transient  distur¬ 
bance  In  the  signals  to  be  recorded  after  the  end  of  the  blanking  Interval. 

In  this  regard,  most  digital  seismic  recording  systems  are  capable  of  starting 
their  digital  recording  after  a  preset  time  delay;  however,  they  are  not 
necessarily  equipped  to  block  input  circuit  overloading  if  it  should  occur. 
Therefore,  detector  circuit  blanking  in  a  hlgh-resolutlon  seismic  recording 
system  Is  a  technique  that  may  require  special  electronic  circuits  to  be 
designed  and  added  to  presently  available  equipment.  Furthermore,  as  will  be 
discussed  later  in  connection  with  seismic  static  corrections,  there  is  an 
auxiliary  need  to  record  direct  arrival  signals  as  an  adjunct  to  successful 
tunnel  detection,  thus  complicating  the  use  of  the  trace-blanking  technique. 
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2. 


Tunnel  Reflection  Signal  Enhancement 


The  amplitude  of  a  seismic  pulse  signal  reflected  from  the 
tunnel  target  is  directly  proportional  to  the  source  signal  amplitude  (or, 
equivalently,  to  the  source  energy)  which  causes  the  reflection.  Thus,  tunnel 
reflections  may  be  increased  by  increasing  the  source  level  provided  that  no 
undesirable  effects  such  as  loss  of  high-frequency  spectral  content  through 
plastic  deformation  of  the  ground  at  the  point  of  source  coupling  or  by 
requiring  unwieldy  source  devices  or  power  supplies.  However,  a  low-energy 
source  may  be  employed  for  the  same  purpose  by  repetitively  operating  it  to 
provide  signal  averaging  enhancement  of  the  weak,  tunnel  reflections.  This 
process,  more  precisely  known  as  repetitive  transient  signal  averaging,  is 
capable  of  enhancing  only  the  signal-to-noncoherent-noise  ratio  since  any 
interference  signals  that  are  coherent  with  the  desired  tunnel  reflections 
will  be  equally  enhanced.  Therefore,  this  process  has  the  effective  advan¬ 
tage  of  combining  the  energy  of  each  repeated  source  pulse  (by  Caking  Che 
repetitive  source  pulse  time  necessary  to  do  so)  to  finally  correspond  to  the 
energy  level  of  a  larger  source.  However,  to  be  successful  in  significantly 
increasing  Che  effective  source  level,  the  pulse  signal  produced  by  the 
repetitive  source  transducer  must  have  an  output  pulse  (radiated  into  the 
geologic  medium)  which  is  particularly  stable  in  pulse  timing  and  waveform. 

A  new  SH-wave  source  transducer  that  operates  on  the  arc  discharge  principle 
is  available  to  meet  these  requirements.  The  details  of  this  device  design 
are  presented  later. 

Spatial  signal  averaging  (or  spatial  filtering)  is  a  means  for 
reducing  the  masking  effects  of  seismic  reflections  that  are  simultaneously 
received  at  the  detectors.  This  method  Involves  adjusting  the  times  of  two  or 
more  detector  output  traces  recorded  from  a  given  source  pulse  but  at  differ¬ 
ent  detector  positions  so  as  to  cause  the  summation  of  the  traces  to  overlay 
and  advantageously  add  the  reflected  signals  from  the  tunnel  target  while,  as 
a  result  of  this  special  trace  time  adjustment,  the  signals  which  come  from 
other  reflectors  located  away  from  the  tunnel  combine  in  a  more  random  way. 
Thus,  by  effectively  focusing  the  detector  array  on  the  reflections  coming 
from  a  given  point  (the  tunnel  target)  in  the  subsurface  medium,  those  reflec¬ 
tions  will  be  enhanced  relative  to  the  other  coherent  signals  whose  origins  or 
reflections  are  located  away  from  the  tunnel  position.  To  accomplish  this 
target  signal  enhancement  data  processing  analysis,  however,  the  detector 
positions  must  be  relatively  close  together  to  eliminate  spatial  aliasing  of 
the  focal  point  mentioned  above  and  the  information  used  to  adjust  the  trace 
times  must  be  very  accurate  since  the  pulse  waveforms  in  the  tunnel  reflec¬ 
tions  must  be  overlaid  with  good  precision  in  order  to  produce  the  desired 
signal  enhancement  effect.  The  Information  used  in  the  trace  time  adjustments 
includes  data  on  the  spacings  and  elevations  of  the  source  and  detector  posi¬ 
tions  and  estimates  or  measurements  of  the  sesimic  wave  velocity  profiles 
along  the  presumed  propagation  paths  for  the  tunnel  reflections,  starting  at 
the  source  location  and  ending  at  the  detector  locations  being  analyzed. 

Since  these  time  adjustments  oust  be  more  accurate  than  the  time  period  of  the 
basic  frequency  characterizing  the  pulse  waveform,  the  realization  of  this 
demanding  form  of  spatial  signal  averaging  must  be  viewed  as  an  investigative 
development  to  be  studied  and  perfected  as  part  of  a  specialized  data  process¬ 
ing  and  analysis  program  for  hlgh-resolutlon  tunnel  detection. 
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Static  Corrections 


Static  corrections  of  seismic  records  pertain  to  time  adjust¬ 
ments  of  the  digitally  recorded  data  to  facilitate  the  spatial  signal  process¬ 
ing  techniques  such  as  the  spatial  filtering  technique  described  above  and  to 
adjust  the  Individual  detector  traces  for  variations  In  source  and  detector 
spaclngs  and  elevations  as  well  as  variations  In  the  seismic  wave  velocity 
such  as  may  be  Introduced  by  the  irregularly  weathered  surface. 

In  the  Intended  military  field  application  environment  of  the 
tunnel  detection  system  In  Korea,  the  terrain  relief  Is  generally  extreme  and 
will  require  careful  static  parameter  measurements  to  yield  high-quality 
tunnel  search  and  detection  results.  Accurate  seismic  propagation  velocity 
profiles  along  the  survey  traverse  In  high  relief  terrain  are  difficult  to 
Infer  from  simple  field  tests  and  may  be  assumed  to  vary  significantly  between 
contrasting  land  forms  such  as  ridges  and  valleys.  To  provide  this  needed 
velocity  data,  specialized  field  measurements  may  be  performed  as  part  of  each 
tunnel  detection  sounding  to  be  separately  analyzed  to  yield  site  specific 
velocity  Information  for  each  new  source  position  and  detector  spread  occupied 
along  the  survey  traverse.  These  measurements  entail  recording  the  direct 
arrival  signals  without  specific  concern  for  acquiring  reflections  from  pos¬ 
sible  tunnel  targets.  For  example,  when  the  geologic  section  being  surveyed 
consists  of  a  surface  soil  layer  overlying  bedrock,  the  tunnel  detection  sys¬ 
tem  will  be  suitable  for  conducting  reflection  mapping  surveys  of  the  depth 
to  the  bedrock  Interface  under  and  along  the  detector  spread.  An  effective 
method  for  mapping  the  depth  profile  of  Irregular  bedrock  Interfaces  has  been 
demonstrated  by  Hunger,  et  al.  (1984)  and  applied  to  low-velocity  soil  layers 
20-30  m  thick  overlying  karst  limestone  bedrock  by  Singh  (1984a;  1984b).  In 
this  method,  designated  as  the  reflection  window  technique,  the  source  offset 
distance  relative  to  the  near-trace  detector  position  must  be  sufficiently 
large  to  allow  the  bedrock  reflections  (traveling  through  the  surface  layer) 
to  be  detected  before  the  direct-ground  roll  interference  arrivals  but  not  so 
large  as  to  allow  the  faster-traveling  refraction  waves  (along  the  surface 
layer-bedrock  Interface)  to  mask  the  bedrock  reflections.  This  same  survey 
layout  can  also  be  used  to  conduct  refraction  survey  measurements  along  the 
soil-bedrock  Interface  to  yield  Information  on  the  velocity  in  the  bedrock 
layer.  Both  compresslonal  and  shear  waves  may  be  employed  In  these  measure¬ 
ments.  The  detector  spaclngs  considered  optimum  for  tunnel  detection  will 
also  be  appropriate  for  both  reflection  window  mapping  and  for  refraction 
surveys. 

In  addition  to  utilizing  direct  arrival  measurements  for  deter¬ 
mining  accurate  velocity  determinations,  supplemental  measurements  and  quick- 
look  field  analyses  can  also  be  performed  on  the  direct  arrival  signals  to 
determine  whether  or  not  the  source  and  detector  transducers  are  coupled  In 
ground  of  sufficiently  high  quality  to  support  SH-wave  operations  in  the 
desired  frequency  range  of  400-1,600  Hz  for  efficient  detection  of  the  tunnel 
targets  of  Interest. 

B.  Use  of  Tunnel  Reflection  Models  In  Data  Analysis 

The  analytical  and  numerical  models  for  SH-wave  reflections  from 
cylindrical  cavities  developed  In  this  study  have  provided  useful  Insights 
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into  the  nature  of  the  tunnel  detection  problem  as  well  as  guidance  for  devel¬ 
oping  an  effective  tunnel  search  and  detection  system.  In  addition,  these 
models  and  their  variations  and  extensions  also  offer  a  powerful  aid  In  ana¬ 
lyzing  the  expected  field  survey  data.  The  value  of  these  models  Is  that, 
given  the  experimental  parameters  of  the  geologic  media  being  explored,  the 
models  can  provide  approximate  predictions  of  any  tunnel  reflection  responses 
that  may  be  present.  Even  though  such  Idealized  predictions  only  approximate 
the  actual  field  data,  they  are  free  of  geologic  noise  and  can  provide  a  use¬ 
ful  template  of  seismic  traces  and  seismogram  patterns  that  will  be  represen¬ 
tative  of  the  tunnel  response.  Then,  by  using  such  trial  templates  containing 
expected  target  responses,  the  model  calculations  may  be  employed  as  matched 
filters  to  be  correlated  and  compared  with  the  field  survey  data  to  extract 
any  similar  responses  that  may  be  present  In  the  experimental  records. 

By  applying  this  forward  modeling  approach  to  the  field  data,  the 
Information  content  In  an  entire  seismic  section  of  data  becomes  Involved  In  a 
unique  process  of  tunnel  target  resoonse  matching  and  extraction  without  the 
problems  of  Instabllty  or  nonconveigence  often  encountered  in  geophysical  data 
Inversion  when  applied  to  low  signal-to-noise  ratio  field  recordings. 

In  order  to  apply  this  model-based  data  analysis  technique,  the 
field  data  must  be  processed  to  a  certain  reference  status  and  format  against 
which  a  similar  ensemble  of  model-predicted  data  for  various  possible  tunnel 
locations  with  respect  to  the  experimental  survey  traverse  may  be  compared. 

The  processing  and  preconditioning  the  field  data  consists  mainly  of  applying 
static  corrections  and  velocity  analysis  for  the  field  site  conditions  associ¬ 
ated  with  the  data  and  application  of  any  data  preprocessing  options  such  as 
the  spatial  signal  averaging  technique  designed  to  Improve  the  tunnel  response 
signal-to-nolse  ratio.  By  means  of  more  advanced  models  developed  to  Include 
compresslonal- waves  and  vertically  polarized  shear  (SV)  waves,  more  complete 
Insights  into  the  complete  tunnel  detection  problem  will  be  gained.  With  this 
Information  available,  additional  advanced  data  processing  techniques  may  be 
considered  which  employ  both  radial  and  transverse  horizontal  particle  motion 
components  recorded  along  the  seismic  survey  line  and  a  more  general  analysis 
based  upon  three-component  vector  displacements  recorded  at  each  detector 
position.  An  important  method  used  In  the  processing  and  analysis  of  oriented 
three-component  seismic  survey  data  is  that  known  as  the  "polarization  method" 
reported  by  Gal'perin  (1985)  in  connection  with  vertical  seismic  profile  data 
analysis.  With  this  method,  reflected  signals  from  the  tunnel  target  can  be 
separated  from  certain  noise  components  by  deriving  spatial  displacement  vec¬ 
tors  versus  time  (depth)  so  as  to  maximize  the  desired  tunnel  reflection 
signal-to-noise  ratio.  This  method  should  be  particularly  applicable  to  hlgh- 
resolutlon  seismic  surveys  In  complex  geological  formations  as  well  as  provide 
an  effective  means  of  minimizing  Che  Interference  of  surface  waves  and  other 
near-surface  coherent  seismic  noise. 

C.  Field  Data  Acquisition  System 

The  seismic  data  acquisition  system  required  to  provide  effective 
tunnel  search  and  detection  Is  a  hlgh-resoluclon  reflection  survey  system 
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equipped  specifically  to  employ  the  SH-wave  techniques  identified  in  the  fore¬ 
going  model  studies.  In  particular,  this  system  must  consist  of: 

(1)  Portable  Seismic  Data  Recording  System 

*  48-channel  digital  data  acquisition  controller 

*  48-channel  preamplifier  and  filter  unit 

*  Digital  tape  recorder 

*  Quick-look  post-acquisition  field  processor 

*  Data  display  plotter 

*  Power  supply 

(2)  Seismic  Detectors  and  Field  Cable 

*  Three-component  geophone  assembly  designed  for  borehole 
clamp  coupling  and  orientation  (48  units) 

*  Multiconductor  seismic  cable  designed  for  use  with  three- 
component  geophones  with  recorder  system  connection  jumpers 
(typical  detector  spread  length:  240  m) 

*  Roll-along  switch  capable  of  being  adapted  for  use  with 
one-,  two-,  or  three-component  detector  stations  (typically 
240  geophone  pairs,  programmable,  manual  operation) 

(3)  SH-Wave  Seismic  Source 

*  Arc  discharge  pulse  transducer  designed  for  operation  in 
water-filled  boreholes  (equipped  with  orientation  capability 
and  bowspring  centralizer) 

*  Surface  power  supply  and  control  unit 

(4)  Borehole  Auger  Drill 

*  Truck-mount  drill  system  (4-inch  diameter,  50-foot  depth 
capacity) 

1 .  Portable  Seismic  Data  Recording  System 

The  high-resolution  seismic  data  acquisition  system  must  be 
capable  of  accurately  recording  multichannel  signals  in  the  frequency  range 
of  400  Hz  to  1600  Hz  with  trace  time  durations  up  to  about  0.3  second  in  each 
channel.  This  frequency  response  will  require  digital  sampling  of  each  of  the 
preamplifier  analog  detector  signals  at  a  rate  of  at  least  4,000  samples  per 
second  and  preferably  higher.  The  dyanmlc  range  of  the  system  must  be  as  wide 
as  practical  but  typically  equivalent  to  15-bit  digital  waveform  amplitude 
resolution  plus  sign  bit,  corresponding  a  dynamic  range  of  90  dB  in  the  digi¬ 
tized  sample  word.  Typical  analog-to-dlgltal  converters  used  for  this  purpose 
provide  a  12-bit  floating  point  mantissa  plus  sign  bit  and  this  72-dB  dynamic 
range  Is  supplemented  with  a  three-or  four-bit  gain  range  control  of  the  ana¬ 
log  signal  amplifier  ahead  of  the  converter  to  provide  an  additional  18  dB  or 
24  dB  overall  system  dynamic  range.  Thus,  for  a  48-channel  recording  system 
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operating  with  such  an  effective  16-bit  word  length  acquired  at  4,000  sample 
words  per  second  and  for  record  lengths  of  0.5  second,  the  digital  data  storage 
requirements  are  1.54  million  bits  per  shot.  If  each  sample  word  is  recorded  as 
two  8-blt  bytes,  each  shot  record  will  require  96,000  bytes  of  digital  magnetic 
tape  or  memory  storage.  A  total  of  1,875  shots  having  96,000  bytes  each  can  be 
recorded  in  SEG-D  format  on  one  10-inch  6,250  bpl  reel  of  magnetic  tape. 

Certain  specialized  features  are  recommended  to  be  Incorporated 
in  the  data  acquisition  system  to  improve  its  high-frequency  response  and  to 
limit  its  response  to  the  most  effective  frequency  range  for  detecting  SH-wave 
reflections  from  tunnels.  The  definition  of  these  features  and  their  techni¬ 
cal  characteristics  are: 

(1)  Low  Cut  Filters:  400  Hz  cut-off  frequency 

-60  dB/oct  stop-band  attenuation 

(2)  Passband  Emphasis:  +24  dB/oct  400  Hz  to  1600  Hz 

(3)  Digital  Sampling  Rate:  Increase  from  4,000  samples/sec 

to  8,000  samples/sec  (goal) 

(4)  Alias  Filters:  2,000  Hz  cut-off  frequency 

-72  dB/oct  stop  band  attenuation 

(5)  Digital  Parallel  Data  Output  Interface. 

A  digital  recording  system  of  this  type  will  effectively  eliminate  low-frequency 
ambient  noise  and  will  compensate  for  much  of  the  frequency-dependent  absorption 
loss  imposed  by  the  geologic  medium.  Such  a  system  configured  to  have  48  chan¬ 
nels  will  accommodate  16  three-component  seismic  detectors. 

Quick-look  post-acquisition  data  processing  will  be  used  to  yield 
information  on  the  quality  of  the  recorded  data  and  the  nature  of  the  subsurface 
conditions  currently  being  encountered  along  the  survey  line.  The  processing 
capabilities  needed  include  digital  filtering,  trace  cross-correlation,  normal 
moveout  corrections,  common  offset  gathers,  and  common  depth  point  stacking.  On 
the  anticipation  that  the  quality  and  usefulness  of  the  high-resolution  seismic 
exploration  data  necessary  for  detecting  tunnels  will  be  difficult  to  determine 
in  practice,  this  quick-look  in-field  processor  can  ensure  that  further  more 
extensive  data  analysis  will  yield  the  desired  results. 

2.  Seismic  Source  Transducer 


High-resolution  seismic  exploration  for  tunnels  requires  that  the 
source  and  detector  transducers  be  coupled  to  competent  ground  so  as  to  convey 
elastic  wave  signals  having  frequency  content  up  to  about  2,000  Hz.  Such  cou¬ 
pling  cannot  be  achieved  with  either  the  source  or  the  detectors  emplaced  at 
the  ground  surface  unless  the  competent  bedrock  tunnel  host  formation  happens 
to  outcrop  along  the  desired  seismic  survey  line.  Therefore,  in  general,  the 
source  and  detector  transducers  oust  be  designed  to  operate  in  boreholes  drilled 
to  reach  subsurface  ground  conditions  which  are  competent  enough  to  support 
high-frequency  elastic  waves. 
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An  efficient  broadband  seismic  pulse  transducer  technique  was 
recently  developed  by  Southwest  Research  Institute  for  use  in  high-resolution 
borehole  seismic  applications  (Peters  and  Owen,  1985).  This  transducer  is 
capable  of  generating  accurately  repetitive  pressure  pulse  waveforms  in  fluid- 
filled  boreholes  to  produce  elastic  wave  pulses  in  the  drilled  geological 
formation.  The  transducer  operates  by  producing  a  high-energy-rate  electric 
discharge  in  a  self-contained  aqueous  electrolyte  solution  to  produce  intense 
pressure  pulses.  Electrical  energy  up  to  approximately  1,000  joules  is  stored 
on  a  high-voltage  capacitor  and  Is  discharged  in  the  arc  chamber  within  a  time 
interval  of  a  few  hundred  microseconds  to  produce  a  positive  pressure  pulse 
having  a  peak  pressure  of  several  atmospheres.  This  arc  discharge  generator 
may  be  operated  at  a  pulse  repetition  rate  of  approximately  10  pulses/min  in 
boreholes  4-6  Inches  in  diameter. 

Figure  V-1  Illustrates  the  general  design  concept  of  the  elec¬ 
tric  arc  discharge  pulse  transducer.  This  figure  shows  the  cylindrical  elec¬ 
trodeless  arc  discharge  chamber  in  which  a  small  channel  of  liquid  electrolyte 
is  vaporized  to  produce  the  plasma  arc  discharge  pressure  pulse.  This  pulse 
is  coupled  to  the  borehole  fluid  through  the  flexible  cylindrical  sleeve  form¬ 
ing  the  electrolyte  chamber.  Such  pressure  pulses  are  generated  each  time  the 
charged  high-voltage  capacitor  is  switched  across  the  electrodes  on  the  upper 
and  lower  ends  of  the  chamber. 


♦HV  (SWITCHED) 


FIGURE  V-1.  ELECTRODELESS  ARC  DISCHARGE  TRANSDUCER  FOR  GENERATING 
ASSYMETRICAL  LATERAL  FORCE  IN  A  BOREHOLE 
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The  high-density  reaction  mass  shown  surrounding  one  half  of 
the  arc  discharge  chamber  in  Figure  V-1  serves  to  produce  the  asymmetrical 
lateral  force  in  the  borehole  necessary  for  exciting  horizontally  polarized 
shear  waves.  The  asymmetrical  force  is  oriented  in  the  plane  A-A'  and  hori¬ 
zontally  polarized  shear  waves  are  radiated  away  from  this  point  source  trans¬ 
ducer  in  the  plane  B-B'. 

A  vertical  array  consisting  of  several  devices  of  this  type 
spaced  at  a  distance  about  0.5  meter  (less  than  one-half  wavelength  at  the 
highest  frequency  of  interest)  and  timed  to  produce  a  downward  end  fire- 
directed  pulse  can  provide  a  very  efficient  controlled  source  of  SH  waves  for 
tunnel  detection.  For  example,  an  array  of  ten  SH-type  arc  discharge  elements 
in  a  4.5-meter  long  array  will  produce  a  30-degree  half-power  radiation  beam 
pattern  at  a  frequency  of  1,000  Hz.  An  array  of  this  type  would  provide  a 
total  energy  of  5,000  joules  per  pulse  in  a  water-filled  borehole  approxi¬ 
mately  4  inches  in  diameter  drilled  to  the  depth  required  to  reach  the  compe¬ 
tent  geological  medium. 

In  the  general  form  described  above,  the  basic  borehole  arc 
discharge  transducer  can  be  designed  to  generate  either  SH-wave  pulses  (with 
reaction  mass  attached)  or  compress ional  wave  pulses  (without  reaction  mass 
attached).  Preliminary  SH-wave  pulse  experiments  were  recently  conducted 
using  the  U.S.  Bureau  of  Mines  500-joule  arc  discharge  source  probe  with  a 
lead  half-cylinder  reaction  mass  attached.  Figure  V-2  illustrates  the  results 
of  these  preliminary  tests  in  which  the  source  transducer  having  the  reaction 
mass  attached  was  placed  in  one  water-filled  borehole  in  a  limestone  formation 
and  a  hydrophone  detector  was  placed  in  a  second  water-filled  borehole  located 
approximately  3.05  meters  away  from  the  source.  Figure  V-2(a)  shows  the  com- 
presslonal  wave  pulse  waveform  transmitted  horizontally  between  the  boreholes 
without  the  lead  reaction  mass  attached  to  the  source  transducer.  With  the 
reaction  mass  attached,  the  source  probe  was  rotated  until  the  strongest  shear 
wave  obtainable  was  observed  superimposed  on  the  compress ional  wave  pulse. 
Figure  V-2(b)  shows  the  superimposed  SH  wave  and  indicates  the  shear  wave 
velocity  in  the  limestone  to  be  about  50.1  percent  of  the  compressional  wave 
velocity.  The  hydrophone  was  then  lowered  in  two-meter  depth  intervals  below 
the  source  to  show  the  Improved  SH-wave /P-wave  ratio  characteristic  of  the 
downward  (and  upward)  directed  SH-wave  radiation  from  the  source  transducer. 
The  SH  wave  pulse  signals  detected  at  downward  depth  offsets  of  2,  4,  and  6 
meters  are  shown  in  Figures  V-2(c)  through  V-2(d).  The  SH-wave  pulse  signals 
observed  in  these  tests  are  judged  to  be  much  stronger  than  the  compressional 
wave  pulses  since  the  hydrophone  detector  used  in  the  experiments  is  not  well 
suited  for  the  detection  of  shear  waves.  These  tests  indicate  that  the  asym¬ 
metrical  force  generated  around  the  circumference  of  the  borehole  by  this 
arrangement  is  effective  in  generating  SH  wave  pulses. 

The  SH-wave  source  transducer  must  be  oriented  so  that  its 
asymmetrical  force  is  directed  transverse  to  the  seismic  survey  line.  For 
this  purpose,  the  SH-wave  source  transducer  must  be  equipped  with  a  clamping 
mechanism  by  which  it  may  be  held  at  a  fixed  orientation.  A  detachable  metal 
suspension  tube  of  square  cross-section  may  be  used  to  handle  and  orient  the 
source  transducer  in  the  proper  direction  and  hold  it  in  place  in  the  rela¬ 
tively  shallow  boreholes  anticipated  for  field  use  while  the  clamping  mecha¬ 
nism  is  operated. 
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3. 


Seismic  Detector  Transducers 


Precision  digital-grade  three-component  seismic  detectors  are 
required  for  use  with  the  high-resolution  seismic  exploration  system.  In 
addition,  each  of  these  detectors  nust  be  equipped  with  a  clamping  mechanism 
by  which  it  may  be  ridigly  coupled  to  the  borehole  wall.  With  such  mechanical 
coupling,  the  detector  boreholes  are  not  required  to  have  a  borehole  fluid 
content . 


The  seismic  detectors  recommended  for  use  with  the  system  are 
acceleration  sensitive  geophones  mounted  in  a  three-component  configuration. 
The  seismic  detectors  should  be  mounted  on  a  specialized  triaxial-sensor 
borehole-conforming  contact  shoe  which  will  be  rigidly  locked  against  the 
borehole  wall  by  means  of  a  pneumatically-expanded  rubber  pressure  boot.  A 
detachable  metal  suspension  tube  of  square  cross-section  or  other  device  may 
be  attached  to  this  downhole  assembly  to  allow  it  to  be  handled  and  oriented 
in  the  proper  direction  and  hold  it  in  place  while  the  clamping  boot  is  being 
inflated.  An  air  compressor  and  pressure  regulating  control  unit  will  be 
needed  to  supply  pneumatic  inflation  of  the  boots.  The  combination  geophone 
cable  and  inflation  air  hose  on  the  downhole  detector  units  will  be  designed 
to  allow  retrieval  of  the  unit. 

A.  Custom  Seismic  Survey  Cable 

Reflection  and  refraction  seismic  exploraton  for  tunnel  targets 
will  require  a  special-purpose  detector  layout  cable  having  a  sufficient  num¬ 
ber  of  electrical  conductor-pair  circuits  and  interval-spaced  take-outs  to 
connect  the  planned  spread  of  seismic  detectors  to  the  48-channel  seismic  data 
acquisition  system.  This  cable  must  have  groups  of  three  conductor-pair  take¬ 
outs  per  detector  station  to  accommodate  the  use  of  three-component  geophone 
sensors  at  each  station  when  desired.  The  three  take-outs  at  each  station 
should  be  spaced  0.5  meter  apart  and  color  coded  to  represent  the  connections 
for  the  x-horizontal ,  y-horizontal ,  and  z-vertical  geophones  contained  in  the 
borehole  detector  assemblies.  Each  borehole  detector  station  must  be  emplaced 
along  the  seismic  line  at  approximately  2.5  meter  spacings  to  prevent  spatial 
aliasing  in  the  field  data,  and  the  cable  interval  spacing  between  take-outs 
should  be  2.75  meters.  The  cable  should  also  be  segmented  in  hanks  of  66 
meters  length  with  24  detector  stations  per  segment  and  color-coded  end  con¬ 
nectors  which  will  allow  each  segment  to  be  moved  forward  and  properly  con¬ 
nected  to  the  leading  end  of  the  cable  on  the  survey  line  after  each  24  source 
shots  are  completed.  Since  the  24  three-component  detectors  will  require  72 
conductor  pairs,  the  total  number  of  conductor  pairs  for  a  three-segment  cable 
layout  (3  X  66  m  *  198  m  length)  will  be  216  pairs.  This  number  of  pairs  can 
be  conveniently  switched  into  the  48-channel  data  recording  system  using  a 
standard  240-detector-station  roll-along  switch.  This  relatively  large  number 
of  conductor  pairs  (216  pairs)  can  be  reduced  to  144  pairs  and  still  allow  the 
use  of  three  cable  segments  if  multiple  take-outs  are  used  in  operating  the 
composite  cable.  With  this  arrangement,  the  multiple  take-outs  on  the  cond  ;c- 
tor  pairs  are  separated  48  detector  stations  apart  so  that  no  two  detectors 
are  simultaneously  attached  to  the  same  pair  to  Introduce  erroneous  signals 
into  the  recorded  data.  The  roll-along  switch  is  programmed  to  switch  the 
conductor  pairs  by  being  back-wired  to  pick  up  appropriate  repeat  connections 
to  the  previously  used  pairs  having  detectors  at  new  forward  take-out  posi¬ 
tions.  The  66-meter  cable  hanks  having  24  three-component  take-outs  and  144 
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conductor  pairs  are  estimated  to  have  a  weight  of  approximately  50  pounds 
using  02S  AWG  conductors  in  each  pair  for  good  cable  strength  and  reliability 
in  the  field. 


Next,  by  operating  the  seismic  data  acquisition  system  in  a 
split-spread  configuration  composed  of  two  three-hank  cables  as  described 
above,  the  same  cable  design  and  240-detector  roll-along  switch  can  be  used 
to  cover  twice  the  total  traverse  line  length  for  each  setup  of  the  data 
recording  system  (truck).  This  is  accomplished  by  means  of  a  patch  panel  by 
which  the  individual  geophones  (conductor  pairs)  are  prewired  to  configure  the 
desired  detectors  (i.e.,  single-vertical,  single-horizontal,  dual-horizontal, 
or  three-component)  to  be  sequentially  selected  by  the  240-detector  roll-along 
switch  without  overflow  of  the  number  switch  positions.  The  use  of  the  patch 
panel  also  allows  the  selected  detectors  to  be  directly  connected  to  their 
analog  preamplifiers  and  low-cut  filters  ahead  of  the  roll-along  switch  for 
maximum  signal-to-nolse  ratio. 

5.  Roll-Along  Switch 


The  system  must  be  equipped  with  a  240-circuit  roll-along  geo¬ 
phone  switch  to  accommodate  the  special  seismic  cable  and  detector  pairs  as 
configured  by  the  patch  panel  used  to  set  up  the  detector  stations  for  field 
operations.  Roll-along  switches  of  this  type  are  conventional  components 
which  are  commercially  available  from  several  suppliers. 

6.  Borehole  Auger  Drill 

A  light-duty  shallow  borehole  auger  drill  is  required  as  a  per¬ 
manent  part  of  the  high-resolution  seismic  system  since  efficient  operation  of 
the  system  is  expected  to  require  the  use  of  boreholes  for  source  and  detector 
coupling  on  a  regular  basis.  This  drill  should  have  a  3-1/2-inch  diameter 
borehole  drilling  capability  to  a  nominal  depth  of  15  meters  using  spade-bit 
and  auger  technology  (dry  drilling)  and  a  4-1/2-inch  diameter  helical  screw 
bit  designed  especially  for  enlarging  and  compacting  the  augered  hole  walls  to 
stand  open  without  the  need  for  casing  and  to  provide  compacted  hole  rigidity 
for  good  seismic  coupling.  The  boreholes  used  for  the  seismic  source  trans¬ 
ducer  will  require  water  to  be  present  or  added  for  proper  operation.  The 
boreholes  used  for  the  seismic  detectors  will  not  require  water  filling  for 
sensor  operation. 

7 .  Spare  Parts 

Spare  parts  are  required  for  the  various  system  components 
which  receive  routine  field  handling  such  as  the  source  and  detector  trans¬ 
ducers,  cables,  borehole  clamping  mechanisms,  borehole  drill,  and  related  sub¬ 
components.  Certain  spare  electronic  modules  are  also  required  to  maintain 
the  digital  data  recording  system  and  the  arc  discharge  source  transducer 
system. 


8.  Field  Vehicles 


Two  field  vehicles  are  required  to  support  the  operation  of  the 
high-resolution  seismic  exploration  system:  a  mobile  data  recording  and  sys¬ 
tem  transport  vehicle,  and  a  dedicated  drill  truck. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  theoretical  analyses  and  computer  evaluations  of  SH-wave  seismic 
reflection  techniques  carried  out  in  this  study  have  yielded  important  insights 
and  useful  results  which  demonstrate  the  feasibility  of  detecting  subversive 
intrusion  tunnels  of  military  interest.  Based  upon  these  results,  the  seismic 
detectability  of  tunnels  2  m  in  diameter  at  depths  of  100  m  or  more  in  hard 
rock  geological  formations  is  judged  to  be  technically  feasible  and  practical 
using  SH  waves,  provided  that  certain  equipment  capabilities,  field  application 
techniques,  and  data  analysis  methods  are  employed.  The  primary  findings  and 
conclusions  gained  from  this  study  are; 

(1)  Seismic  Detectability  of  Tunnels:  Horizontally  polarized  shear  (SH) 
waves  offer  certain  intrinsic  geophysical  advantages  for  detecting 
tunnels  based  upon  their  shorter  wavelength  at  practical  seismic 
frequencies,  their  reduced  tendency  to  excite  surface  ground  roll 
interference,  and  the  fact  that  SH  waves  do  not  suffer  wave-type 
conversion  losses  when  reflections  occur  at  geologic  interfaces 
oriented  parallel  to  the  shear  wave  polarization. 

Plane  wave  scattering  analysis  was  used  to  show  the  frequency  depen¬ 
dence  and  angular  dependence  of  SH-wave  scattering  from  cylindrical 
cavities.  The  effective  frequency  range  for  detecting  the  tunnels 
of  interest  is  the  two-octave  band  of  400-1,600  Hz  since  spectral 
components  in  the  incident  wave  below  400  Hz  provide  only  a  marginal 
response  from  the  small  tunnel  cavity  cross-section  and  spectral 
components  above  about  1,600  Hz  are  attenuated  significantly  and 
tend  to  be  scattered  from  the  many  natural  inhomogeneities  in  the 
rock  formation  that  are  smaller  than  the  tunnel. 

(2)  Coherent  Interference  and  Masking:  Interference  caused  by  direct 
seismic  waves  traveling  from  source  to  detector  was  illustrated  in 
synthetic  seismograms  derived  using  homogeneous  half-space  and  two- 
layer  half-space  models.  For  the  SH  waves  modeled,  the  amplitude 
ratios  of  reflections  from  a  lOO-meter  deep  tunnel  to  the  direct 
arrival  ranged  from  about  -32  dB  for  a  lossless  half-space  to  -55  dB 
for  the  case  where  the  seismic  detectors  are  located  at  the  surface 
of  a  low-quality  layer  overlying  the  tunnel  host  bedrock  medium. 

Table  VI-1  summarizes  these  reflection-to-direct  arrival  ratios  for 
the  various  half-space  models  and  tunnel  cavity  depths  evaluated. 

These  results  clearly  show  the  importance  of  providing  a  wide  dynamic 
range  in  the  field  data  acquisition  system.  The  results  also  present 
quantitative  information  on  the  influences  of  tunnel  target  depth  on 
the  SH-wave  reflections  as  well  as  the  effects  of  a  lossy  surface 
layer  on  the  detection  effectiveness. 

(3)  Effects  of  Low-Quality  Surface  Layer:  The  simplified  field  proce¬ 
dures  afforded  by  the  use  of  surface-coupled  seismic  detectors 
greatly  diminish  the  effectiveness  of  the  seismic  tunnel  detection 
concept.  In  the  model  cases  where  the  detector  responses  were 
located  at  the  top  of  a  thin  low-quality  (Qi  =  10)  surface  layer, 
the  layer  attenuation  introduces  reductions  in  the  tunnel  reflection 
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TABLE  VI- I 


RATIO  OF  SH-WAVE  TUNNEL  REFLECTIONS  TO  DIRECT 
SH-WAVE  INTERFERENCE  FOR  THE  LAYERED  HALF-SPACE  MODEL 


Model  Description 


^ SH-Wave  Tunnel  Reflection 
Direct  SH-Wave  Amplitude 

_ (dB) _ 

Tunnel  Depth  Tunnel  Depth 
=  50  m  =  100  m 


Two-layer  half-space  with  source 
and  detector  in  the  bedrock. 


Reference 
Figures  in 
Volume  II 
Appendix  G 


Surface  layer  thickness  h  =  5  m 
Bedrock  quality  factor  Q2  =  100 
Surface  Layer  Qi  =  50: 

Surface  Layer  Qj  ■  10: 


Surface  layer  thickness  h  ■  10  m 
Bedrock  quality  factor  Q2  “  100 
Surface  Layer  Qj  *  50: 

Surface  Layer  Qj  *  10: 


-26 

-26 


-26 

-26 


-32  G-1  &  G-5 

-32  G-3  &  G-7 


-32  G-2  &  G-6 

-32  G-4  &  G-8 


Two-layer  half-space  with  source 
with  source  in  the  bedrock  and 
detector  at  top  of  surface  layer. 


Surface  layer  thickness  h  =» 
Bedrock  quality  factor  Q2  * 
Surface  Layer  Qi  ■  50: 
Surface  Layer  Qi  »  10: 


Surface  layer  thickness  h  =• 
Bedrock  quality  factor  Q2  ** 
Surface  Layer  Qj  =  50: 
Surface  Layer  Qi  =  10: 


5  m 
100 


10 

100 


-26 

-44 


-34  G-9  &  G-1 3 

-49  G-11  &  G-15 


m 


-32 

-50 


-40  G-10  &  G-14 

-55  G-12  &  G-16 


Notes:  (1)  Cylindrical  SH-wave  source  (3)  Surface  layer  vsi  =  200  m/sec 

(2)  2-meter  diameter  tunnel  (4)  Bedrock  layer  vs2  “  2,500  m/sec 


amplitudes  of  approximately  24  dB  compared  with  the  reflection  ampli¬ 
tudes  detectable  in  the  top  of  the  bedrock.  The  practical  influences 
of  surface  geophone  coupling  efficiency  and  variability  among  conven¬ 
tionally  deployed  detectors  can  be  anticipated  to  introduce  addi¬ 
tional  losses  in  the  detection  of  tunnels.  The  combined  effects  of 
surface  layer  losses  imposed  on  the  tunnel  target  reflections  and 
the  stronger  direct  arrival  interference  when  the  detectors  are 
located  on  the  surface  lead  to  the  conclusion  that  effective  tunnel 
search  and  detection  cannot  be  achieved  using  surface-coupled  seis¬ 
mic  detectors. 

(4)  Specialized  Data  Acquisition  and  Processing  Techniques:  Field  data 
acquisition  and  analysis  of  SH-wave  reflections  from  tunnel  targets 
in  field  terrain  such  as  the  DMZ  of  Korea  require  special  consid¬ 
erations  and  techniques  not  normally  used  in  conventional  seismic 
exploration  practice.  These  requirements  stem  largely  from  the  high- 
frequency  seismic  signals  needed  for  resolving  the  small  tunnel  cavi¬ 
ties  of  interest.  They  include;  (1)  the  use  of  a  subsurface-coupled 
orientable  SH-wave  seismic  source  and  surbsurf ace-coupl ed  three- 
component  detectors  capable  of  being  oriented  relative  to  the  source 
orientation  and  survey  line  layout;  (2)  special  site-characterizing 
measurements  to  provide  detailed  information  on  the  behavior  of  the 
very  short  wavelength  signals  and  corrections  for  static  site  condi¬ 
tions  such  as  topographic  elevations,  seismic  velocity  variations  in 
near-surface  layers,  and  the  thickness  and  subsurface  shear-wave  cou¬ 
pling  quality  of  the  surface  materials;  and  (3)  the  use  of  optimum 
detector  layouts,  source  locations,  and  repetitive  signal  averaging 
to  enhance  the  tunnel  reflection  signals  during  the  field  data  acqui¬ 
sition  process  so  as  to  provide  sufficient  quality  data  for  subse¬ 
quent  processing  and  analysis. 

Givei  high-quality  field  data,  processing  techniques  for  analyzing 
these  data  must  be  applied  carefully  and  with  an  objective  view 
toward  the  results  because  of  the  complex  geological  structure  of 
the  tunnel  host  medium.  That  is,  no  substantial  experience  is  avail¬ 
able  to  indicate  whether  conventional  data  processing  algorithms  will 
be  effective  when  applied  to  the  high-resolution  seismic  survey  for 
tunnel  detection.  Therefore,  while  standard  methods  of  velocity 
analysis  and  spatial  filtering  are  planned  to  be  used  to  their  full¬ 
est  extent,  supplemental  concepts  designed  to  respond  specifically  to 
the  tunnel  target  seismic  responses  are  also  recommended  for  inves¬ 
tigation  and  potential  use.  For  this  purpose,  the  computer  models 
developed  in  this  study  can  be  used  to  provide  certain  predictions  of 
tunnel  target  responses  which,  through  proper  data  comparison  meth¬ 
ods,  may  be  used  as  correlation  references  for  selectively  extracting 
tunnel  reflections  from  any  spatially  complex  structural  responses. 

(5)  Field  Survey  System  Specifications;  The  hardware  components  needed 
to  implement  the  high-resolution  SH-wave  reflection  seismic  tech¬ 
nique  were  defined  based  upon  the  general  findings  of  the  model 
studies  and  analyses.  These  components  can  be  obtained  by  adapting 
certain  commercial  equipment  and  prototype  developments  to  the  indi¬ 
cated  needs.  For  example,  a  48-channel  digital  data  recording  sys¬ 
tem  capable  of  providing  a  4,000  sample/sec  sampling  rate  and  15-bit 
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digital  resolution  is  specified  as  the  basic  field  survey  system. 
Then,  by  adding  special-purpose  components  such  as  borehole-coupled 
detectors  capable  of  providing  three-component  seismic  measurements 
and  sharp  cut-off  low-cut  filtering  in  each  channel,  the  recording 
system  can  be  made  to  respond  primarily  to  high-resolution  SH  waves 
in  the  frequency  range  of  400-1,600  Hz. 

Another  specialized  component  part  of  this  system  is  a  borehole  SH- 
wave  source  capable  of  generating  horizontally-polarized  pulse  sig¬ 
nals  covering  the  required  frequency  range.  This  requirement  can  be 
met  by  an  A  prototype  arc  discharge  transducer  recently  developed 
for  another  borehole  seismic  application  using  compressional  waves 
has  now  been  demonstrated  to  be  capable  of  producing  the  desired  SH 
waves.  This  component,  because  of  its  important  features  of  preci¬ 
sion  pulse  waveform  timing,  pul se-to-pul se  waveform  repeatability, 
broadband  seismic  frequency  spectrum  up  to  about  2,000  Hz,  respect¬ 
able  energy  level  of  1,000  joules/pulse,  and  efficient  operation  as 
either  a  compressional -wave  or  shear-wave  pulse  source,  make  it  an 
ideal  device  for  this  application.  The  additional  system  components 
needed  to  make  the  SH-wave  reflection  seismic  survey  system  complete 
and  fully  self-contained  were  also  specified,  including  a  customized 
seismic  field  cable,  a  programmable  roll-along  switch  capable  of 
handling  either  one-,  two-,  or  three-component  seismic  detector  sta¬ 
tions,  and  a  light-duty  borehole  drill. 

(6)  Recommendations ;  The  results  of  this  theoretical  and  numerical 

model  study  of  SH-wave  detection  of  tunnel  cavities,  suggests  sev¬ 
eral  recommendations  for  continuing  the  studies  and  applying  the 
results.  These  recommendations  are: 

(a)  Implement  the  specified  SH-wave  seismic  exploration  system  in 
prototype  form  for  purposes  of  experimental  deployment,  evalu¬ 
ation,  and  further  development; 

(b)  Develop  and  demonstrate  field  application  procedures  and  data 
processing  methods  which  are  capable  of  acquiring  high-quality 
seismic  survey  data  and  extracting  tunnel  target  responses 
under  practical  field  operating  conditions; 

(c)  Extend  the  use  of  the  SH-wave  numerical  models  to  other  study 
cases  which  may  be  more  accurately  and  more  completely  repre¬ 
sentative  of  the  field  conditions  encountered  in  the  tunnel 
search  and  detection  environment  of  Korea; 

(d)  Conduct  advanced  development  model  studies  in  which  compres¬ 
sional  waves  and  vertically-polarized  waves  are  incorporated  in 
the  analyses.  These  models  will  allow  the  complete  responses 
of  tunnels  to  be  studied,  including  the  conversion  and  inter¬ 
change  of  the  seismic  wave  types  between  compressional  and 
shear  waves,  the  characterization  of  seismic  waves  reflected 
from  tunnels  which  are  oriented  obliquely  with  respect  to  the 
seismic  survey  traverse,  and  the  use  of  these  numerical  models 
to  produce  target-matched  responses  for  possible  use  in  spe¬ 
cialized  tunnel  detection  data  processing. 
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